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The fact that the observed motion of the Moon is apparently affected by inequal- 
ities of long period not yet satisfactorily accounted for is well fiiown to all who have 
interested themselves in the subject. The study of these inequalities was undertaken 
by the author about the year 1870. The discussion of such observations as seemed 
best fitted to throw light on the subject was made at the U. S. Naval Observatory, and 
published in the year 1878, as an appendix to the Washington Observations for 1875. 
The result of that discussion was that the inequalities in question could not be repre- 
sented by any existing theory, and might, therefore, be attributed to changes in the 
Earth's axial rotation, if it could be shown by other independent observations that 
such changes really existed. In the absence of such independent proof an attempt 
was made to represent the outstanding differences between theory and observation by 
an empirical formula; and the curious result was reached that by a quite silnple 
change in what is now known as Hansen's first inequality of long period a fairly good 
agreement with observation could be brought about. 

The changes indicated by this formula were made in the mean longitude of Han- 
sen's lunar tables, and the observations up to the year 1 863 show that the corrected 
motion agrees remarkably well with the fact. The simple way in -which this agree- 
ment was brought about would seem to point very strongly to the theoretical correct- 
ness of the proposed change. Yet the most rigid examination into the subject has 
failed to show the possibility of any eiTor in Hansen's result. 

Meanwhile an independent test of the constancy of the Earth's axis of rotation 
has been made through the discussion of the observed transits of Mercary, published in 
the first volume of the present series. The results of this test were inconclusive, in 
that, while it seemed proved very clearly that the apparent inequalities in the motion 
of the Moon were not wholly due to changes in the Earth's rotation, yet such changes 
were indicated in a somewhat inconclusive way, having a certain general analogy to 
those previously inferred from the motion of the Moon, though much smaller in amount, 
and therefore not very certain. What is most remarkable in this connection is that a 
study of the motions of the Moon, Mercury, and the first satellite of Jupiter all agree 
in indicating that, during some years previous to 1863, there was an increase in the 
speed of the rotation of the Earth which, however, rapidly subsided about that epoch. 
Added importance was thus given to the theoretical side of the question. 

Among the problems of celestial mechanics which admit of rigorous treatment, it 

is well understood that that of the action of the planets on the Moon is the most diffi- 
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cult which presents itself. It was therefore necessary, in commencing the work, to 
devise a method by which it was supposed these difficulties could best be surmounted, 
and especially, by which every possible inequality could be detected and investigated. 
It occurred to the author that a more rigorous and general treatment than that hitherto 
adopted would be to consider the arbitrary constants which enter into the solution of 
the problem of tlu'ee bodies, the Sun, Moon, and Earth, as simultaneously varied by 
the action of a fourth body, a planet, using La Grange's method of the variation of 
elements. The determination of the La Grangian coefficients necessary in the theory 
proved to be a work of more simplicity and elegance than would have been supposed. 
Accordingly he took advantage of a visit to Paris in Api-il, 1871, to read, before the 
French Academy, a preliminary chapter on the subject, which was soon afterward 
published in LiouviUe^s Journal de Mathematiques Pures et Appliquees. In the course of 
the next two years the results of Delaunay's great work were utilized for the purpose 
of forming the actual values of the expressions which entered into the theory, and the 
greater part of the numerical developments pertaining to the action of the individual 
planets was also effected. The computation of the inequalities of long period produced 
by the action of Venus and Mars was also supposed to be nearly completed. 

At this stage of the work a more careful investigation was made of its foundation, 
leading to the discovery that the author had wrongly measured the importance of 
quantities in the theory of three bodies by the numerical value of the perturba- 
tions to which they gave rise. The motion of the center of gravity of the Earth and 
Moon around the Sun, as determined by the action of the Sun itself, deviates from an 
invariable ellipse by quantities of no*practical importance in ordinary astronomy. It 
was, therefore, supposed at the outset that the elements of the orbit of the Moon need 
not be regarded as contained in this motion. An examination of the problem showed, 
however, that this conclusion was ill-founded, and that, so far as mere order of mag- 
nitude went, the consideration of the lunar elements in this motion was as necessary 
as in the motion of the Moon itself It therefore seemed to be necessary to form expres- 
sions for the solar elements in terms of the lunar elements, a work which, if completely 
executed, might be very laborious. Under these circumstances the work was laid aside 
until the author should find the leisure necessary to carry through the required inves- 
tigations. Attempts were made from time to time to overcome or evade the difficult}^, 
but none of them were brought to an entirely satisfactory conclusion. Indeed, the 
author is quite prepared for the conclusion that the rule found to hold in the theories 
of the planets, that the more general and abstractly rigorous a method of considering 
perturbations is, the less it is adapted to their actual numerical computation, may hold 
true in the present case. He has, however, at length deemed it best to publish the 
work, notwithstanding its lack of completeness, for the following reasons: 

I . Its analogy with some of the modem theories of celestial mechanics, and its 
possible value as illustrating those theories, especially in the form developed by 
P01NCAR6 in his recent treatise Les methodes nouv^lles de la Mecanique celeste. The path 
on which it proceeds to extend the solution of the problem of three bodies to that of four 
may perhaps lead to the method by which the most general solution of the equations 
of planetary motion is hereafter to be investigated. 
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2. The numerical developments for the action of the planets, from Mercury to 
Jupiter, inclusive, are equally necessary in applying any other method. Indeed, one 
circumstance which led to the question of publication being considered at the present 
moment, was the remark that these numerical developments might perhaps have been 
of use to Radau in his recent elaborate researches on the subject. 

3. Many of the analytic developments of the present method will be found closely 
analogous to, or identical with, those of other methods, and may thus serve as a control 
on the latter, or as a basis for their immediate application. 

4. We have also to consider the usefulness of even an abstract method of consid- 
ering the theory in its most general form, in that a more luminous view of the scope 
and principles of other methods may be obtained. 

5. Finally, the author is not without hope that some way of avoiding or over- 
coming the difficulties he has encountered will be devised by an abler investigator; or, 
at least, that such an investigator will find the work useful. In any case the problem 
of applying the method is now merely that of completing the analytic theory of the 
motion of three bodies. 

The first part of the work, relating to the general dynamical theory of the subject, 
was reconstructed at various times between the years 1873 ^^^ 1888, with a view of 
putting it in such a shape as would best admit of farther development and wider appli- 
cation. This remark applies especially to §§ 2 to 4, which, though first in order, are 
practically the last that have been prepared. All the rest of the work is now published 
as originally written in the years 1871-1873, except ordinary revision of details. 
Being once shaped into a logical and continirous whole, it was scarcely possible to 
make any extended changes or additions, except by a process which would have been 
equivalent to a complete reconstruction. The author was the more willing to refrain 
from such a course because he might then be more readily excused if it should appear 
that he had not duly considered the advances made within ttie last twenty years in the 
theory of the subject By this course he has made plain an oversight in his conclusions 
as to the possibility of inequalities being produced by Jupiter, whereby what is now 
known as the Jovian evection was left to be discovered by Mr. Nevill. 

It may not be amiss in this connection to call attention to the fact that the problem 
has lost none of its importance during the time which has elapsed. There are three 
possible sources to which may be attributed the apparent inequalities in question: 

(a) Defects in the gravitational theory. 

(>ff) Inequalities in the Earth's time of rotation. 

(y) Action on the Moon of other forces than the gravitation of known bodies. 

The second of these causes seems to be not at all improbable, in view of recent 
discoveries respecting variations of latitude, which can be accounted for only by minute 
changes in the Earth's axis of rotation. But the independent tests of this cause offered 
by the observed transits of Mercury do not show that it suffices to explain the observed 
inequalities in the motion of the Moon. Two other tests remain to be applied: the 
motion of the first satellite of Jupiter and observations of Mercury, Venus, and Mars. 
So far as evidence on these points has been brought out, the conclusion indicated is 
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that inequalities of the Earth's rotation are real, but not large enough to produce the 
observed effects on the Moon's mean longitude. The subject is therefore still involved 
in doubt, and the first step toward the removal of this doubt is such a determination 
of the inequalities in the mean motion produced by the action of the planets, as 
shall place the subject on an undoubted basis. 



GENERAL THEORY. 



§1. 

FUNDAMENTAL DIFFERENTIAL EQUATIONS FOR THE DISTURBED AND 
UNDISTURBED MOTION OF THREE BODIES. 

For convenient reference, a part of the notation used throughout this paper is 
gi^Kjn here. We put 

mi, W2, W3, m^j the masses of the Sun, Earth, Moon, and planet, respectively; 

X, y, z, the rectangular coordinates of the Moon relative to the Earth; 

a/, 1/', y, the rectangular coordinates of the Sun relative to the center of gravity 
of the Earth and Moon; 

M = ' — 

^W2 + ^ 



^1 + ^2 + ^*3 

/ 



X4, Y4, Z4, etc., the rectangular coordinates of the planet relativel}' to the center of 
gravity of the Sun, Earth, and Moon; 

x^j ^4, ^4, r4, the rectangular coordinates and radius vector of the planet referred to 
the Sun; 

r, the radius vector of the Moon, so that 

r^=rr^+2^+^^; 

r', the Sun's radius vector from center of gravity of Earth and Moon; . 

p, the distance of the planet from the center of gravity of the Eartli and 
Moon; 

a, e, ;r, etc., the elements of the Moon's orbit round the Earth, the major axes 
being determined by the condition a^t^^=sum of masses; 
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a', 6', 7r\ etc., the elements of the relative orbit of the Sun round the center of 
gravity of the Earth and Moon ; 
w, w', the mean motions of the Moon and Sun, respectively, as actually 
observed; 

/i, the potential, or sum of the products of each pair of masses divided by 
the distance, when only three bodies are considered; 

R, the terms of the potential added by the planet; 

Pzz/i+R, the entire potential; 

^i> Vij -2^11 ^'\i 6t^-y t^® derivatives of x^ y, z^ x\ etc., as to the time. 

If we suppose, for the present, ^<, 77^, and ^< (i having in succession the values 
I, 2, 3, 4), to represent the rectangular coordinates of the four bodies referred to any 
arbitrary fixed origin, and ^o> ^o? and ^0 to represent the coordinates of the center of 
gravity of tlie Sun, Earth and Moon, the differential equations of the motion of the three 
bodies will be 



nii 



-=.f^^^^' i^O 



de ~d^,' * de ~drfi 

Both sides of these equations admit of being completely expressed in terms of the 
three relative coordinates already defined. The expressions for the relative coordinates 
X, y^ZjX^y^j z' in terms of the absolute coordinates are 

X ZZ ^3"~" c?2 

(2) 



1 =^3-?2 



a;'=^i— (i— )u)^2— // ^3 

y = '7i_— ( I — /<) % — /^ '/s 

DiflFerentiating theae equations twice, and substituting the values of the second 
derivatives iu (i) we find 

rf^_ i^[P_ I dV 
dfi t»3 d ^3 W2 d ^2 

(3) 
^a/_i^rfP I {dP dV\ 

d fi ~«ii d^i Wj-I- W3 Kd H2 d f.3) 



with similar equations in iffi/,^ and /. 
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The four absolute coordinates ^i, ^2> ^8> ^a admit of being completely expressed as 
a function of the four quantities x^ x\ ^o, and X4. If, therefore, we suppose P expressed 
in terms of the last four quantities, we may put, in (3) 

d^,~dxd^,dx' ?^idX,d^,d^odSi ^^^ 

the last term, however, being identically zero, because ^0 disappears from P. 
The complete expressions for the coordinates ^ are 

Representing the distance of the bodies numbered i and j by p<j, the complete 
expression for P is 

the general value of p<j being derived from 

The expressions for the six values of Pij in terms of the relative coordinates are 
from (5) 

p\s=(a;'-(i-//)a;y+(y-(i-Ai)yy+(^-(i-/^)0' 

p\,=(X,-M' x7+(Y«-/i' y'y-^.(Z,-M' z'y 

p^3=r^=^+/+^ (6) 

+(z*-\-(i-M')z'-(i-M)^y 

From (2), combined with the value of X4 derived from the expression (5), namely, 

^ _^ _ ^h ^1 + ^2 ^2+^3 ^3 

we find the partial derivatives 





dx 


dx _ d X _ 
dS~ d^3~ 






d3f 

dar 


dx' dx' 

di;,-^ dar 




dX^ 


w, 


dXi_ m^ dX^__ 

d^2~ »»! + »»2 + W»3 d^3~ 


rth 


dHx~ 


Wi + Wj+Wg 


»»1 + Wl! + W3 
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These values being substituted in the general expression (4), we find for the 
values of the derivatives in (3) 

d^i~dx' mi+wi2+wi3rfX4 

dV dV , , vrfP w. rfP 



d S2 dx dx' nii+m^+rn^dlL^, 

dV_dV_ rfP_ ^?3 drP^ 

dSz dx dx- Wi+^W2+^^X4 

By the substitution of these values and the corresponding ones for the other two 
coordinates equations (3) become 

d'^_dV 
^""dt^^dx 

^_dV 
d <2 dy 



d^'z dV 

(7) 



'^df-dz 



d^x'_dF 
^'dfdx' 

(Py'_dV_ 
^'df-'dy' 

d^z'^dF 
^'dt'-'dz' 

where P is to be cx)nceived as expressed in the form 

P=/ (x, y, z, x\ y\ z', X4, Y4, Z4) 
by the substitution (5). 

These are the fundamental differential equations of the problem. As thus 
presented, the coordinates of the planet X4, Y4, and Z4, are supposed to be given in 
tenns of the time. The differential equations for these coordinates, supposing the 
planet to be acted upon only by the Sun, Earth, and Mo©n, would be of the same 
form, and if joined to (7), the whole would represent the differential equations of the 
problem of four bodies. In reality, however, we should add nothing to the astronom- 
ical completeness of our problem by considering it as a problem of four bodies, because, 
if we regard the coordinates of the i)lanet as depending on the attraction of the Earth 
and Moon, we must equally regard them as depending on the action of the other 
planets. We shall therefore regard our problem as that of the disturbed motion of 
three bodies, the position of the disturbing body being supposed to be given inde- 
pendently of the integrations by whicli the relative coordinates of the tlu-ee bodies are 
determined. 
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§2. 

EEDUCED FORM OF THE POTENTIAL. 
Let us first consider the principal portion /2 of the potential P, namely, 

Pi, 2 Pi, 3 P2, 3 

the three denominators having the values given in equations (6). The coordinates x^ 
y, and z being very small in comparison with af^ i/j and /, the first two terms of /i may 
be developed in powers of their ratios. The first approximation will be the value 
which these terms assume when we suppose Xj y, and js to vanish. Putting 

r^ =^ +/ +^=P\3, 
we find for a portion of ifit, which we shall call D.q, 

fl^—'^^h (^>2+^%) _|_ rth ^3 /g\ 

while the remaining terms of £1 will be 



/2,=.^.{m,(-L_^,)+„,,(l_^,)| 



Now let us put 



W, the aiagular distance of the Sun and Moon, as seen from the center of gravity of 
the Earth and Moon. We shall then have from (6) 

p\, = /^jl+2/l^^COsW+/i^^,| 

p\3=r-{i-2(i->u)^cosW+(i-/i/^,} 

Developing the power of — - of these expressions to terms of the second order in 

-. we find 
r 

-i=;|i-/.^cosW+V^*:',(3cos='W-i)} 
Substituting these values in /ii, and noticing that by the expression for /i, namely, 

OT2 + »t3 
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we have 

tn3/.2+m3(i-/.)^=-^^=>u2 
mj + ma 

we find 

^ inwfhr^ (3 008=^ W- 1) (9) 

The expressions (8) and (9) make up the complete value of *Q when terms of the 
third and higher orders with respect to — are omitted. 

The remaining terms of P are those depending on the mass of the disturbing 
planet, the expression being 

V1.4 ^2,4 Pzy 

For convenience in subsequent developments we shall substitute for the first tenn 
the integral of its derivatives with respect to the variables x\ y\ ^ ^ which we do because 
only these derivatives enter into the fundamental differential equations, and p^^^ does 
not contain any of the quantities x^ y, or z. This term of R may then be replaced by 
the expression 




d — d — d — 



Substituting the value of Pi, 4 in (6), and putting after differentiation 

^l Va, z^i U 

for the coordinates of the planet relative to the Sun, or 

0^4 = X4 — /// a;', etc., 
this integral becomes 

^ P 1,4 

and may be expressed in the form 

^,m,/.-^^^+y/+"«^ (10) 

4 

provided that in forming the partial derivatives we regard the coordinates of the planet 
relative to the Sun as independent of x', y\ and /. 

The remaining two terms of R may be developed in the same way that we have 
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just developed /ii. If we put 

p, the distance of the planet from the center of gravity of the Earth and Moon, 
and for brevity, 

5, Tfy ^, the coordinates of the planet relative to the same origin, we have 

Taking the powei' — and putting, for brevity, 
we find for the principal terms 

Multiplying the first of these by w^ Wa, and the second by m^ m-^, and adding, we 
find for the corresponding terms of R 

M p 2p^^ 2p^ ^ ) ^ ^ 

The first term is subject to the inconvenience that p is there to be regai'ded in 
differentiating with respect to jd as a function of the coordinates of the Earth and planet 
relative to the center of gravity of the Sun, Earth, and Moon, and therefore as a func- 
tion of a;'. To avoid this inconvenience, and introduce the coordinates of the planet 
relative to the Sun, we shall make the same transformation we made on the first term 

of R, putting instead of - the expression 



f d- d- A- \ ' r\-ij . 



When we introduce instead of X4, Y4, Z4 the coordinates of the planet relatively to 
the Sun we have 

X4=a?4— (i — //') a/, etc., 
and hence 

5=:C4+a/, 7=^4+!/, ?=^4+^ • 

Substituting these values in the integral it becomes 
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SO that the first term of (i i) is simply to be multiplied by i — /<', and p to be regarded 
as a function of the heliocentric position of the planet only. 

Collecting the various terms of R in (lo) and (ii) as modified by (12) we have 
its complete transformed expression. The fundamental differential equations (7) now 
become 

dy _dn dR 
^'dt'^'-dx'^daf 

dW_dn dB. 
^'df~dy''^dy' 



d^_dn dR 
dfi d/ dsi 



d'^x _dD.,dR 
^""dt'-dx^dx 



(13) 



d^y _dn,dB. 
d^ ~ dy dy 



d^z_dn dR 
^^de~ dz'^ dz 



Where 

/2 =/io+/2i 



d^o— -? 1 — - 



/ 



(H) 



/2i=-mi//2^3(^3 cos^W— ij 
in which we liave put 

W2+W8 

It will be remarked that the variable part of f^i is identical with the force function 
for the action of the Sun in the lunar theory, and that R is closely analogous to the 
disturbing function for the action of a planet in Delaunay's theory* It is also to be 
noted that Hj and R are developed in powers of the ratio of the mean distance of the 
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Moon to that of the Sun or planet, and that only the first terms of the development are 
included in the expression (14). 



)^ 



§3. 



COMPARISON OF TWO METHODS OF INTEGRATING THE PRECEDING EQUA- 
TIONS. 

The differential equations (13) of the preceding chapter, or some transformation 
of them, form the necessary starting point of all methods of treating the theory of the 
Moon's motion. It may also be held that the most rigorous method hitherto developed 
of treating the problem analytically, and now well known, is that resting upon the 
results of Delaunay's system of integrating the equations in question, or upon some 
process equivalent to it. 

The problem as actually considered by Delaunay was that of three bodies only, 
considered as material points. Even with this restriction he did not fully complete the 
solution, because he considered the relative elliptic orbit of the Sun to be performed 
around the center of the Earth, instead of the center of gravity of the Earth and Moon. 
So far as the Moon's motion is concerned, the additions required to reduce his solution 
to the actual case are extremely simple, so that for all practical purposes we may con- 
sider Delaunay's theory to be complete as regards the Moon. The deviations of the 
path described by the center of gravity of the Earth and Moon from an elliptic orbit 
around the Sun are, from a theoretical point of view, necessary to complete the solution, 
but are too minute to be of importance in observational astronomy. 

To complete the solution of the problem of the Moon's motion subsidiary actions 
have to be taken into account, such as those arising from the ellipticity of the Earth and 
from the action of the planets. It is to the two methods of effecting these subsidiary 
solutions that the attention of the reader is now called. One of these was doubtless 
that intended by Delaunay himself, and consists essentially in continuing the method 
adopted by Delaunay, considering the actions in question simply as perturbations 
additional to those produced by the Sun. This method has been followed by G. W. 
Hill on Delaunay's lines in Vol. Ill of these papers. The methods of other investiga- 
tors, as Radau and Nevill, may also be considered as founded substantially upon the 
same idea. This I shall designate hereafter as the first method. 

In the second method, which it is the object of the present paper to develop, all 
the constants of integration of the problem of three bodies are considered as simulta- 
neously variable under the action of a fourth body, a planet. In the actual applica- 
tion of this second method there is so close an analogy to the first that the two must 
be mutually helpful to each other. Their relations can be clearly shown by a compara- 
tive statement of the general plan followed in them. 

In both methods the solution of the last three equations (13) is obtained by three 
successive steps. In the first step, which is common to both methods, we have, instead 
of these equations, those which result from putting il^ for il-^-R^ the value of the force 
function then used being 

A. p. — ^VOL. V, PT. Ill 2 
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The results of this step are the theories of eUiptic motion; the moon revolves around 
the Earth in an ellipse, as does the center of gravity of the Earth and Moon around 
the Sun. Thus in the solution the Moon's coordinates, x, t/, and ^, are represented as 
functions of six arbitrary constants, say 

^1,02, . . . ae 

and of the time; while the Sun's coordinates, rz/, }/, and /, appear as functions of six 
other arbitrary constants, say 

a'l, a'aj . . . a'o 
and of the time. 

In this step the effect of the solution is that the twelve arbitrary constants of 
motion are completely separated into two classes, six appearing only in the expressions 
for the Moon's coordinates, and the other six only in those of the Sun. 

In the second step /ioH-*^! i» substituted for /i^ in the second members of the differ- 
ential equations which determine the Moon's motion, /ii is here a function of the 
coordinates of both bodies, and the value of its principal term is given in (14). The 
complete solution of this step is effected by regarding the twelve elements as variable. 
Since at the first step the Moon's coordinates contain on)y unaccented elements, which 
we may call lunar elements, and the Sun only the accented elements, which we may 
call solar elements, it follows that the I^a Grangian coefficients which pertain to the 
Moon will contain only the six lunar elements, and those which pertain to the Sun 
only the solar ones. The solution might then be made by successive approximations, 
the lunar elements being first regarded as constant in integrating the equations for the 
Sun's motion, and the solar elements as constant in integi'ating the equations for the 
Moon's motion. The values of the twelve elements obtained in this first approxima- 
tion may then be substituted in the second members of all the equations, a second 
approximation made, and so on until a solution of tlie requisite accuracy is obtained. 

In Delaunay's theory, as in all others hitherto published, the attention being con- 
fined to the Moon's coordinates, the solution is effected only for the lunar elements, and 
moreover the solar elements are regarded as constant, it being practically unnecessary 
to consider the variations of the solar elements due to Oi. Successive approximations 
are, however, made in the case of the lunar elements. These elements, and hence the 
lunar coordinates, thus become functions of twelve arbitrary constants, six of which 
pertain to the Sun alone, and are not constants of integration, but are considered as 
entering into the differential equations of the Moon's motion. The other six are the 
constants of integration, and may be called lunar constants. We may represent these 
six lunar constants by the symbols 

«i, «2 • . . . «e 

So far we have only the solution of the problem of three bodies, which, in the second 
method, is considered to be solved in advance. The divergence of the two methods 
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now commences in the continuation of the sohition. As a matter of fact, the solar 
elements, say 

a'l, a'2 • • • • ^e 

become variable through the action of the planet. The determination of the variations 
oflFers no difficulty, their value being the same whichever method be adopted. But in 
the first method, when the coordinates of the Moon are considered, these solar elements 
are not regarded as constants of integration, but as quantities given a priori^ and appear- 
ing in the differential equations. The La Grangian coefficients 

*' ^ da^daj duida^ da^daj 

are therefore formed only for the combinations of the six lunar elements, and the per- 
turbative function R appears in two parts, the one being the value of R, as found 
in (14), the other the value of S fl^ arising from the perturbations of the solar elements 
through the action of the planet, which we may conceive expressed in the form 

Thus in .the final solution obtained in this way the Moon's coordinates will not 
appear as functions of the actual variable elements of the Sun, but of mean elements 
at some epoch, the entire effect of the perturbations produced by the planet being 
thrown upon the lunar elements. 

In the second method, all twelve elements are regarded as simultaneously variable, 
and the La Grangian coefficients are formed accordingly, so that the lunar coordinates 
are in the first plane considered as functions of all twelve varying elements. 

The final forms of the solution by the two methods may now be seen and com- 
pared. 

In the results of Delaunay's work, where only three bodies are considered, the 
Moon's coordinates are expressed as a function of the twelve constants 

«i, ^2 • . . • «6, «'i, a'2 . • • • a'c 

When the action of the planet is included by Hill's continuation of Delaunay's 
process, the Moon's coordinates are regarded as still the same functions of these twelve 
elements, the accented ones remaining constant, and the others becoming variable. 

In the new method the coordinates are also regarded as still the same functions of 
these elements, but all twelve of the elements are regarded as variable. 

We may now construct the equations of condition, which must be satisfied in order 
that the two methods may give the same result. Let us represent by the symbol <S' 
the perturbations of the elements by the first method, and by S those by the second. 
The corresponding perturbations of any coordinate v will then be 
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duy^ da^ da^ 

TT i^ dv J. , t dv J. . dv J. , . I ^^ jt ^ 

da^ aa^ aai doc^ 

The equality of these perturbations gives the relation 

da,y / V da\ * (15) 

(iz=i,2, .... 6) 

which is the required relation. 

§4. 
FUNDAMENTAL EQUATIONS OF THE" METHOD. 

The present paper being devoted wholly to the development of the second of the 
preceding methods, we have to consider further the general method of proceeding. 
We recall that in the general La Gbangian method, which alone will be used, it is pre- 
supposed that a solution is obtained giving the values of the coordinates in terms of the 
requisite number of arbitrary quantities, such as shall completely satisfy the diflferential 
equations when certain terms of the potential are omitted. A solution of the case when 
the omitted terms are taken account of is tlien obtained by a system of ordinary differ- 
ential equations, in which the variables are the arbitrary elements introduced by the 
preceding solution. There is no limit to the number of times by which this process 
may be successively repeated. The effect of each new term of the potential introduced 
is taken account of by the integration of a system of differential equations in which the 
arbitrary constants of the preceding solution appear as the variables. Delaunay's 
classic work forms a beautiful example of the operations by which the equations of the 
Moon's motion, as affected by the action of the Sun, may thus be integrated. 

The problem with which the present paper is ultimately concerned is not that of 
the motion of the Moon as affected by the action of the Sun, but as affected by the 
action of the planets. This action is involved only in the term R of the potential. 
Hence the sole object of the present paper is to to determine the variations of the lunar, 
and incidentally of the solar elements produced by the term in question, the equations 
being completely integrated for the case when R is omitted. The general nature of 
the process is so well known that no development is necessary. We shall therefore 
merely present a r(^sum(5 of the necessary equations. 

Firstly. We presuppose that the values of the six coordinates x, y, z, x\y\z\ are 
expressed in terms of the time and of twelve arbitrary constants, which we may rep- 
resent in the most general way by the symbols 

«!, a2> • • • • ^12> 

in such way as to completely satisfy the six equations 
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dt^~dx ^'dt^~~dx' 



d^y dn d^i/ d£l , ^. 

^z _dn dP/_dn 

^^dt'-dz ^'dt'~dY 

Let U8 represent these expressions in the form 

x=f, y=Fi 

y=A y'=F, (17) 

z=zA ^'=F3 

where / and F are functions of the twelve arbitrary constants and of t. 

Secondly. We construct the conjugate functions r„ //j, ^„ x'l, 1/1, z\ by the equa- 
tions 

dfi , dY^ 

x,=^,/; x,=Mr ^/ 

rf/3 / d F, 

The variations of the twelve elements when R is taken into account must then satisfy 
a simultaneous system of twelve differential equations, of which four are as follows: 

dx day d X d a^, , d Xi2d a^z _ 

da^dt da^di • • • • ^ ^^^ dt '^ 



(19) 



dxida^ dxid a2 . . dx^ da^^ rfR 

da^ dt da^dt ' ' ' * da^z dt dx 

dx' da^ dx' d a2. , dxf da^^ 

da^ dt da^ dt ' ' ' " da^^ ^' ~ 

dx'idai dx'ida2. dx'i d a^z ^R 

da^ di da2 dt . . • . ^^^^ ^^ ""rfj-' 

while the remaining eight are formed by writing y and z in the place of rr, and need 
not therefore be actually written. 

From these equations are derived, by a well known linear combination, the 
La Grangian equations for the variation of elements in the following form. We put 

r T ^/^dxdxi dxdxi\ z^^x 
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the sign 2 indicating a summation of the temis which arise by writing ir, y, Zjx\y\ / in 
succession for x. 

We also suppose that R becomes a function of the elements a„ and of the time by 
the substitution of the values of the coordinates as functions of tliese quantities. The 
La Grangian equations for the variation of the elements then take the following form: 

r T^^iir 1^*^2 1 ir -t d a^2 ^^ / \ 

[«.-.] ^/+K-.]^/+ +K..]-^-=3- (21) 

(i=i; 2; 3; 12) 

In determining these La Grangian coefficients I make use of the known theorem that 
they are constants. This theorem rests, however, upon the supposition that the expres- 
sions for the coordinates actually satisfy the differential equations. As a matter of fact 
they are developed in series, so that it is impossible that this condition shall be rigor- 
ously fulfilled. It is therefore important to determine finally to what degree of approx- 
imation the constancy of the coefficients in question is assured. For our present purpose 
I have deemed it sufficient to develoj) the constant portions of the coefficients so far as 
the degree of approximation to which the solution is carried would admit. 

This theorem of the constancy of the coefficients enables their values to be formed 
with much greater facility than if it were necessary to actually perform all the multipli- 
cations which enter into the equation (20). 

In considering the general theory a yet farther simplification may be made- It is 
known that if we take six of the elements, which we may call 

fri, r2, . . . . /(s 
at pleasure, we can always substitute for the other six certain functions of them, 

Cij C2 . . . . Tg 

such that the La Grangian coefficients shall reduce to zero in all the combinations 
except those six in which an /< is coordinated with its corresponding r<, in which case 
the combination will be reduced to unity. We shall make frequent use of this theorem 
in the course of the present investigation. 

Such a selection of elements I shall hereafter call canonical. 



§5. 
THE PROBLEM OF THREE BODIES. 

The preliminary step now is to express the six coordinates, x, y, z^ x\ yl ^ /, in terms 
of twelve arbitrary constants, so as to satisfy the differential equations (16). 

We recall that two processes of integration are necessary in the preliminary inves- 
tigation, namely: 

1. We integrate omitting £1^ in the equations (13), writing /2o for /i + R. 

2. We include n^, but omit R, as in the equations (16). 
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The first step involves no difficulty, since it gives rise only to the equations of 
elliptic motion — 

(a) Of the Moon around the Earth ; 

(fc) Of the center of gravity of the Earth and Moon around the Sun. 

Putting in a general way 

Ci, C2j C3, fj, Uj '2 

for a set of canonical elliptic elements pertaining to the Moon, and 

r' f' r' V V V 

for those pertaining to the Sun, we shall then have, as the result of putting ilji^o 

a\ y, and z z=:f (cj, c^, ^3, h k, h ^^^^ 

(22) 
x', y , and / =/ (c'l, r'2, r'a, V,, l\, }\, t) 

Next, to take account of /2„ we substitute in .O^, for ir, y, ^, a;', y\ z\ their expres- 
sions (22). Then, using canonical elements, the form (21) will become 

dt ~ dli dt~ dCi 

(i=i,2,3) (23) 

d(^^_dn, dl\__d£\ 

dt ~dr\ dt ~ dc\ 

These equations are now to be integrated by successive approx;iraations. We have first 
to integrate the equations for c and / on the supposition that c' and l^ are constant. 
Next, by substituting these values of c and I in the equations for c' and Z' and integrat- 
ing, we obtain approximate expressions for these quantities. If necessary these new 
values of & and I' can be employed to reintegrate the first set of equations, and so on. 
The first integration has been performed by Delaunay in his great work Theorie 
du mouvement de la Lune, We may consider that, in this work, Delaunay has implic- 
itly formed expressions c and I which substituted in (22) give values o( x, y, and z that 
completely satisfy the equations 

^'~df- — di — ^^- ^^^ 

on the supposition that c'i and /\- are treated as constants. More exactly, if we put, 
in (14) 

^ ~ r 

(25) 



llZ Q -^ 



(30) 
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§6. 

CONSIDERATIONS ON THE FORM AND ORDER OF MAGNITUDE OF THE 
EXPRESSIONS WHICH ENTER INTO THE PRECEDING THEORY. 

The fifteen coefficients [a», aj which enter into the equations (21) are each divis- 
ible into two parts, one of which is formed by differentiating the lunar coordinates Xy y, 
and ^, the other by differentiating the solar coordinates x\y\ and /. If we put 

*- * ^ d a^ daj d a^ d a^ da^d a^ 

dx dxi dy dyi dz dzi 

da^da^ da^d a< da^da^ 

^ * ^ da^ da^ da^da^ d a^d a^ 

_dx' dxfi_dy' dyfi_dii^ dz'^ 
d a J da^ d a^da^ d a^d a^ 

\da.^J\~ d xd a^ dyda^ dzda^ ^'^ ^ 

fd'R\ _d'Rdx' dRdy' dRd/ . . 

\d aj^^d a/da^'^dydai'^d/dai ^^^^ 

we may form by means of the expression (17) the two equations 

^^'r -. da, /rfR 



(30 



^=^\ ■ da. ^dK\ , . 

f^\ ^daj ^dR\ . . 



and the sum of these two equations will make up the equation (21). 

We may divide the twelve elements a into two classes, which may be designated 
lunar and solar elements. The first class comprise those which, when the Sun's dis- 
turbing force vanishes, merge into the elements of the Moon's elliptic motion around 
the Earth; the second, those which, when the mass of the Moon vanishes, merge into 
the elements of the Earth's elliptic motion around the Sun. 

Let us niBxt investigate the order of magnitude of these several coefficients and of 

(^ — j and ("j— ) • Ii^ tl^^ expressions (17) for the solar coordinates we are to dis- 
tinguish two portions, one a purely elliptic part containing only the six elements of 
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the elliptic motion of the Sun, the other part containing the terms added when we take 

into account the distance of the Moon as compared with that of the Sun. The first 

part arises from the integration of the differential equations for rr', y', / when only the 

part *f2o of the potential is considered, and is equal to the mean distance of the Sun, a', 

multiplied by purely numerical ratios. 

In the expressions for x, y, and z, all the terms containing solar elements will be 

nf 
multiplied by — . 

We shall use the symbols S, S', etc., to represent in the most general way the sum 
of a series of sines and cosines of angles, each multiplied into purely numerical ratios. 
These symbols will not twice represent tlie same quantity, and may, in special cases, 
reduce to purely numerical constants. If we put for brevity 



n' 

mzzz — 

n 



7 «' 



(36) 



the expressions of the form (17) and (18) for x^ x^^ etc., and their derivatives may be 
put in the form 

x =a (S+mS') 

:r^ =a' (S'+AS) 

Xy :zi//2«w(S+wS) 

etc. etc. 

If, in forming the derivatives of these expressions as to the elements a^ , , , ajg, 
we choose for the latter purely numerical quantities, which we may always do by 
replacing a line element by its logarithrti, then the derivatives in question will be of 
the same form and the same order of magnitude with x, x\ etc., themselves, or more 
exactly with the terms of those quantities which contain the elements. If, then, as 
before we distinguish the solar elements by accents, the order of magnitude of the 
largest terms of several derivatives will be as follow: 

a a 

Ct X\ c% 

-y-^ = /Yg a ^ S 
a a 

(38) 

da 
a a 
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-^-7 =maS 
da' 

d Xt cii 

da 

(38) 

d xl i Q, 

a a 

The various classes of coefficients [a<, aj formed by (30) and (31) will therefore 
be of the following orders of magnitude: 

[a, ai]i=:>U2«^wS 

\a^ a'^^-zimma^n^ 
\a\ a' J, = //j m^ a^ w S 

(39) 

[a, a'0,=/^,A;a'^«'S 

Let us next revert to the perturbative function R as found in (14). We see that 
it is composed of two parts of very diflFerent magnitude. Considering all linear helio- 
centric coordinates as of the same order of magnitude, a', and the coordinates of the 
Moon as of the order of a, the order of magnitude of the two terms of R will be 

f 

73 



p^_/W4jii j^ ;W4 /i2^2 g 



or, if we use coordinates instead of mean distances, the orders of magnitude will be 
represented by 

where the first term does not contain the Moon's coordinates. 

From this equation we shall have, using the sign = only to indicate identity of form 
and order of magnitude, 

rfR m^M2^(;i W4//2«G 
dx X' a^ 



PRODUCED BY THE ACTION OF THE PLANETS. I 25 

and thence, from (32) and (33), having reference to the orders of magnitude of -j^, etc., 

dot 

\daJi~ ~(]P ~ 'cT' 

rd^\ _m^fjLik^ 
\da )%~ a' 

(40) 

\.do^J\~ a'^ ~ a' 

\daO2- a' 
The approximate values of >Ui and yug are 

/il = W2 

/^2 = W3 

by which we shall replace them in R. The relations, between the mean motions and 
mean distances give, approximately, 

or n zn — - 
an 

a n 

Substituting these values in the expressions (39) for the coefficients [^i, aj, and 
representing by Si, Sg, etc., Ui, U2, etc., numerical quantities like those we have hitherto 
represented by the general symbol S, we have 



an 



an 



an 

K, <.]i = -^m2Ul^> 
an 

a n 
a n 
a n 



(t=rl,2 ... 6) 



(41) • 



a n 
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The expressions thus formed are to be substituted in the pair of general equations 
(34) and (35) which, it must be remembered, may belong to either of two classes, accord- 
ing as a^ there represents a solar (accented) element, or a lunar (unaccented) one. Let 
us consider the former case first. The right-hand members of these equations, taken 
from (40), may then be put in the form 

( —-J ) iz: ni -A_* ]c \](^>—^l—^M -7-, U^** 

(ll).=T*(S'»+*S") 
The pair of equations then become 

an I at at at V at aty] 

(Z 

(42) 

an [ \ at at y \ at at y ] 



n,3 



a a^ 

The sum of these two equations gives one of the rigorous La Gbangian equations 
for the variations of the elements, in which each coefficient of the differential of an 

element is a constant. The coefficient of —=-^ in the sum is 

at 

m,m,g,,, m^tHs^^, Ur,^[< a',]=constant. 
an an l 7 j 

Here the coefficient of S^ is greater than that of U< in the ratio — -f- -^^o, or about a 

1^3 an 

hundred millions to imity. The constant terms of Si^^ must therefore be to the periodic 

terms in at least this ratio. The coefficient of Ui*^ may therefore be entirely neglected 

along side that of S|^^ 

Again, the approximate value of the coefficient k is 



Wja'^ 12,000,000 



the terms multiplied by this quantity may therefore be entirely neglected. The varia- 
tions of the solar elements may therefore be determined entirely from the first of the 
pair of equations by a system of six equations of the form 
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Sf)^+S?>^H . • . +Si«^«=^n'S<« (43) 

at at at mi ^ 

which will, in fact, be the usual equations obtained on the supposition that the masses 
of the Earth and Sim are combined in their center of gravity. This hypothesis is 
known to suffice even for the action of the Sun, a fortiori must it suffice for the disturb- 
ing action of a planet. 

Let us next consider that in the pair of equations (34) and (35), «j represents a 
lunar element. If we then express the first two equations (40) in the form 

\daJi a' 

(44) 

f^^\ _ ^2^4 J^ g(7) 

\daJ2 a' 

if, also, we substitute for the coefficients their values in (41) the pair of equations (34) 
and (35) assume the form 

an \ at at - y an \ at at J 

a' 

^^2;i^rS{*>^^^+S^*>^^^+ . . V^,^A:rSf>^'HSr^,^^^ . . ^ 
a n \ ' at at y a n ^ at at / 



Mzni^ 



a' 



The sum of these two equations forms the final equation for the vaiiations of the lunar 
elements after the variations of the solar elements are determined from (43). If we 
multiply their terms by a n-i-^Wamg, putting for k its value 



JTf2 



Wg a* 



ni2+nt^a^^ 

using the relations 

a^ n^ zumg+mg 

a'^ n'^ zn Wi + ^^2 + Wb 

and regarding as unity the ratio 

the equations become 
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the sign 2 indicating a summation of the terms which arise by writings, y, ^, x^j y\ / in 
succession for x. 

We also suppose that R becomes a function of the elements a», and of the time by 
the substitution of the values of the coordinates as functions of these quantities. The 
La Grangian equations for the variation of the elements then take the following form: 

[a.-.]-^/+K-J^/+ +K..]-^--=^^ (21) 

(*=i; 2; 3; .12) 

In determining these La Grangian coefficients I make use of the known theorem that 
they are constants. This theorem rests, however, upon the supposition that the expres- 
sions for the coordinates actually satisfy the differential equations. As a matter of fact 
they are developed in series, so that it is impossible that this condition shall be rigor- 
ously fulfilled. It is therefore important to determine finally to what degree of approx- 
imation the constancy of the coefficients in question is assured. For our present purpose 
I have deemed it sufficient to develop the constant portions of the coefficients so far as 
the degree of approximation to which the solution is carried would admit. 

This theorem of the constancy of the coefficients enables their values to be formed 
with much greater facility than if it were necessary to actually perform all the multipli- 
cations which enter into the equation (20). 

In considering the general theory a yet farther simplification may be made- It is 
known that if we take six of the elements, which we may call 

at pleasure, we can always substitute for the other six certain functions of them, 

^1,(^2 . . . . C^ 

such that the La Grangian coefficients shall reduce to zero in all the combinations 
except those six in which an l^ is coordinated with its corresponding r^, in which case 
the combination will be reduced to unity. We shall make frequent use of this theorem 
in the course of the present investigation. 

Such a selection of elements I shall hereafter call canofiical 



§5. 
THE PROBLEM OF THREE BODIES. 

The preliminary step now is to express the six coordinates, x, y, z^ x\ y\ ^, in terms 
of twelve arbitrary constants, so as to satisfy the differential equations (16). 

We recall that two processes of integration are necessary in the preliminary inves- 
tigation, namely: 

1. We integrate omitting £1^ in the equations (13), writing H^ for /2 + R. 

2. We include /ij, but omit R, as in the equations (16). 
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The first step involves no difficulty, since it gives rise only to the equations of 
elliptic motion — 

(a) Of the Moon around the Earth; 

(b) Of the center of gravity of the Earth and Moon around the Sun. 
Putting in a general way 

^1? ^*2> ^3* M' ''2i '3 

for a set of canonical elliptic elements pertaining to the Moon, and 

r' o' p' V V ]' 

for those pertaining to the Sun, we shall then have, as the result of putting /2i=zo 

X, y, and^=/(c^i, e^, c^, /„ /g, 4, and 

(22) 

x\ y\ and / =f (c\, c\, c',, l\, l\, 1\, t) 

Next, to take account of /2|, we substitute in .O^, for x, //, Zj x\y\ /, their expres- 
sions (22). Then, using canonical elements, the form (21) will become 

dCi rfif2i dli diQi 

dt ~ dli dt~ dCi 

(i=i,2,3) (23) 

dc'i^da^^ dVi__d£\ 

dt - dl\ dt ~ dc\ 

These equations are now to be integrated by successive approximations. We have first 
to integrate the equations for c and I on the supposition that c' and V are constant. 
Next, by substituting these values of c and I in the equations for c' and /' and integrat- 
ing, we obtain approximate expressions for these quantities. If necessary these new 
values of c^ and V can be employed to reintegrate the first set of equations, and so on. 
The first integration has been performed by Delaunay in his great work Theorie 
du 9nouvement de la Lune. We may consider that, in this work, Delaunay has implic- 
itly formed expressions c and I which substituted in (22) give values of x, y, and z that 
completely satisfy the equations 



M2 



d^x rf(/2o+/2i) . x^ X 



dt^~ dx 



on the supposition tliat c'< and I'i are treated as constants. More exactly, if we put, 
in (14) 

rj^f ^2 ^^^ 

r 

(25) 
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Delaunay's values of the coordinates, x, y, and z, of the Moon satisfy the differential 
equations 

^Ut'- dx 



with similar equations in y and -2^, provided that we regard /2'o and O.^ as functions of 
Xj y, z^ rr', y\ z'^ before differentiating as to x^ and afterwards substitute for x^ y^ z^ their 
expressions in terms of c< and Z<. 

Since fif'^^ when regarded as a function of the six coordinates, does not contain Xj 
y^ or ^, we may add it to if2'o+»^i i^ the equations without altering the result. We 
thus find that the equations 

(Px_dn 
^Ue- dx 

will be identically satisfied with regard to all twelve constants in Delaunay's theory on 
the supposition above made. 

It must be noted that c and I in these equations have not the same signification as 
in (22), but are now the arbitrary constants mtroduced by the integration or*(23). 

To complete the solution of the problem of three bodies we have next to integrate 
the second set of equations (23). We therefore substitute in D.^ the values oi x^ y, z^ 
od ^ y , / as functions of all twelve elements c, c', ?, V . But it is to be noted that at this 
gtage we are to differentiate /2i with respect to d and J! only so far as these quantities 
are introduced through x\y\ and /, because, from the mode of forming the equation? 
(23) we have 

dfi^_d0.^dx'.dflidy*dflidi^ . ,x 

rf? ~7Z d7^T^ dJ'^T7 dd ^^ ^ 

c' being any accented element. 

Since i^i contains the lunar elements, c and I. the integration of the equations 

~dt-JT, ^*''- ^'7^ 

will give c' i and V ^ as functions of r<, /< and six new arbitrary constants, which we may 
again represent by c' i and I' i. By substituting these values of d i and I'i in (22) our 
three equations 

d'^xi d£l ^ 
^ df^d-^^^''^ 

will then be identically satisfied with respect to all twelve constants, provided that in 
representing the identity we first regard /2 as a function of the coordinates expressed 
in terms of the twelve constants. 

By this process we have for the first time introduced the lunar elements c and I 
into the expressions for the Sun's coordinates. Hence, in strictness, the equations (23) 
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will now no longer be satisfied with these new values of a;', y', etc. It may be shown, 
however, that the eiror thus introduced is too small to affect the solution of the prob- 
lem which now concerns us. 

The result of integrating the equations in question may be expressed in the fol- 
lowing way: By completing the solution of the problem of three bodies in the manner 
here proposed we shall implicitly obtain values of the six coordinates rr, y, ^, x' y\ and 
s/ as functions of the quantities 

^Ij C2j ^3? 'ij *2> '3? C 1> ^ 2j ^37 ^ 1? ' 2? ' 3? ^ • (^} 

which completely satisfy the differential equations 

d^x _d£l 
^"^ df - dx 

d^x' _da 
^'d¥-dx' 

etc. etc. 

the two members becoming identical if, after differentiating /i, we substitute for the 
coordinates their expressions in terms of the elements (a). 

Regai-ding iCl as a function of the constants c, /, c\ /', and of the time, we may indi- 
cate the satisfaction of the equations by expressing them in the form 

^^df- \dc,dx'^dUd'x^~dc\dx'^dr\dx] ^ ^ 

with the five additional equations formed by writing successively y, z^ x',y\ z' in place 
of X in this equation, the seventy-two partial derivatives of c^, /., etc., with respect to rr, 
a?', etc., being determined by the seventy-two equations each of one of the forms 

%M dx dCi.dx dli dx dc\. dx ^?\|_ 
i=i[T7, dx '^ dl- dx^dcU'd^'^ drV(TI]~ 

(29) 



i=3 



^ ^dXl dCi.dxi dJi.dx 
i^i\d c'i dx dli (l-^ dc, 



\ dCj^dx, dl'i\ 
?i dx dli ^f"^ i 



the remaining seventy equations being formed by writing in succession x\ x^ x\, ;y, etc., 
in place of x or x^ in the numerators found in these equations, and 1/, z^ x\ y\ and z' for x 
in the denominators. 
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§6. 

CONSIDERATIONS ON THE FORM AND ORDER OF MAGNITUDE OF THE 
EXPRESSIONS WHICH ENTER INTO THE PRECEDING THEORY. 

The fifteen coefficients [a^^ orj which enter into the equations (21) are each divis- 
ible into two parts, one of which is formed by differentiating the lunar coordinates Xy y^ 
and ^, the other by differentiating tlie solar coordinates x' j y\ and /, If we put 

(1 a^ daj d ai daj da^d a^ 

__dx dxi dy dy^ dz dzi 

d oL^ da^ d a J da^ da^da^ 

r^ ^1 _dx'dx\ dy^dy'i dz" d/i 
dcti daj da^da^ d a^d a^ 

_dx! dx\_dy' d}/i_dz* dz^^ 
daj da^ d a^da^ d a^d a^ 



(30) 



(31) 



\daj i~~dxdrai ~dy rfoj ~dls Joj ^^^' 

/rfR\ _ d^doi^ dUdy' dUd/ . . 

\d aji'-d x'da^'^dy'da.'^dg'da, ^^^^ 

we may form by means of the expression (17) the two equations 

i^\ ^ da. ^dR\ . . 

i=^f. -, da. ^dR\ . . 

and the sum of these two equations will make up the equation (21). 

We may divide the twelve elements a into two classes, which may be designated 
lunar and solar elements. The first class comprise those which, when the Sun's dis- 
turbing force vanishes, merge into the elements of the Moon's elliptic motion around 
the Earth; the second, those which, when the mass of the Moon vanishes, merge into 
the elements of the Earth's elliptic motion around the Sun. 

Let us next investigate the order of magnitude of these several coefficients and of 

(^ — J and (-7—) • Ii^ the expressions (17) for the solar coordinates we are to dis- 
tinguish two portions, one a purely elliptic part containing only the six elements of 
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the elliptic motion of the Sun, the other part containing the terms added when we take 
into account the distance of the Moon as compared with that of the Sun. The first 
part arises from the integration of the diflFerential equations for rr', y, / when only the 
part *f2o of the potential is considered, and is equal to the mean distance of the Sun, a', 
multiplied by purely numerical ratios. 

In the expressions for Xj y, and z, all the terms containing solar elements will be 

multiplied by — . 

We shall use the symbols S, S', etc., to represent in the most general way the sum 
of a series of sines and cosines of angles, each multiplied into purely numerical ratios. 
These symbols will not twice represent the same quantity, and may, in special cases, 
reduce to purely numerical constants. If we put for brevity 

n 
n 

(36) 

the expressions of the form (17) and (18) for x^ x^^ etc., and their derivatives may be 
put in the form 

X -a (S+>hS') 

:r'=a' (S' + AS) 

Xx ZZ>U2«^(S+WS) 

a/izz>Uia'n'(S'+fcS) 
etc. etc. 

If, in forming the derivatives of these expressions as to the elements otj . . . ajg, 
we choose for the latter purely numerical quantities, which we may always do by 
replacing a line element by its logarithrti, then the derivatives in question will be of 
the same form and the same order of magnitude with x^ xi ^ etc., themselves, or more 
exactly with the tenns of those quantities which contain the elements. If, then, as 
before we distinguish the solar elements by accents, the order of magnitude of the 
largest terms of several derivatives will be as follow: 

dX c! 

-^ — zzab 
a a 

(m X\ c\ 

i-^ = /Yo « >^ S 

a a 

i, ^^'^ 

da 
a a 
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dx 



da' 



j=maS' 



dxi 



da! 



J =r;/2»tonS' 



(38) 



(39) 



da' 

da 

The various classes of coefficients [a^, a,] formed by (30) and (31) will therefore 
be of the following orders of magnitude: 

[a, tjr'Jizz/i2^«^wS 
[«', ar'<], zz //2 w^ a^ w S 

[a, a^JaZZ/^ii-^a'^w'S 
[a, a' J2 = >Ui ^" a'^ w' S 

[a', a^]2Zz>Ula'^^^'S 

Let us next revert to the perturbative function R as found in (14). We see that 
it is composed of two parts of very diflFerent magnitude. Considering all linear helio- 
centric coordinates as of the same order of magnitude, a\ and the coordinates of the 
Moon as of the order of a, the order of magnitude of the two terms of R will be 

or, if we use coordinates instead of mean distances, the orders of magnitude will be 
represented by 

R— ^^J S+^-?^^— g 

where the first term does not contain the Moon's coordinates. 

From this equation we shall have, using the sign = only to indicate identity of form 
and order of magnitude, 

d R m^jUzX ^ m^ A'2 « a 

~dl^-1^'~ af''~ 

:;7^— ":.72~ ^"1 — ^4 " ^— ~;r'2 ^"i — V4 — ^ 
ax X X a a 
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and thence, from (3 2) and (33), having reference to the orders of magnitude of — ^, etc., 



rd^\ _m^lUik^ 
\daJ2~ a' 

\d aOi ~ a'^ "" a' 

\daO2- a! 



(40) 



The approximate values of yui and //g are 

/^2 = ^ 

by which we shall replace them in R. The relation* between the mean motions and 
mean distances give, approximately, 

2 ^^h 

or n =, — - 

an 

a'^n' = -j^^ 
a n 

Substituting these values in the expressions (39) for the coefficients [«!, aj, and 
representing by Si, S2, etc., Ui, Ug, etc., numerical quantities like those we have hitherto 
represented by the general symbol S, we have 

an 
an 

[a, «',-],="']'"? ksr 
a n 

a n 



(izzi,2 ... 6) 



(41) 



a n 
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The expressions thus formed are to be substituted in the pair of general equations 
(34) and (3 5) which, it must be remembered, may belong to either of two classes, accord- 
ing as a^ there represents a solar (accented) element, or a lunar (unaccented) one. Let 
us consider the former case first. The right-hand members of these equations, taken 
from (40), may then be put in the form 

\da'Ji a a a' 

The pair of equations then become 

an [ at at at \ at aty] 

a 

(42) 



a a^ 

The sum of these two equations gives one of the rigorous La Gbangian equations 
for the variations of the elements, in which each coefficient of the differential of an 

element is a constant. The coefficient of -^ in the sum is 

at 



an an 



nu . an' 



Here the coefficient of S< is greater than that of U< in the ratio — -7- -, -t», or about a 

W3 an 

hundred millions to unity. The constant terms of Si*^ must therefore be to the periodic 

terms in at least this ratio. The coefficient of Uf ^ may therefore be entirely neglected 

along side that of Si^\ 

Again, the approximate value of the coefficient k is 



iWga'^' i2,ocx:),ooo 



the terms multiplied by this quantity may therefore be entirely neglected. The varia- 
tions of the solar elements may therefore be determined entirely from the first of the 
pair of equations by a system of six equations of the form 
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S»'^+S?>^^+ . . . +Si'>^«=^*n'S'« (43) 

at at at nil 

which will, in fact, be the usual equations obtained on the supposition that the masses 
of the Earth and Sun ai-e combined in their center of gravity. This hypothesis is 
known to suffice even for the action of the Sun, a fortiori must it suffice for the disturb- 
ing action of a planet. 

Let us next consider that in the pair of equations (34) and (35), «;• represents a 
lunar element. If we then express the first two equations (40) in the form 

\daJi a' 

(44) 
(^R\ _m2m4^g(7, 
\day2 a! 

if, also, we substitute for the coefficients their values in (41) the pair of equations (34) 
and (35) assume the form 

an \^ ^ dt ^ dt ' '. J an \ ^ d t dt 'J 

a' 

a n \ ' dt dt ^ a n \ at at y 

a' 

The sum of these two equations forms the final equation for the vai'iations of the lunar 
elements after the variations of the solar elements are determined from (43). If we 
multiply their terms hy a n-^nu^fihj putting for k its value 



k= 



Wg a' 



ni^+m^a'^ 

using the relations 

(X^ n^ =^2+^3 

a'^ n*^ = m^ + ^i2 + ^3 
and regarding as unity the ratio 

mi'7'mi+ni2+nh 
the equations become 
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dt dt dt dt m^-\-ni^a^ 

(45) 

^ dt dt ^ dt dt ft^+m^a 

The sum of these two equations will give one of the La Grangian equations for 
the variations of the elements, in which the coefficients 

and 

m(Uf>+sfo 

are functions of the elements alone, and do not contain the time explicitly. 

In these sums, however, the coefficient of S^*^ in the second equation is less than 
that of U* in the first in the ratio of m k to i , and since the periodic terms in the 
two coefficients destroy each other, the quantities \J[^\ U^*\ etc., must be constant to 
their 1 50,cxx),ooo part, and may therefore be regarded as absolute constants. In fact 
these coefficients Ui*' are the same as those which determine the variations of the lunar 
elements in what I have called the first method, and which is that employed by Hill. 

Thus the final equations for tlie variations of the lunar elements may be written 
in the form 

•Q(4) ^ ^1 I jj{Ay^j_f^i , ^ _j_ u<®> da^ ___ m^n a' mc?) i g(7)\ 

^ dt ^ dt • ' ' * dt m2+ni^a'^ 

- (Uf > + Si^>) m ^^ - (Uf + SrO ni ^ - etc. 
those of the solar elements being first derived from (43). 



(46) 



§7. 
METHOD OF EXPRESSING THE COORDINATES IN TERMS OF THE ELEMENTS. 

Our next step is to investigate the formation of the coefficients [a^, a^ in (20). 
This presupposes that the rectangular coordinates x^ y, z^ x\ y\ / are expressed in 
terms of the twelve arbitrary constants a^ , . . ^12 and of the time, in such a 
form as to be a complete solution of the differential equations (16). This requires 
that the coordinate planes and lines of reference be arbitrary, but no loss of astronom- 
ical rigor will be incurred if we suppose the fundamental plane to differ from that of 
the ecliptic at some epoch only by a quantity which may be treated as an infinitesimal. 

The theory of Delaunay gives the Moon's longitude and latitude referred to the 
plane of the ecliptic and expressed in terms of the arbitrary constants 

a, e, y, a',e',h, g, I, h!,g', V, 

of which, however, one is superfluous. 
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The notation which I shall use is, however, different from that of Delaunay, the 
symbols given at the end of each line being those of Delaunay. I put 

v; the Moon^s true longitude, from an arbitrary origin . . . (V in Delaunay) 
/?; its latitude above the ecliptic (U *^ /' 

I] its mean longitude (^h+g-^-l " 

g; its mean anomaly . . . (I ** 

6; the longitude of its node on the ecliptic (h " 

y; the sine of half of the inclination of the Moon's orbit to the 

ecliptic (y " ). 

The corresponding accented quantities refer to the coordinates a/, y', / of the 
Sun. 

Delaunay's work now enables us to express the spherical coordinates of the Moon 
in the form 

vzul+Sk sin N. 

rzz2J(f cos N. (47) 

/3=2c sin N'. 

The coefficients k, kf and c are functions of five independent elements 

a, e, y, a' e\ 

while N and W are linear functions with integral coefficients of four angles, 

. ff, /, ^j ^' 

CO and oof being the mean distances of the lunar and solar perigees from the Moon's 
ascending node. Thus the coordinates of the Moon are expressed in terms of ten 
independent quantities, of which, however, one is only a constant additive to I or to v. 

It is to be remarked that, so far as the problem of three bodies is concerned, the 
position of the ecliptic is to be regarded as invariable. 

If w^e represent by ^, 9, ^, the rectangular coordinates of the Moon referred to 
the ecliptic, and put, for brevity 

Sv=v—l=2k sin N 

the expressions ^, 7 and ^ may be written 

Sznr cos /? cos v=:r cos /? cos Sv cos l—r cos /S sin Sv sin I 

rf=ir cos /? sin v=ir cos /? cos <J v sin l-\- r cos ft mi 8 v cos I (48) 

^=:r sin ft. 

Since ft and dv are small quantities, the limit of whose values is little more than 
o"i, the three quantities r cos ft cos <Jt;, r cos ft sin dv and r sin ft each admit of bemg 
developed in a convergent series, as follows: 

A. p.— VOL. V, PT. in 3 
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We first substitute for the sines and cosines of yff and S v their developments 

sinot;=ot7( I — — 5 f-+etc. . . . ) 

V 3' 5' 7' ^ 

cos ov=i r + — T FT + ^tc- 

2! 4! 0! 

sin /?=/?fi — ^,-+-i — etc. . . .^ 
cos /?=i— -, +- i — etc 

2! 4! 



Forming the powers of <J t' and ^ and substituting for the powers and products of the 
sines and cosines their expressions in terms of the sines and cosines of the multiple 
angles, we see that the three quantities in question can be represented in the same 
general form with 6 1;, r, and /?, that is, in the form 

r cos fi cos 8 v=.2h cos N 
r cos ^ sin S vzuSK sin N 
r sin I3z:z2 d sin N' 

A, /^' and N being of the same form as A; and N in the values of ^;, yff and r m (47). Thus 
the values of ^, 7, and ^ in (48) assume the form 

^=i-SAcos(/+N) 

7=-2A;sin(/+N) (49) 

^=-2csinN' 

where 

k=h±h' 

and N is still a linear function of the four quantities 

Sly ff^j ^y <^'' 

The corresponding coordinates of the Sun may be expressed in the same general 
form, which we may write 

^'=-SA/cos(i'+N) 

7'--2A/sin(r+N) (50) 

^'=-Sc'sinN' 

It is however to be remarked that, in the case of the Sun, the principal values of }^ 
will correspond to cases in which N depends on g' alone, and that when g^ co or co' 
enter into N or N^ the values of Id and c' will contain the minute factor A; of § 6. 



PRODUCED BY THE ACTION OF THE PLANETS. I3I 

Let US now suppose the positions of the Sun and Moon referred to any system of 
rectangular coordinates whatever. If we represent by 

a, b, c; a^ b', c'; a", b", c", 

the nine direction cosines for transforming ^, 7, ^ into a;, y, z^ and if we put 

&j the ascending node of the ecliptic on the plane of XY counted from the axis 
ofX, 

y'^ the sine of half the inclination of these planes, 

T, the distance of the same node from the axis of B^ 

we shall have 

a zzCi—x'^) coses' — r)+/2 cos (e' + r) 

b =_(i_/2)8m(©'-r)+/28in(0'+r)' 



c =2 )/ Vi— /^sin S' 

a' =(i-/^sin(©'-r)+/2sin(0' + O 

b' =(i-/2)eos(0'-r)-/*-*cos(©' + r) (51) 



c' =— 2 yW I — y'^ cos & 



a'' =— 2 / V I —y^ sin r 



b" = 2 / V I —y'^ cos r 
c" = i-2/« 

By means of the values of ^, 7, and ^ given in (49), and the above values of a, b, 
etc., we find for the general rectangular coordinates 

a;=-SA(i-/?)cos(^-r+/+N) 
+ 2ifc/^cos(©' + r-/-N) 



-2c/Vi-/'cos(©'+N0 



+2c/Vi-/'cos(©'-N0 



(52) 



-2 c / V I -/' sin (0' +N') 



+2c/ Vi-/'8m(0'-N') 
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z=2Sky Vi-/'8in(Z-r+N) 
+2(i-2/^)cBmW 



-Sc'yWi- y^ cos (e'^W) 



+2cywi-y^cos(e'-w) 

i/=2I(f (i-/^) Bin (&-t+r+^) 
+2kyBm(e' + r-l'-JS) 



(52) 



-2&yWi-y'sm(&+^') 



+2c'rWi-y'sm(&-W) 



y = 22^/Vi-/'sm(r-r+N) 
+2(i-2/2)c'8inN' 

Examining these equations, we see that the coefficients of the varying angles are 
now functions of the six constants 

a, e, r, a', ^^y' 

and that the angles under the signs sin and cos may be expressed in terms of the time 
and six other arbitrary constants. For these angles contain only & ^ ^— ^ N, and N'. 
The two latter can be expressed as a function of the four quantities already mentioned, 
T—l need contain but one constant, and thus, with & ^ we have the required number 
of six. 

To reduce the preceding values of x*, y, and ^ to a common form exhibiting these 
six arbitrary constants, let us represent the four varying angles of which N tod N' 
are functions by A, A', oo, od\ These quantities are not to be considered as having a 
perfectly defined geometrical signification, being only certain analytical expressions 
which substituted in the equations (47) will give numerical vahies of the Moon's coordi- 
nates. But, we may consider that they represent 

A, the mean distance of the Moon from her node on the ecliptic. 

G?, the distance of the perigee from the same node. 

A', Go\ the same quantities for the Sun, counted from the same node. 

♦ We shall now express all the quantities found under the sign sin or cos in a;, y, and 
z in terms of the six arbitrary quantities 

i?, n, e, e\ re, & 
the first five of which are determined by the conditions 
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e—ff + l—T 

e=s—x 

c' = A' — A+€,z=(9' + r — r 
7r'=G+a}' 

They are therefore introduced by putting in the expressions for a;, y, and £, 

l—T=e—& 

A =e—d 

, . (53) 

A'zze'— © 

w=.n — B 

m'zun'—e 

The expressions for x, y, and z thus take the form 

x=2k cos (ie+i'ir-\-i!^e-^i'^e'-\-i'^7r'+f&) 

yz=2k sin (ie+i'a-+i«5+i"'«'+t^;r'+i^0') (54) 

jizzSc sin (je-\.f7r+j^e^j"U'+r7r'^j''e') 

the coefficients * and j being positive or negative integers satisfying the conditions 

and A; and c being functions of the six quantities a, of, e, e', y, y'. 

The values of a;', y', and i/ in (52) can be reduced to the same form, so that it is 
unnecessary to write them. 

§8. 
METHOD OF FORMING THE LA GEANGIAN COEFFICIENTS. 

It will be seen by what precedes that our problem requires as a preliminary work 
the complete solution of the problem of three bodies by obtaining not only the rela- 
tive coordinates of the Moon and the Earth, but the coordinates of the Sun relative 
to the center of gravity of the Earth and Moon, which has not yet been done. Instead 
of effecting the*solution of this preliminary problem at the present stage, I shall pro- 
ceed on the supposition that this part of the work is done, and proceed with the 
investigation growing out of it, in order that more exact information may thus be 
obtained as to the degree of approximation to which the solution must be carried. 
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It is shown in my paper of 1875, on the General Integrals of Planetary Motion,* 
that when the solution in question is carried on by successive approximations, the 
expressions for such coordinates as x^ y, z^ od ^ y', /, in terms of the twelve arbitrary 
constants which were required in the solution of the problem of three bodies, could 
always be reduced to a general form of the kind found in the last section, namely 

X or x'-=z2k cos N 
y or y^=2 k sin N 
z or /rz-2 c sin W 

in which each value of N and of N' is of the form 

W=ji X1+J2 A2+ . . . +j\ Ae 
ii, i^, . . ie, Ji . . . je, being integers satisfying the conditions 

J1+J2+ . . +J3=l , 

and Ai . . . Ag expressions of the form 

X,=h+h, t 

liy Zg) • • • ^6 being arbitrary constants, and 61, 62? • • • h as well as the 
various values of k and of c, functions of six other arbitrary constants, which, if used 
as in the theory of Delaunay, may be called 

a, e, y, a\ e', / 

The symbol 2 indicates a summation extended to all combinations of integral values 
of the indices ii, ^2, etc. 

As special values of the six constants / we may take those which in the ordinary 
theory are designated as follows: 

0, the longitude of the Moon's node, or Delaunay's h\ 
TTy the'longitude of the Moon's perigee, or Delaunay's h-^g] 
ly the Moon's mean longitude, or Delaunay's h+g+l-, 
& y 7t\ I'y the corresponding quantities relating to the Sun. 

It is to be remarked that the assignment of the usual astronomical designation 
to these quantities does not adequately define them, and may therefore be misleading. 
No precise definition is possible except by their analytic expressions as integrals, each of 
them being regarded as a function of the coordinates, their first derivatives, and the 
time. This definition might, however, be considered as replaced by any complete 
solution expressing the values of the coordinates in terms of the arbitrary constants 
and of the time. We shall therefore simply suppose that expressions of the preceding 

'Smithsonian Contributions to Knowledge, Vol. XXI. 



PRODUCED BY THE ACTION OF THE PLANETS. 



^35 



form for the rectangular coordinates are obtained in terms of two sets of six arbitrary 
constants each, which satisfy the diflferential equations for the case of three bodies to 
the necessary degree of approximation. 

The danger of error in ai^signing geometrical significations to the elements will 
be seen by considering the special quantity which Delaunay has called the longitude 
of the Moon's node. His entire lunar theory is constructed on the supposition that 
the plane of the ecliptic is invariable, an hypothesis which is sufficiently correct when the 
action of the planets is omitted, so far, at least, as any secular variation is concerned. 
But when we consider the complete problem of three bodies we shall find that there 
are two planes, either of which may be called tliat of the ecliptic. One of these is 
the plane ^f maximum areal velocity, which is really invariable. The other is the 
plane determined at each instant by the relative motion of the Sun and of the center 
of gravity of the Earth and Moon, which plane is variable. 

Since all general equations between physical quantities are necessarily homogene- 
ous when each quantity is expressed in units of fundamental or arbitrary elements, it 
will generally conduce to clearness to write them in such a way that this homogene- 
ousness shall always be manifest. In physics the three fundamental elements whose 
oyiits are arbitrary are mass, length, and time. But in astronomy it will generally be 
simpler, as it is more common, to regard the unit of mass as a derived one, and 
determined by the condition that two bodies of unit mass shall attract each other with 
unit force when at unit distance. If we represent by the letters M, L, and T the 
dimensions mass, length, and time respectively, tlie equation for M to which we are 

thus led is 

MzzL«T^2 

In a theory like the present it will conduce to clearness to write all equations in such 
a way that the homogenity shall be apparent. Since N and N' are angles, and there- 
fore pure numbers, it follows that A, I, and b t are also pure numbers. Hence h is of 
dimensions T~\ as are also n and n\ the mean motions. Since the coordinates are lines 
all the coefficients k are of the form 

a or a' X a pure number 

while 61, 62 • • • &6 are each of the form w 5, or n' 5, s being a pure number. 
Hence each conjugate variable of the kind 

dx 

is of the dimensions M L T~^ or L* T^ according as we do or do not retain M, and 
may be written in the form 

m a w X a pure number. 

The six pairs of conjugate variables from which the La Grangian coefficients are 
to be formed may therefore be written in the form 

x-^iS ki cos N^ Xiz=:—2 k'j sin N^ 

y=2kiBm}!ii yi±z 2kfjC0&^j (55) 

;s =: 2 Ci sin N^ ^j zz -5" c'^ cos N'^ 
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where the indices i and j are used to distinguish the terms belonging to the various 
linear combinations which make up the values of N. We have, in each case, 

pi representing either yui or /ig, according as the coordinates are those of the Moon or 
of the Sun, and b and V being the coefficients of t in N and N', so that in general 

(56) 
V=jibi+J2b2+ . . . +jeb9 

Now, let a and /S be any two elements for which we are to form the La Gbangian 
coefficient [a, /?]. Putting no restriction as to the class to -which either of these ele- 
ments may belong, we have from (55) the following values of the derivatives from which 
[«i, /?] is to be formed 

P=2 (-^-Kj sin N,-A:', ^-^J cos N,) 
da \ d a ^ ^ da v 

dy ^/"^^i • XT 1 I ^Ni XT \ 

-,-^-=-2 ( — -* sin J^i+ki ^ -* cos N< ) 
da \d a da J 

dyi ^ fdli . XT 1/ ^N. . XT \ 
V^=-2 ( -—J cos Ni— A i — ^^ sin N, ) 
dp \d/S ^ ^d/3 V 

^=2 (p sin W,+c, -"i:- cos N',') 
da \da da J 

dz^ ^fdc!^ XT/ f d^\ . xT/^ 
3--ZZ-2 ( — - ^ cos N ,.— c . -, -^ sm N '. ) 
da \da ^ ^ da V 
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Substituting these values in tile expressions of the form (30), p. 122, and putting 
for brevity 



P ^1 dq)dd dq)d6 



we find 



(57) 



[a,/?],=z 2''{[*„A:'J+A;,A:',[N„N,]}sin(N,-N,) 
+ 2' \h [N„ hf,]-k/jlk, N,] |co8 (N,-N,) 
+i ^M [Ci, <^'j] +Q C, [N'„ N',] } sin (N',-N',) 

+i 2^ { [c„ c',]-c, c'j [N'„ N'J } sin (N',+N',) 
+i^{c,[N'„c',]-c',[Q,N'Jlco8(N',-N',) 

+i^-^{c,[N'„ r',]+c',[c„N'J}co8(N',+N',) 

It will be understood that the two indices i, j, have each the same system of values, 
and their combination under the sign -5^ indicates that each value of N is to be com- 
bined with all the other values, itself included. 

In order to find the function [a, /?] of which we are in quest, it is necessary to omit 
all the periodic terms in (57). Substituting for N^ and Nj their values /+Ai< and 
fj-{-bjt, we see that all the terms will be periodic, unless we have either 

bi —bj =zo 
or 

b\±b'j=o^ 

The values of bi and bj being from (22) 

fc< =«i.J>i+H<^2+ . - . +\ibe 
bj = iij 61 + ^2,; &2 + • • • + hj be • 

b\=Ji,ih+J2,ib2+ . . . +Mibe 
Vj=jijbi+J2jb2+ . . . +J6jbe 

it is necessary that we have either 

(hi—hj) bi + (h,i—i2j) 62+ . • . +(\i—kj) &6=o 
or 

But 61 . . . fee being incommensurable quantities, these equations can not be satis- 
fied except we have 
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from which follows 

Therefore we retain in the equation (57) only those terms which satisfy the condition 

i=j 
We have, therefore, omitting the index i, 

[a,/3l= ^{A[N,*<]-A/[A;,N]} 

+i2{c[N',c']-c'[c,N']} 
or, placing for f its value /i^ ^ k, and /J2 b' c for c', 



1/3 da da d /3 . 

2b'c')dW _d( M2b'c^)dW 
Ifi da da ' d/S 



"'"2 I d0 da da del 



Omitting in - — the terms which contain < as a factor, we then have 
\a B^-u ^VdfdQ>T^)_dfd(hl^)-\ 

i u ^rdf d{b'<») _df rf(6V)-| 
2^^ \-da dp d 13 da'J 



(58) 



It will be remembered that a and /3 here represent any pair of the tw<?lve arbi- 
trary constants which enter into the expressions (55) for x^ y, and z. These con- 
stants are divisible into two classes. 

Class (i) comprises the elements 

a, e, ;/, a\ e\ /, 

which enter into, 6, 6', A;, and c, but not into/ 
Class (2) comprises the elements 

^1? ^2? • • • • ^6 

which enter into/ and/, but not into 6, &', c, or k. 

It is evident from the form of the preceding expression for [a )ff] that it will 
vanish whenever the elements a and /? belong to the same class. We. have therefore 
only to consider the combinations in which they belong to different classes, of which 
the number is 18. 
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Let us now conceive the general quantities A^ . . . A, to be replaced by the 
six elements 

e, TT, e, £', ir', &, 
already de5ned, and let us put 

n; TT^, . . . . ©1', 

for the motions of the corresponding elements in unit time. Each value of 6 or V may 
then be written 



We then have 



h =i n + i'n-i + i" e,+i'"n'+i'^ ^i+i^ ©'1 

b'=j n+r ^x+y" Oi+rn'+r ^\+f &x 

df_d^_. df _dW_. 

3T"~ de ~~ de ~ de 

i/_<^N_., df_dW_, 

dv d TT dTT dx 



df_d^_^ df _dW_. 

de'~d6'—'^ dQ'~d&~^ 



We have therefore, supposing ^ to represent in succession the six elements of the 
second class and restoring the sign 2, 

r T v-.<^(&^^ I . ^d(h'<^ , . 



,.,e,.=-^..llP-l^.r'-f/> 



The sign 2 liere refers to the several terms which enter into the expresssions (54). 
The entire number of sensible coefficients [a, /?]i is now reduced to thirty-six. 

To complete the required expressions we still require the values of [a /?]2. The 
rigorous values of af^ j/, and / can be developed in precisely the same way with 
those of X, y, and z. As we have seen all the terms containing a, c, y, «, ;r, or 6 will 
contain the minute factor k in (36), but for the sake of generality we shall suppose 
them all retained. If we put 

x'zz^k cos (I e+V ;r + P0+r"£' +rv;r' + r0') 

j/=2ksm (le-^VTT+re+Ve' + V^Tr' + re') 

/=S€ sin(j£+J';r+J"0+J"'£'+J^^';r'+J^©') (60) 

B=l n+V 7r, + V^, + V'W + J''' 7r\ + r e\ 

we have as in the case of the lunar coordinates 
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r T vT dCBK') I ^. dCB'C^ 



Finally, we have for the complete La Grangian coefficient 

[a, /?] = [«,/?], + [«, /?]2 

We conclude that the values of the La Grangian coefficients may be obtained as fol- 
lows. Let us determine the six quantities A., fc,, k^, A;/, kj, k^ from the equations: 

K =pix^l BK^+%i-SJ B'C'+/^2-^'i bk^+^-M22j h' c^ 

2 2 

2 2 



2 2 

A;/ - /i, ^^ I'" B K^'^- - /^i -^^ J"' B' 0^+tii2 1'" ?> A-* + I ^^2fn' <? 

2 2 



V = >«l^I''HK*+Vl2J^B'C='+yUj^r/>A-*+lyMj^VV>V 
2 2 



We shall then have 






[/.]=-J|,-; [/']=-^*;- 



(62) 



(fi3) 



PRODUCED BY THE ACTION OF THE PLANETS. I4I 

$9. 

EBDUOTION OP THE EQUATIONS TO TAB CANONICAL FORM. 

Let us now resume the equations (21), supposing the symbols a to represent in 
succession the twelve arbitrary elements 

a, e, 7, o', e', /, «,, ^To, ©o,«'o, ^\, ^0 

where the six elements £01^0 represent the constant parts of e, tt . . . . 

Observing that all pairs of constants, both members of which ai'e found either among 
the first six or the last six of these constants, will make the combination (a^, a) vanish, 
and that we have 

rfR rfR dR dR ^ 

a tf) ae aTC^ dTr 

the equations (21) by substituting the values of the coefficients in (63) become 

dk,d a.dk, de .dk,dy ,dk, da' dk, de' ,dk, dy' _d R 
da di d'e dl'^ dy dt dl,' di "^ di di'^dy' It ~ de 

dk, da .dk, de .dk, dy.dk, da' .dk, de' .dk, dy' d R 

daUt Je dt^dy di^ da' dl ~^d7 dt'^dy' di ~dlr 

dkfda .dkt'de .dk, dy .dk, da' .dk, de' .dk, dy' _ d R 
Ja d't^de "dt'^dy dt^ da' dJ '^ de' dJ~^dy' ~dt~~de 

dk-l da^ dk.' de dk,' dy dk^ da' dk^ de' dk^ dy' _dR 
da dt'^ de dt'^ dy dt'^ da' dt "*" rf e' dt^ dy' d i ~ de' 

dk,' d a dk^ d e dkj d v dk,' da' dk,' de' dk,' dy' _dR 
da dt'^ de dt^ dy dt^da' dt^ de' dt "*" dy' dt ~'d7r' 



dk^da (Ik,^ de dk,' dy dk^ da' dk^ de^ dk^ dy' _dR 
da "di'^ de dt^ dy dt^ da' dt '^ de' dt "'' dy' dt ~ d& 

dk, d e^ dk, djTo.dk, dO^ dk^ de'„ dk,' dTr'^ dk,' d '>'o__ d R 
da di da dt da dt da dt da dt da di ~ da 

dk^dsp dk,d7To dk,de„ dk/ dt'^ dk,' dTr'^ dk,' dO'o __dR 
de dt de dt de dt de dt de dt de dt ~ de 

dk,d£„ dk,d7ro dk,d9^ dk.'de'p dk,' drr',, dk,' dO'o__dR 
dy dt dy dt dy dt dy dt dy dt dy dt ^ dy 

dk.de^ dk,d7r^ dk,d9o . dk,'de'^ dk,' drr'^ dk,' d6'^ __dR 
da' di'^Ja' di '^dd' ~dJ'^da'~dt^Ta' df^ da' dt~ da' 

d,k.d £o,dk, djTj dk, (1% dkj^ d^.dkjdjr'^dk^ ^^'o—_d^ 
de' dt '^de' dt^de' Tt^'dV dt "^ de' dt '^ de' dt~ de' 

dk,d€o dk,d^o dk,ddo dk,' de'p dk,' drr'^ dk,' d9'Q __dR 
dj^'lU'^dy'l^t^dy'llt-jy'Tt^dy dt '^dy' dt ~ dy' 



(64) 
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Since A^ A;,, etc., are functions of the six elements a, e, y, a', e\ / the left hand 

d k dk 

member of the first of these equations is equal to - -j that of the second to -ry and 

so on through the first six equations. The latter thus reduce to the six left hand equa- 
tions in (65) below. 

The quantities k,, A,, etc., being six in number, we may conceive the six elements 
a, 6, X, a\ e\y\ and R to be expressed in terms of these six new elements. If, then, 

we multiply the seventh of the preceding ecjuations by . - the eighth by -=-j- and so 

on to the twelfth, which we multiply by -^^ and notice that the theory of differentia- 

tion of inverse functions gives us 

dk^da dk/de dk^dy djc, d a' d k^de'dk,dy'_ 
dadk^ dedk, dydk, da'dk, de' dk, dy'dk,^ 



d K da , dk, de 



J dk,dj^.dk„da'dk^de'.dk^dy'_ 
V^'^ cf^dF.da!dK^Ti^dK^d^ 



da dk, de dk, 



with similar equations for every combination of two k elements, we find 

d €q___ d R 



dt d k. 



In the same way we obtain 



rf/To __d R 
dt dk, 

etc. etc. 

Tlie equations (64) may therefore be reduced to the canonical form 



dk,_dR 


d e^ 


dR 


d t~ de 


dt- 


dk. 


dK_dR 


dTTo 


dR 


dt ~dn 


dt ~ 


dk^ 


dK_dR 


dd. 


dR 


dt ~dd 


dt~ 


"dk, 


dK' dR 


de'. 


dR 


dt~ de' 


dt ~ 


dK' 


dkf,_dR 


dn'. 


d R 


dt ~d7c' 


dt ~ 


dk'. 


dk',_dR 


de'. 


dR 


dt de' 


dt ~ 


dk'. 



(65) 



We may regard these as the rigorous differential equations for the solution of the 
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problem of four bodies after the solution of the problem is effected for the case of three 
bodies. 

In deriving the equations (64) from (21), we have limited ourselves by no 
hypothesis respecting the signification of the elements f, ;r, etc., which appear under 
the signs sin and cos. We have only supposed the variable angles to be linear func- 
tions of six variables each of the form ao+6 f, a© being an arbitrary constant and b a 
function of the six arbitrary constants which form the coefficients A;. We may there- 
fore substitute for any of the quantities €, tt^ 5, €% 7r\ ©', any linear function of other 
independent quantities completely determined in terms of f, ;r, etc., by equations of 
the first degree. For every such system of new variables we shall have a system of 
canonical equations corresponding to (65). The number of these canonical forms is 
infinite. 

The transformation in question can indeed be performed immediately on the 
equations (65). Let us determine six new variables /i, ^, ^, Z^, /gj ki by tjie conditions 

Pi k +P2 k +Pz h +Pi h +P5 h +P6 ^6 = f 

2i ^i+«2 4+?3 k+qi h+q5 k+qti k=^o 

n ^ + ^2 k+r^ h + r^ h+r^ U+re k=0^ 

. (66) 

^1 ^1 "T~ ^2 ^ ~r ^ ^ I ^4 ^4 "r '5 ^5 1 ^6 'e ^ ^ 
Wi ll+'U.^ k + U-i k + '^i, ^4 + % ^5 + ^6 h = ^\ 

Pj ii ^» ^tc, being any constants whatever of which the determinant does not vanish. 
Let us put also 

c^zupi A\ + ji K + ri ke + Si K' + ti kJ + Ui kj 

C2=Pi K + q., k, + r^ K+s.^ kj + t^ V + W2 V * (67) 

c^znp^ K+q K + u K+Sfi k^ + k K +w« ki . 

Multiply the first of tlie equations (65) on the left by ^1, the second by ji, etc., and 
add the six products thus formed. The sum reduces by (67) to 

dc^_ dR, rfR, dR, dR,,dR , dR 



The equations (66) give 



d Co d tTq . 



Substituting these values of ^1, q^, etc., in the last equation if becomes 

dci_d^ 
dt dli 
We have in the same way 

dc2_dR d C:i_dR 
Tt~dl~' Ti~~dk^^' 
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\ 



To obtain the values of ^- let us suppose that we solve the equations (66) with 

(I V 

respect to /i, /g? ^tc, so as to obtain the values of these quantities in terms of e^ ttq, 0^ 
etc. These values will be of the form 

/jzzPi £o + Qi ^0+Ki «o+Si f'o+Ti ;r'o+U, ©'o 



P, Q, etc., being functions of p, j, etc., which are readily formed by the theory of 
determinants. From the theorem that any two determinants are identical when the 
horizontal lines of the one are identical with the vertical lines of the other it follows 
that the values of k,, k„, etc., derived from (67) will be 



kj = Ui ri .+ U2 6'2+ U3 6'3 -f U4 tJ4+ Us ^5+ Uc Ce 

Let us now multiply the first of the right hand equations (65) by P<, the second by 
Qi, etc., and add the products. We shall have 

dk__dJK 
dt d Ci 

We may therefore replace the equations (65) liy the more general ones 



dc,_dR 


dh_ dR 


dt~ dh 


d t~ dCi 


dCi dR 


d L d R 


dt~dh 


dt~ d c'2 


dCi^dR 


dh_ dR 


dt~dh 


dt~ dc^ 


dc^_d R 


dh_ dR 


dt~ dh 


dt~ dc^ 


d Cf,_ d R 


dh_ d R 


dt~ dh 


dt~ dCi 


dce dR 


dh dR 


dt~dl- 


dt~ dc. 



(68) 
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Let US, for the sake of clearness, recapitulate the process by which we have arrived 
at these equations. Considering first only the mutual action of the Sun, the Earth, 
and the Moon, we obtain the rectangular coordinates of the Moon relatively to the 
Earth (x, y, z) and those of the Sun relatively to the center of gravity of the Earth and 
Moon (x\ y'\ /) in the form 

xzzSk cos (A +6f) 

y=z2 k sin (A +fei) 

js=2c sin (X'+b't) 

(69) 
a;'z=^Kcos(^+B 

y'=SK sin (^+B 

z'zzSC sin (^'+B'0 

k, c, 6, b\ K, C, B, and B' being functions of the six constants 

and A, A', A, and A' being linear functions of six other constants 

M> »2l '3 J Uj 'bi Ui 

which fix the mean longitudes of the bodies, and the positions of the lines of apsides 
and of the nodes by any system of linear equations like (66) which we choose to adopt. 
Of the infinite series of values of A, A', A^ A\ represented by the sign 2 each will be 
of the form 

A orA = l, k + L /2+I3 /3 + I4 h+hh+h k 

(2) For each term of x or each value of k we fonn the product 

:?^f ^J_ h k' 

which, for the moment, we represent by k'. 

For each term of z, or each value of c we form the product 

which we may represent by c\ 

For each term of a;', or each value of A we form the product 

For each term of /, or each value of A' form the product 

2 1^1+^2 + ma 

(3) Multiply each k' thus formed by the corresponding value of ii, or the coeffi- 
cient of li in the corresponding value of A; each & by the coefficient of /jin the corre- 
sponding value of A'; each K' by the coefficient of h in the corresponding value of A^ 
and each C by the coefficient of /i in the corresponding value of A\ Represent the 
sum of all the products thus fonned by Ci. 

A. p!— VOL. V, PT. Ill i 
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In the same way form the value of c^ from the coefficients of 1^^ and so on to c^. 
These six values of c< will then be functions of the six elements 

a, e, y, a!, e', y\ 

and may therefore be conceived to replace them in the expressions for the coordinates. 
(4) The disturbing function R being a function of the coordinates x^ y^z^x\j/^ /, 
let these quantities be replaced by their values in term of the elements Cj, C2, . . . Cg, 
Zi ^2 • • ■ ^6- Then determine the values of these elements as variables by the mte- 
gration of the equations 

dCi_rfR rf/i__rfR 

df^dli df^ dCi 

The values of c^, li thus obtained being substituted in the expressions for the coordi- 
nates we shall have the values of the latter under the influence of the disturbing force. 

§ 10. 

TRANSFORMATION OF THE PRECEDING EQUATIONS. 

The values of the elements deduced from the equations (65) or (68) are to be 
substituted in (69) to obtain the coordinates of the Moon relatively to the Earth under 
the influence of the attraction of the planet. If we take the equations in their literal 
acceptation, we shall find terms in which the coefficients of the sines or cosines of the 
varying angles have the time as a factor. These terms^ are introduced in the following 
way: We have assumed that in the disturbed motion the coordinates are expressed in 
terms of the elements in the same form as in the undisturbed motion, only that the 
elements are variable. In the undisturbed motion each of the variables 

€y TTj S, €\ TT' 

is expressed in the form 

l + hU 

I being an arbitrary constant of the second class, and h a function of the arbitrary con- 
stants of the first class. In the disturbed motion the derivative of this expression 
will be 

dl. ±dh . . 

The expression - being formed by taking the derivative of R with respect to 

quantities which are multiplied by the time in the varying angles, will contain the 
factor t outside the signs sin and cos. 

We shall now show that these terms are destroyed by the term t -^ . Let us, as 

before, represent the angles by which replace f, ;r, etc., in the general equations (68) 
by the general symbols 

Aj, A2, A3, A^, A5, A(, 
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where 

etc. etc. 

We then have for any one A and the corresponding c 

dCi__dR dli_^dR . . 

Ji-Jx, Tt- dc, ^^^^ 

and 

The terms of R in question are of the form 

Rr=A cos (ii Ai-j-ij ^z+h h + h K + h h + h K) = f^ cos N. 
Consider any one of them. Since Ai, Ag, etc., are functions of c<, through 6<, the 

value of -^— will contain the terms 
dCi 



dRdX^ dUdX, 
dXi dCi rf Ag d Ci 



But 



a R . 7 . XT rf R 



zz — iiAsinN; ^^^rz — i2^sinN, etc. 



and 



dc^ dCi^ d Ci dCi 



Therefore the value of — ,— , and consequently of -^, will contain the term 

dCi ^ '^ dt 

z' . dbi , . dbo t . ' d b^\ . r . xt 

\ dCi d Ci d c^y 

To find the second term of -^ in (71) we have 

dbi d b^dci.d b^d C2 , ,d bidc^ 

dt~' dcxdt dc2 dt dc^ dt 

Substituting for -y^\ -^^ etc., their values in (70), and for -=^ the expression just 
dt dt d A. 

found, weliave 

,dh__f. dbi . dbi ,: ^^^sinN 
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Each terra of R will therefore produce iu -,— * the term 

'^ at 

l*-d7;-";)+<fr^'-^)+ • ■ ('"'"»'■ 

It is known that the values of b are so related to those of c that for all combina- 
tions of the indices i and j we have 

dbi dh. XX 

^—. ^—o (72) 

dCj dCi ^ 

The general proof of this theorem is found in my paper on the General Integrals 
of Planetary Motion, so that it is not necessary to repeat it here 

It appears that we may omit the term t — \Vi the value of jj in the disturbed orbit; 

provided that in taking the derivatives with respect to c, we suppose the values of 6 
under the signs sin and cos, to be invariable, and vary only the coefficients k. We 
then have in the disturbed orbit 

~dJ-dt^ ' 



\,=k+ fb.di 



This result may be expressed in the elegant form adopted by Delaunay in his 
theory of the Moon. Taking A,, Ag, etc., as the variables instead of /j, I2, etc, the 
above equations reduce the forms (68) to the following: — 

d_Ci__dR lA*— — -— +ii (73) 

dt dXi d t dCi 

Consider the diflferential function 

bi dci + b<idc2+ • . - +bQdCfi 

The equations (72) show that this expression is an exact differential of a certain func- 
tion of Cij c«, etc., which we may represent by 0. If, then, we find 9 by integrating 
this differential, we have 

. de 

d Ci 
Substituting this value of />;iu the above expressions for -=-j, and putting 

R'=R-© 
the equations which give Ci and Aj become 

~dt~dX^ 'dt~ dCi ^^^^ 

The fact that the mean motions of the Sun and Moon and the secular variation of the 
nodes and apsides of the solar and lunar orbits under the influence of the mutual 
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attractions of the three bodies can all be expressed as partial derivatives of a single 
function lends an interest to the form and properties of this function. As a prelimi- 
nary lemma in its investigation we remark that 6i, 62? ©tc , k and may all be expressed 
as homogeneous functions of Ci, C2, C3, C4, Cg, rg. If we express the coordinates x% y\ z\ 
x^ y, ^, of the Sun and Moon by the general expression 

X, y, z,x\y\i^=2k ^^^^Q + h t,) (75) 

the quantity k will necessarily be a homogeneous function of the first degree in a and a\ 
the mean distances, and therefore in the language of physics, of dimensions L, while, 
since 6 Hs an angle, and so a pure number, b itself is of dimensions T~* = M*L"'* and 

therefore homogeneous and of the order — - in a and a' when unit mass is defined as 

that which attracts with unit force at unit distance. Now, d is of the general form (62), 
that is, we have 

Cizz- 2 /di^b Jc^ 
2 

C2ZZ- 2 jLti^bl^ 
2 

(76) 

2 



ii, t'aj . . . i« being the coefficients of Aj Ag, . . . Ag respectively in the angle 
l+bty and /j being the functions of the masses employed in (62). Hence r^ consid- 
ered as a function of 

Wi, W2, W3, a, e, y, a\ e\ y\ 

will be homogeneous and of the degree - in a and a\ Hence, conversely when we 

express these elements as functions of Ci, c^^ etc., these functions will be homogeneous 
and of the second degree in Cj, ^2, etc. Hence, 

a, a\ and k are functions of the degree 2 in c^ C2, etc., 

and, since 

. dOzubi dci + b2 dc2+ .... -ffeg^^c 

it follows that is homogeneous and of the order — 2 in c^ Cz, etc. 

The same reasoning may be expressed in another way. Let us suppose the 
linear dimensions of the system increased in the ratio v:i, the masses, inclinations, and 
eccentricities all remaining unaltered. Then the analytical forms of the integrals of 
motion will remain unaltered, while the following changes will result in the values of 
the principal constants: — 

All the values of k will be increased in the ratio v:i 

, b y~^:i 

c . . v*:i 
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These conditions can be fulfilled without changing the form of b and k only when h 
is homogeneous and of the order —3 and k homogeneous and of the order 2 in c. It 
follows that each value of ky 6, and & must satisfy the partial differential equations 

d k , d k , d k , dk . dk , d k , 
dCi dc2 dC2 dCi ab^ dc^ 



Since 



d0 , d , d , d , d , d0 , ^ 

d Ci d c^ d c^ d C4 d c^ d c^ 

. d0 I. d0 ^ 

d Ti d C2 



we have, from the last equation 

= — ~(piCi+b2C2 + hc3+b^c,+b^e5+bfiCe) 
If we substitute for c„ C2, etc., thair values (76) this equation becomes 

= 2 JU (ii fei-i-l2 ^2 + ^3 ^3 + ^4 ^A + h h + h K) i> ^ 

4 
or since 

fezztifc,+t2'>2+etc. 

0zz—-2/JLb^]c' 
4 

From this formula it can be shown that is the negative of the function which 
Clausius has named the viriel of the system, that is, the constant term in the expres- 
sion for the living force of the system of three bodies. Let us represent this living 
force by T. Then 

2T=m,(5'^+7'?+n)+^(^i+7i+r20+^(^i+7i+Q 

the accents indicating derivatives as to t If we differentiate the first three equations 
(5) with respect to the time, with 5o a constant, and substitute for 5, ^, ^ their values 
in Xy y, etc., we have 

Differentiating the expression (75) we have 

^^(x, y, z, X', ^, or ^)=2Jthh'^^{f+ht) 

If we square each equation of this form, each term under the sign 2 will give in 
the square the constant term 

2 



/ 
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Hence the sura of all the constants m 2 T, which we shall represent by 2 To, will be 



2 



the sign 2 being extended so as to include all the terms in the six coordinates x',y\/, 
X, tfy z. Comparing this with the above value of & we have 

0=-To 

To being the viriel of Clausius. 

This result may be expressed as follows: Taking for the elements Ci, Cg, etc., the 
functions A., A„ fc^, etc., and expressing the constant of the living force of the sys- 
tem composed of the Sun, Earth, and Moon as a function of &„ A;,, etc., the mean 
motions of the Sun and Moon, and the secular motions of their nodes and apsides will 
be given by the equations 

d e d To_ 
dt '^'dK- 

^ . ^To . 
dt^ dK 

d9 ,dTo 
dt '^ dh.~ 

(77) 

dt '^ dk', 
di^dk'~ 

m 

d_e:_ dT 

dt ^ dk'," 

These equations can not be directly used to obtain the values of these secular 
motions of «, ;r, etc., because we can not form To until we have solved the equations of 
motion. But, the solutions of these equations being once obtained, they may form a 
valuable check on their accuracy. 
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SPECIAL CASE OF THE MOON'S MOTION 



§11. 

COMPUTATION OF THE ANALYTIC EXPRESSIONS A:,, K, ETC., EMPLOYED IN THE 

PRECEDING THEORY. 

In the preceding chapters we have shown, in a general way, how we may deduce 
expressions for the simultaneous variations of the elements of the orbits of the Earth 
and of the Moon which are due to the action of a planet. 

We now proceed to the actual calculation of these expressions. The coordinates 
of the Moon relatively to the Earth, and of the Sun relatively to the center of gravity of 
the Earth and Moon, are supposed to be given in advance, in terms of the time and of 
arbitrary constants. As already indicated, we shall take these data from the theory of 
Delaunay. In this work, however, the parallax is carried only to terms of the fifth 
order, while, for the present purpose, it is desired to carry the approximation to terms 
of a higher order. 

During my visit to Paris in April, I871, Mr. Delaunay was so good as to calcu- 
late the constant term of the parallax to terms of the sixtli order, and afterward Pro- 
fessor J. C. Adams kindly communicated his results for this same term to the seventh 
and higher orders, I am thus enabled to present the actual analytic expressions for 
k^j Jc^j etc., to terms of the order ni' and e'^ m^, and to terms of the sixth order in e, y, e' 
and m. 

ITie notation of the constants is as follows : 

The elements of the orbit of the moon, a, w, e, x, have the same signification as 
in the theory of Delaunay, n being the actual mean motion of the Moon, and a the 
mean distance, derived from n by the equation 

e is the eccentricity of the lunar orbit, and / the sine of half inclination to the eclip- 
tic. The coordinates of the moon relatively ta the ecliptic are expressed in terms of 
g, y, A, and A.', which represent respectively the mean anomalies of the Sun and Moon, 
and the mean angular distances of these bodies from the node of the Moon's orbit on 
the ecliptic. 
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These quantities are used for convenience instead of the /, /', F, and D of Delau- 
NAY, with which they are connected by the relations 

D=A-A' 

The required formulae are presented in the form I have found most convenient for 
manipulation in the following pages Instead of writing in full each expression such as 



^'Z('^+''^+^^+^'^') 



I use the plan of the German astronomers, simply writing the indices i, i',i, /, 
opposite the coefficient C of the term. 

We commence with the value of sin P or - given by Delaunay. If we put 



we have 



r 



^=i-c5+c5^-6^+etc. 
a 



The several terms of , so far as is necessary for our present purpose, are given in 

the first of the following series of tables. Only those terms are given which are necessary 
for determining A:,, A„ and k^ as high as m^ and e'^ w*, and as high as the sixth dimen- 
sion in all the other small quantities, as well as their combinations. Now, the mode of 
forming the quantities sought is such that it is not necessary to carry any coefficient 
to a greater degree of approximation than is necessary for obtaining its square to the 

necessary order of accuracy, with the single exception of the coefficient of cos g* in -. 

All the coefficients exceeding the third order are therefore omitted, and those retained 
are only given to such a degree ofapproximation that the sum of the lowest and high- 
est orders of the terms of each separate coefficient shall be 6 or 7, according to circum- 
stances. This remark applies to all the following expressions for the Moon's coordi- 
nates. 

The several powers of - — i to the sixth, inclusive, being formed, they are combined 

according to the preceding formulae to obtain the value of -, the coefficients for which 

(I 

are given in the next table. 

The value of )ff, or the Moon's latitude, as given by Dllaunay, follows. The 

several powers of yff being formed, sin /3 and cos /3 are obtained by the formulae 

sin /?=/?-! /?«+-i- /?^ 
o 120 

cos /?= I -i /?^+ -/?*-— >ff*, 

2 24 720 ' 

for which the expressions are given. 



154 



THEORY OF THE INEQUALITIES IN THE MOON'S MOTION. 



These quantities are then multiplied by -, and the products are given in the next 

two tables. The expressions for cos d v and sin S v are formed in the same way from 

S V, and their products by - cos /S are given. We thus have the coefficients of cos / 

and sin / in the expressions (48) for ^ and 7, the rectangular coordinates of the Moon 
referred to the ecliptic. 



cos 



^■+ 



0000 



■^ /»H - e'* >«*- i|i m*-^l »»»+ - ' - e*m'-{- A ef* m^-\- 2 e" y^ m^ 
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COS 

9i /» ^1 ^' 
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r 
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From this we obtain 



cos 

9^ 9\ \ ^' 
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COS 
9^ 9\ \ ^ 



- (continued). 
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We have from Delaunay 



157 



sin 
-2010 
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From this we obtain 
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sin /3z=, 

2 

— 2eY—\eY^+^^Y+^eynt^ 
4 32 

2 y — 2 c* Y—y^+ — e* y +e^ y^— - r*— — y^ ^ 
32 4 32 

2ey—^€? y—ey^ eym^ 

2 2 



$eym 
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^eym 
4 
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C08)ff=I + 
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128 64 32 2048 
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COS 
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COS >ff (continued). 
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- sin pzz 
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sm 
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- sin )ff (continued). 
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COS 
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We have also from Delaunay 
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Sin 
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6 V (continued). 
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From this we derive 
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COS S V (continued). 
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Multiplying by - cos ^ we find 
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77 e^^^^^SSe^m' 
16 16 



15 / 



16 48 



0002 



— 9.y^fn 

4 



O I I— I 



e' a 



O O I— I 



15 
— ^ am 



ThiB values of A;., A,, etc., may be formed from these expressions without going 
through the process^of forming the rectangular coordinates x, y, Zj in full. Let us 

consider any term in the value of- cos /3 cos.^t;, the corresponding term in the value 

of - COS 13 sin S v, and any term whatever in sin /?. Let these terms be 
a a 



cos 13 cos d v-zzk cos N 



cos >ff sin S vzz. hf sin N 
sin >ff =c sin W 



k and kf representing the several coefficients given on the preceding pages. If we 
substitute these values in the second members of (48) p. 129, we find 



a 2 



=i (k+k') cos (/+N) + i (k-}(f) cos (/-N) 



?=- (k+k') sin (/+N)+i (k-k") sin (Z-N) 

^ =: c sin N' 
a 

The general expressions for the rectangular coordinates of the Moon, referred to 
any system of axes having their origin at the center of the Earth are 
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the direction cosines a, b, . . . having the values given in (51) p. 131. If we 
put 



2 



h=\(k-1d) 



(78) 



we find 



■=*i(i-x'Ocos(0'-r+^+N) 

+^(i -/^) cos (9' _T+?-N) 

+A:,/^co8(^ + r-/-N) (79) 

+Ai/^co8(9' + r-Z + M) 

-cyV "T^cos ((9' + N') 



+c/ V I -x" cos (©'-N') 

?'=fc,(i-/^)8in(9'-r+/+N) 

+ ^(i_/2)sin(9'-r + /-N) 
+A;,/''8in(0' + r-Z-N) 
+ ifc2/*8in(0' + r-/+N) 
— c / V 7^^* sin (^ + N') 

+c / VF^'" sin ((9'— N') 

-=2 Ai/ vr^/^sinG— r+N) 

+ 2 A:, ;/' Vi— /" sin (i— r — N) 
+ c(i-2/^8inN' 
Representing the coefficients oi g,g', A, A' in N by*, i',y,/, respectively, we have 

N -i g+i' /+> A+/ A' 

(80) 
W=Hg+i\g'-\-j,X+j\X' 

Substituting these values of N and N', after effecting the substitutions (53) to 
reduce the expressions to the form .(54), remembering that 

g zzX — CO =:« — tt 

g'zzX' — co' = e' — 7r' 
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the values of x and e become 

?=*, (I -O COS {(I +i+j) e-i n-ij+j') ©+(*'+/) e'-i' 7r'\ 
+ A, (i -/O cos {(i -i-j) e-\-i rr + (j+j') e-(i'+f) «'+i' rr'] ' 
+A^/*co8 {-(i +i+j) £+i 7r + (j+j') e-(i'+j') e'-\-i' n' + 2&\ 
+A^/*co8 \-{i-i-j) e-i 7r-(j+j') 0+(i'+j') e' -i' ;r' + 2 e'l 
-cy' vr^^cos {(*,+.;,) €-i, ;r-(j,+j\) 9+(i\+j',) e'-i\ a-' + ©'| (8i) 
+ cr' Vr^^cos {-(i.+ii) « + «! ^+Oi+/i) '9-(»',+/,) e'+t'i T' + e^i ■ 
-=2k,y' vr^^sin Ki +i+» e-i 7r-(j+j') e+(i'+j') e' -i' n'-ef\ 

Or 

+ c(i-2 /^)8in {(ii+iO *-^i ^-Oi+/i) ^+(i'i+/i) ^'-i'l ^'\ 
It is unnecessary to write the value of y, as it is the same with that of Xj except 
that cos is replaced by sin. The entire values of x^ y, and^ will consist of the sum of 
an infinite series of terms like the above, formed by the successive substitution of each 
of the different values of N in the general expressions (79), each term of N giving rise 
to four terms in x, and each term of W to two. Should the argument be the same in 
any two of the terms thus originating they must be combined into one. But it is 
easy to see that this can not be the case. For, as already shown (§ 8 p. 137), we 
must then have the coefficients of e, tTj <9, e\ 7r\ and 6^ the same in both terms. Now, 
let us represent any one system of values of t, i\ j, f by 

h ^ yJiJ ^^^u '^2> *3» ^4> 

and any other system by 

^» * ^J^J ^^'^11 hi Hy *i' 

In order that the substitution of these two systems of values in (81) should give 
the same coefficients of e, ;r, etc., in any two terms we must have 

S3 + S^ = dtz(k + tA) 

the upper signs alone or the lower signs alone being used. 
But these* equations give 

Si^zizh'j ^a^^it^J ^s^^i^J ^4^^ it ^4 
and consequently 

But, all the pairs of terms in which the values of N were the same, with opposite 
signs only, have been combined into one as we went along, so that there are none 
fulfilling this last condition, as we may easily see by simple inspection of the indices 
corresponding to the preceding values of k and k\ 
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Let lis represent that part of A;,, A,, etc., formed from the expressions of the coordi- 
nates of the Moon by 

a;^ A.; (^f^K etc. (82) 

Let us also put 

^, the coefficient of the time in t g-\-i' (f +j' A +/ A', or in N 
^1 that in N' 

so that the coefficients of the time in the first four terms of x and y, formerly repre- 
sented by fe, are now it ^ it P- Referring to the formulae (62) we see that each value 
of N and of N' in the preceding tables gives rise to the following terms in A, 

(i+i+»(»+l>)A^(i-/7 
+ (i-i-j)(«-J»)*l(i-y'*)* 

+ (i-i-;)(n-|,)A|/* 

+ 2(i+»+»(n4-p)A?/*(i-/') 
+ 2{i-i-j) (n-p)Ai/«(i-0 

This expression reduces to 

(, +e+j) (n+p) *?+ (I -i-i) (n-p) *i+ J (i,+j,)px c» 

The values of h,, A„ etc., are fonned in the same way by replacing i -\-i-\-j, 
I — i— 7, etc., by the corresponding coefficients of a; 6, etc. 

Completing the expressions by the sign 2, which refers to the different systems 
of values of i, i',j, f, we have. 



A.=^|(i+t+i) (n+p)ki+(i-i-j) (n-p)li+L(i,+j,)p,<^^ 
K=2\ —i(n+p) 14+i (n—p) Ai— ^ ti/>, A 

h,=s\-(j+j') (n+p) A^+ 0+/) (n-p) At- -^ (.n+j\)p, c»} (82) 

h:=2i^(i'+j'){n+p)l4-(if+j')(n-p)li+l(i\+j\)p,<j 
h,'=2\—i' (n+p)ici—i' (n—p)l4— -' i\p, (^ 
h:=2\- 2 {n+p) k\- 2 (n -p) A| i /" 
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We now give the expressions of A?, AJ, (?^ and />, corresponding to each term of 

the expressions for - cos P cos rf v, ^ cos P sin 6 v, and r sin /S^ except the first, or con- 
a a 

stant term of - cos /S cos S Vy to which the above expression does not apply. If we 

consider this term alone, and call it Aq we shall have 

h,z=:nI4 

while all the other values of A will vanish. The formation of h^, A:J, AJ, is too simple to 
render any explanation necessary. To form the values of p, it is necessary to have 
the expressions for the secular motion of the perigee and node of the Moon, or of its 
mean anomaly and argument of latitude, under the action of the Sun, in terms of the 
elements of the Sun and Moon. These are obtained from Tome II, p. 237, of Delau- 
nay's Theory, -by making the change of elements indicated on p. 800, 

We have 
dg 



dt \ 4 32 ^8 8 ^ ^ 128 



dg'_ 
dt 



mn 



(83) 



dt 



1 4 32 2 ^8 2^ 128 ) 



Thua arise the following values of the quantities which enter into the equa- 
tions (82). 



*, i, J, J 
0000 



3 32 ' 192 96 

— 5 e" /«='+ ' 591 ,„4 _ 3315 ,,u „,s — 57^* ,«» + 347 ,„.. _ J_ ^t 
1 152 64 64 96 144 



y*m' 



.425 
32 



e'*m^ 



32 ^ ' 768 ^96 

I 2567 , 2 2 1679 2 /» 2 , 1585 ,2 2 2 6053899 , . 

96 64 32 18432 

_69_43? =,„,.^LL56i9,.2^.^3.9365^._2J^,a^ 

9216 768 10368 16 

_3_^«2 ,„2_^. 63.651 ^.,„._262 ^, 

256 64 32 



p=o 



174 

1 o o o 



2 O O O 



3000 



— I 100 



0100 
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482 ^4 , 25^ '9^ ^ 

2 16 384 8 9216 

42 o 256 4 120 

128 2 1024 






32 



Af =^?- e*- ?- e'— ^ e* /+ ^ e* m^ 



64 32 32 



128 



A4= j- e* + -L e«+^ e* y^ |_e« m'' 



"64 96 32 



384 



2-^w'') 



9 

576 
i>=3» 

256 



256 



^r= — n 



81 



4 4 16 4 16 



.^iile"'m*-i^^Me"'-* 



32 



256 



jK=? e'* »«*4-5 e'* m'+^ ^ e'* m"-?? «'« y^ tn*+^-i e'* m^ 
^4 ^4 16 4 ' ' 16 

221 1 ,, 4 140 877 ,, . 

e nr — -i-^ em 

32 256 

pzzmn 
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», *, Ji J 



O Q 



256 



256 



pz=in 



0200 



64 



64 



p=:2 mn 



— I 020 



10 



ii=%^y 



0020 






Ai=z/~?eV*-^/m' 



12 



/>z=w(2 + ^m2) 



1020 



A;if=:o 



4 



— 2 O 2 — 2 



1024 



— I O 2 — 2 



)K=5625,.^2 
■ 1024 

;)=o 

ifc?="fe''»?^+-5°5_5e2^3_225^4^2_675g2 2^2__i 125^^,2^, 
256 512 256 128 256 

+ 553 743^2^4 
8192 



176 
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^ ^ J J 
-^I O 2—2 



;o 2025 2 2 , 33075 2 3 202s 4 2 5265 2 2 2 

256 512 256 128 
256 8192 



O O 2 — 2 



I O C— 2 



— I— I 2—2 



O— I 2 — 2 



— I I 2—2 



7" 9 4. 135 2 3 . 21 6 , 2025 4 2 . 5085 o 4 27 2 4 

256 256 64 1024 1024 128 
.45 ^/2 ,^4_^447 .^,j«_93 ^'2 ,^5^ 83 ^^^ 



256 



256 32 



12 



256 256 32 192 1024 32 



256 



1024 



64 



16' 



779 /2 6 . 106 513 7 
1 2 864 

p=:n(2 — 2 w) 

64 



^-256 

i)=3w 



Jfc?='^5g2g/*^2 
256 

256 



^_44Lg'2^4^ l8^g/2^,» 

1024 250 

1024 16 

p=n{2 — zm) 

256 



A^= 



256 



p-=n 



768 



■ »» 
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t t J J 
O I 2—2 



1024 256 



1024 



584 



pz=,n{2—ni) 



0002 



4 



o o I — I 



16 



p=o 

A:2--A2_5.a*w2 
1024 

1024 
pzzn 



o I I — I I A?zz:^5e'2a2 

64 



Hzr^e'^a^ 



64 



^niw 



2010 c^z=.e^Y^ 



Pi=—n 



4 32 



0010 



2 



325 



24 



12 



72 



y" W 



«,=nri+^m2-9 m3^3g2,,,2_^9^/2^2_3 2^2_273^4\ 
^ ^4 32 2 8 2 128 ^ 



A. P.— VOL. V, PT. Ill- 



^8 
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i, i', J, f I 



2 6 



2 O I O 



Pi=2n 



^,2_. ? 4 2 

I i6 



i>i=3 w 



o— I I o . c ^ ^^e y m 



Pi=n 



Olio (rzz^e'^y^fn^ 

lo 



jh=n 



—I o 3-2 



64 
.pi=z2n 



e^ y^ w* 



O— I 1 — 2 (r^zz^-^ r/2 y2^2 

16 



— I o 1 — 2 rzz- c^y^w 

i ^ 

I i>,=io 



o o 1 — 2 r^iz: \ y^ 'm^+ -^^ y^ m^'\-^- y* m^+ ^ f'^ y^ w^— "^^ ^'^ r^ w^ 

16 32 16 32 16 

. 8655 ., , 
^ ^ y^ m^ 



4 ^ 



I o 1 — 2 r^=^^ e^Y^m^ 

64 

I |?i:=2W 

o I 1 — 2 c^-=z-^^e'^y^m^ 

16 



' i>i=z-i 
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From these expressions we form the values of -^^ A;? and — ^A|, which, when 

multiplied by n enter into the expressions (82) for A., A„ etc. The results are shown 
in the following table. On the left margin are given the indices i,^, i',/, as heretofore; 

then follow the corresponding expressions for ' ^ A:? and ~^ A| which are given in 

the middle of the page. For brevity the former is represented by the symbol (i) and 
the latter by (2). On the right margin are given the coefficients by which these 
quantities are multiplied in the expressions for A„ h^, A^, etc. To save space, each set 
of multipliers is arranged in two lines, and the several sets are separated by a dash. 
Each term being multiplied by the corresponding coefficient, and the simis taken^ 

and, in the case of A., the value of A:? being also included, we have the values of— ^, 



n 



-. etc. 
n 



^ ^y Jj J 



1000 



2000 



3 O O. O 



— I 



O O 



Terms arising from the Longitude. 



^ "^ 2 4 48 64 96 

+ - e' x'+3 e' y'-^-^ e' m^- ilii ^ y^m" 
4 4 192 

-305^2^/2^^2^9359 ^2 ^4 



32 



1152 



(2) = ^1 (? m^- ^1 e'm^+^^-^^ e^ m^-i^ ^2 2 ^2 
^ ^ 16 32 128 8 

32 512 

(i)zzf7e*-?7g6_27^4 2^9 ^. 
64 32 32 32 



'm^ 



(2) = - ' e*-- Le«-9,4 2 

64 96 32 192 



7_c*w* 



0): 



:4.« 



«=-.88'' 



(0=0 

(2)= '""^^ ^2^/2^2 
^ ^ 128 



Coefficients 
for 
£, ;r, ©, 

+ 2 — 1 O 

o 0—2 



O+I O 

o 0—2 



+ 3- 


-2 








0- 


-2 


-1 + 2 








0- 


-2 


+ 4- 


-3 








0- 


-2 


— 2 + 3 








0- 


-2 



+ 2 — 1 o 

— I + I — 2 



i8o 
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Terms arising from the Longitude — Continued. 



*, », J, 3 



o I o o I (i)=5e'^«2+l|e^e'*»i*-99e'2^2,„2+8i^„^2 
^ "^ 4 i6 4 16 

2283 ,, A 1 18 oa'^ /, 5 

^ e ^ m* — --^ e »r 

256 



32 



81 



I I o o 



32 

(2)=0 

(2)=|ie-m^ 



1020 (i)=^e^ y^ 

o 



^ ^ 4 16 4 16 

2283 ,2 4 12^ 189 ,0 5 



256 



'm' 



0200 



0020 



1020 



(2)=0 
(0=0 

(2)=-y'-\-^e^y*+l^r*m^ 
4 12 

(0=0 

(2)=-^-e*y* • 



— 2 o 2 — 2 



(0= 



225 4 2 



1024 



(2) = ^-^e'm' 



Coefficients 
for 

+ 1 00 
+ 1 — I — 2 



+ 1 00 
— 1 + 1 — 2 



+ 2—1 o 
+ I — I — 2 



+ 1 00 

+ 2 — 2—2 

+ 1 00 

— 2 + 2 — 2 



+ 2+1-2 

O O — 2 



-1 0+2 

o 0—2 



■2+1+2 
O O — 2 



+ 1 + 2 O 
— 2 O — 2 



1024 



+ 1—2 O 
+ 2 0—2 



h ^, Jy J 



-I O 2—2 



O O 2 — 2 



I O 2—2 



— I 1 2—2 



O — I 2 — 2 
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Terms arising from the Longitude — Continued. 



^ ^ 128 256 128 64 

^iL25,2,.2,„2^475 563,2,^4 

1 28 4096 

, N 2025 , 3 , 126 225 , 4 

(2)= ^e^nr+ — ^rm* 
^ ^ 128 1024 



27 



405 



( I ) = JL/ „,* + ^P f 2 m' + '--L „,» + 



117, 



256 256 



28 



6075 
1024 



e* wi* 



+ 



14 175 2 4 81 2 4 135 '2 4 1 1173 6 



1024 



128 



256 



256 



m" 



IO7I ,2 6 1 6627 7 

128 ^ 384 



(2)=-3^^«^+-^^5 ,2 „,.+57 2„,3_26o3 ,,^, 
256 256 32 384 

-.225 e*„,2_45 e^y2„i2_ 9 ^4,„2_^9447g2„,4 
1024 64 16 1024 

+ 365 2^4^ 1805 ^,^„,4_i5217„,. 
32 256 768 

384 3456 
^ "^ 16 

(2) = -i3.^«.* 



^ ' 128 

(2)=0 

^ 1024 1024. 



^W'^ 



>. V ) 7 680 ,0 4 142 700 ,2 5 

(2):=.— ' — v" w* — ^ ' -^e^m^ 

^ ^ 1024 1024 



181 



Coefficients 
for 

+.2+1 O 
— 2 0—2 



O— I o 

+ 2 0—2 



+ 300 
— 2 0—2 



— I GO 
+ 2 0—2 



+4-1 o 

— 2 O— 2 

— 2+1 O 
+ 2 0-^2 



+ 2+1 O 
— 3+1 — 2 



+ 300 

— 3+1—2 



— I 00 
+ 3 — 1—2 



l{{2 
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Terms arising frmn the Longitude — ^Continued. 



^ *, ^, ; 



— I I * 2 — 2 



O I 2 — 2 



O O O 2 



O O I — I 



O I I— I 



128 



=:0 



--^-e''m'+—^-e'^m' 
1024 1024 



3072 



1024 



e'^ nr* 



zn^yUn' 



16' 

225 2 2 
512 

rzo 

32 
=0 



Coefficients 
for 

+ 2 + 1 O 
— I —I —2 



1 

+ 3 





— I - 


- I —2 


— I 





+ 1 + 1-2 


+ 1 


0—2 


1 +2 


0—2 


+ 1 
— 2 


+ 2 
*o— 2 


+ 2 





— I 


0—2 


+ 2 


j 


0- 


-1 — 2 ' 



h i\ i, / 



— 2 O I O 



Terms arising from the Latittide, 

or values of -^ c^ 
2 n 



--e'y' 
2 



* 27 2 2 2 

— I o I o -Le^y^m^ 



Coefficients 
for 

f, ;r,. 6 

^^ ^', & 

— 1 + 2 — I 

000 



O+I — I 

000 
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Terms arising from the Latitude^ 
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or values of - ^ c^ — Continued. 
2 » 



^> «, 3, J 



0010 



2 >'^ — 2 c^^— 2 y*+ 1 y^ m^ + \ y^ m^—^- i'^ y^ 
6 16 16 

24 4 48 



I O I o 



^2 y2 _ 3 g4 j/2_g2 ^4 _^ g2 ^2 ^^! 
2 6 



27 4 2 

2010 ' ^* y^ 

! 32 



O— I 1 O -9_ e'2 y2 ^2 



32 



O I I O i -9 6'2j/2^2 



32 



I o 3—2 ^^^e^y^m^ 

64 



2 7 

O O 3—2 —^1/2^* 

I 128 



o— I 1—2 j 12 e'^ 1/2 ^2 
32 



—J O 1—2 



O O 1 — 2 



9 ,,2^2_, 87^2^3 , _9_.2,,2^2_45 / 



32 



J/2 W-2+ ^^- y2 ^3 + ^ ^2 ^,2 ^2_ -to ^,2 ^2 ^1 



64 



64 



32 



+ 9.^,4^2^4503 2 ^^^4 
32 1024 



I O 1 — 2 



O I 1 — 2 



64 

^-e'^y^m^ 
32 



Coefficients 
for 
f, r 

+ 1 o— I 
000 



+2- 




-I- 



- I 



+ 3- 



— 2- 



- 1 



+ 1 
— I- 


0- 

4-1 


- 1 



+ 1 
+ 1- 


0- 
— I 


- 1 



+ 2- 
— 2 


■fi- 



-I 



+ 3 
— 2 


0- 



-I 



+ 1 


o+i 
+ 1 




— 2 


+ 1 + 1 




+ 1 O+I 
— 2 00 



+2-I+I 

— 2 00 

+ 1 O+I 

— I — I O 
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We have now all the data for forming, by means of the equations (82), the final 
valufes of A„ h,, etc., so far as they depend on the coordinates x, y, z of the Moon, which 
parts are represented below by adding the symbol 5- We shall, however, for conve- 
nience in our subsequent work, change the form of the common factor //g <^ ^ l)y which 
all these quantities are multiplied in order to obtain A;., A:,, etc. We recall that, in 
Delaunay's theory, the quantities a and n are connected by the relation 



Then, since we have put 



a*w^zrw2+^ 



_ m^ ma 

A*2 — - 



we have 



t>*2 + W8 



an 



which we shall adopt for the common factor in (juestion. 

We now find from the preceding tables for the terms in question. 

aw 1^38 2 2' 32 16 8 44 

— y* — ;^a^^w^ 

o 



+(-f^-l'") 



w* 



-< 895 23255 ^, 3995. /2_ 1585 y2\„« 
^^288 512 ^ 192" 128 ^ ^ 

, /1 261 , II ?49 ,2x 5 , 3028 . , 76741 7I 

, ^945 2 3,2 CK)i 623 2 J 
256 36864 J 

. _»V^3 j_ 2 2^^2 2_35 ^2 j^^43 ^i ^2 

a« 1 4 ' ' 32 

+ ( _53 ,,2_2ioi^2 ^2_45 ^/2 ^2^ H.3 4A 
' ^ 96' 192 32 32 ^ 

+ 75 2„,,^i6oo7 2^J 
64 9216 J 



w^ 
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•" an [V 64 ^16 •^ ^ 32 ^ J 

an [ 3 

We have next to consider the corresponding quantities formed from the solar 
coordinates, x\y\ /. Owing to the minutenesss of the perturbations of these coordi- 
nates in the problem of three bodies, the largest terms of this kind will be those aris- 
ing from the values of a/yt/j and / in the theory of elliptic motion. Into this theory 
the lunar elements g^ tt, and 9 do not enter; it will not therefore add any terms to A:., k^ 
or ke. However from the form of the canonical equations (64) it appears that A:/, A;/, and 
k/ are the conjugate eleinents of «', 7r\ and ©', that is to say, functions of f', ^r', and & 
satisfying the conditions 

while all other combinations of the same kind vanish. Now in the theory of elliptic 
motion the values of these conjugate functions are known to be 

V0=/k.(i-O* (86) 

K^=kJ cos i'=kj {1—2 /^ 

which are therefore the principal terms in question. 

We have next to consider the order of magnitude of the terms which will be 
introduced into the values of a/, y', and / by the perturbations of these quantities due 
to the separation of the Earth and Moon. From the process of integration developed 
i^ §§ 3> 5» ^^^ 6 when .Qi is taken account of, it follows that any term of x' in (60) 
which contains either f , r, or 9 in its argument will contain the factor 

and -will be of the form 

rc'rzAw^a' cos (I f+P ;r+ . . .) 

s being a numerical factor of the order of magnitude unity or less. The principal 
term of A;, thus arising will be by (62) 

and the complete expression for A;, will contain the factor >Ui n A;^ m^ a'^ the ratio of 
which to the common factor ja^ ^ ^ of \^ is of the evanescent order of magnitude ir^ h 
These terms are, therefore, entirely unimportant, ^nd the complete values of A;,, A;,, and 
k^ may be taken as equal to A;, j, , etc. 

The terms of ocf which depend on the arguments e\ n* ^ & exclusively will receive 
increments of the order of magnitude of the terms in f, ;r, and 0; that is to say, the 
terms in question will be of the form 
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x! z=.a' (s' -{-k m s) cos (i^ f' +h ^' + h ^) 

5' and s being numerical factors. Each tenn of this kind will give rise in A/^to the 
term 

ria'^n' (s'+kmsy 

The principal term ij a'^ n' ,9'^ of this expression will be contained in the value of kj^ 
already ^iven in (86). The principal part of the remaining terms will have as a factor 
the factor of the former terms, yuj a'^ n\ multiplied by m k. This gives rise to the 
common factor 

a n 

in the terms of A:/, kj ^ and kj. The terms introduced by the coordinates x^ y^ and, z, as 
found in (85), contain this same factor. We conclude, therefore, that this portion of A:/ 
etc., must remain undetermined until the solution of the problem of three bodies is 
complete. 

§.2. ^ 

FOKMATION OF THE DIFFERENTIAL EQUATIONS. 

We have now to differentiate the preceding values of k„ k^ etc., with respect to 
the elements of which they are functions. To find the derivatives of the powers of m 
with respect to the logaritlim of a, which I call », we remark that its value is 

^^ _f^ __ I Wi + »^2 + W3 a' P 






n I ^2+^3 
From this we derive the general rule that if 

u z= a* \ C(^ + (^1 m + ('2 m^ + ^3 m^ + etc. } 
then 

^-^-=aN^Co+(i+-)^lW + («+2 )^2Wi' + etc.j 

d u i\ % 6 o < 



(87) 



- , , . , „■ ^2 nr — <Cz w — etc. 

a 70 I 2 2 2 



We thus find, dropping the symbols }) and in (84) and (86), 

(lk,_m2ntJi 7 21 21 7 7 4 , 35 2 2 7 « 63 63 

Clio- arr\2+^-6-v6' + ^y-^' -6/+^'-y-2y-'i2^' +^' y 



^'-W-] 
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a e an[^4 8^2'^ 8 y 16 256 J 

dy an [^ *^ ^ ^2 ^ ^' '2 ^ ^ 4^ 64 I * 

ae a w [^ 8 ^ 2^^ 2 ^ 96 32 I 

^* m,m,j_i I . s_^ 5 . . -«I97.,,_'4 483^._595,.,. 

a» o« ( 4 16 32 32 -^ ^16 "^ ^^ 768 768 16 

216^ , ,,. , 19725 , , 26021000 o J 

256 ' ^ 256 ^ 73 728 I 

dk, »M2»»3| I J , 5 i 53.. 3 5 , /i'7i 2069 85 , 309 „. 2 

de an \ 2 ^4 -^ 4 -^ 8 ' ^ 192 96 4 • 64 -^ 



, ■?94'> s I 2 001 623 4 
^ 128 ^ 18432 



dy a n \2 ' 8 ' 4 



(ffc rWaWaJ 309^2 ^,^a[ 
rf e' ~ a n I 64 • I 

(fA;, m^nh] , , 1.2.,2 35 o. 43.„2,/' 371., H 707 ;$ , , J.ooi . 



df» aw(''2^ 8^ 64 '^'^v 192' 384 "^ ' 64 

64 ^ / ^ 64 '^ ^ 18432 ) 



dk0_m2 Wg 



rfe aw 



28 06 I 



di- WJ2^^3( .0 9 , 43 A ^ S'? 2101 „ 14^ s 

J-/-^\H>'+^^>'-35''V'+t|«V+(-fiK--^-e^r + -|^y' 

_4i,.^N^.+ 7J , .i6o^^,„. 

16 '^ ^^ 32 '^ ' 4608 

rf&,_m2;/i3J_45.,2g/^2l 
rf e' a n ] 16 I 

The derivatives with respect to »' may be found by the condition 

dk , dk I , 
a» dv 2 

A rigorous solution of the problem will now be obtained by taking the deriva- 
tives of A/, kjj and kg, with respect to the six elements which they contain, substituting 
all the derivatives in general equations of the form (64), and then solving these equa- 
tions with respect to - -, - , etc. But we have seen that in the values of A,', */» etc., 
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the lunar elements are multiplied by factors so minute that they may be wholly neg- 
lected. The perturbations of the solar elements may therefore be obtained from the 
values of A/, hj and V, and the results substituted in the expressions for the lunar ele- 
ments. The process will be as follows: 
From the equations 

dio' _ dv' dR dio' d'R d v\d R 
d f "■ dkj d e' '^dkj dTT'^dV, d 6' 

(88) 
de^ _d^ rf R de' dTR d_e' dU 
d t -dk\ d t' +rf^-'. d7r''^d k\ d & 

we are to derive the perturbations of »' and e' . Those of y' are at present omitted, 
because this quantity does not enter into A;„ A:,, or A^, so far as we have written them. 
From the equations 

dj'_,_dv'd^ d^dB^ rf/rfR 
~dt ~^ Ik: dio'~dk:~d€''^dk: dy' 

dTt' _ _ dio' d^jU rfR dy' rfR^ 
dt ~ d kj dio' ~dkj d e ~ dk; d y' 

where 

w' rz { Wj + ^^^2 + ^^3 1 * ^~^ — (^h + 'W'i! + ^Wa)* C""^"' 

c being the Neperien base, we are to derive e' and Tr\ or their perturbations. 

Putting the first three of the equations (64) in the form 

dk^d X? .dk^ de dkg d y_dR_dk, dv' _dk^ d e' 
dio dt de dt d y dt de dio' dt ' de' dt 

dk„ d X? dk„ de^dk^ d y dR dk^dv' dk^de' ,^ . 

Tio ~di^Te dl'^Ty U^Ttt rfV "rfi ~de' "di ^ ^^ 

d ke d X? ^d ke d e ^d k^ d y __dR __d k^ d v^ dk^ de' 

dl^ U'^de (ft'^'dy Tt^dO d'l^' ~dt 17 Tt 

We are to solve these equations with respect to - ■, -, and -VV- Let us put 

dt dt dt ^ 

^, the time introduced through the varying solar elements, so that 

dk^_dk, dv' dk,d e' 
dt'~~Tv' 'dt^^d? Ti 

etc. etc. etc. 
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We shall then have 



dt~^dkj ' \de dif) 

, /rf^'N fd/R_dk^\ 
'^^dkj'\d7r dif) 

x(dx>\ rd^dkjx 
'^\dk\) ' \d~9 dt'J 

de _/^de\ /'dR_dk,\ 
dt~\dk\J ' ^de (ft J 

xfde\ /'dU_dk^\ 
'^^dkj ' \dO dr) 

dy_ — (dY\ fdU_dk,\ 
di^\dkj ' Vrf7 dlO 

x(dy\ /'dR_dk,\ 

'^\dkj ' ^dTT It J 

i(dy\ fdR_dk^\ 



The derivatives ( — \ etc., are here inclosed in parentheses to indicate that they 

are to be taken subject to the condition that »' and e' shall remain constant; in other 
words, the derivatives with respect to A, are to be formed by solving the equations 

d k, /'d X? ^ d k, / de x d k, /d v \ 

djo VrfT./"^r/T KdlJ'^dr ydkj~ 

dv\dkJ^de ydkj^dr^dkj ^^ ^ 

dk^ / d 10 \ dkg /'de \.dk^ rd y \__ 
dv\7lX)'^de ydkJ^dyVdkJ"^ 

and those with respect to k„ and k^ by solving similar sets in which those quantities 
enter. Having found the perturbations of », e, and x? those of f, ;r, and 9 are to be 
obtained from the equations 
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dk, (I fo dk„ d tTq d k^ dO^ rfR 

d 70 dt d X) ^d t d 70 dt d 70 

dk^d e^.dk^ d tTq dke dO^ ciR 

de dt ~de d t de dt~ de 



dk^ d €q dk^ d tTq dk^ ddQ__dR 
d Y dt d Y dt d y dt ~ d Y 



which give 



rf_^_a 6^ g _^ _fd7o\ dli_^ de\ dH _/d Y\d\i 
dt" dt^ '~ ' \Jk,) dx7 \JkJ~d^ \dkJ~dY 

The values of n, tt^^ and 6^ are the expressions for the mean motion and the secu- 
lar motions of the perigee and the node in the theory of three bodies, which last are 
taken from Delaunay {Comptes Bendus^ 1872, I, VoL 74, p. 19), as follows: 

7r^=M^n\^—le'-6Y^+%e'''- ^^e'+^^e^Y'-'^^-y'-gv^e'^- 9-e^e''+'^-^e'' 
14 8 8 32 8 4 16 32 

/225_675,2_l89 2+82^ ,.2+81^, ,ip^ 

V32 64 8 32 32 16 4 64 y 

/407i_3i 6o5^2_3963 2 . 6ij79^,2\ 2 

V 128 512 32 256 y 

. p65 49 3_i 483 665 ^2_335 403 2 . I 7 67849 ^/2^^^3 

\ 2048 4696 512 1024 y 

12 822 631 4 . i_273 925 965^5 .71 028 685 589 ^J 
^ 24576 ^ 589824 ^ 7077888 ""M 

I 4 2 2 8 04 o- 4 4 32 

'^V32 32 16 32 256 16 16 16 16 y 
^Vi28 128 128' ^ 256 y 

+r9797_i6j4ij,._7i_85 ..73 423,.2V„3_^i99^73^*+1657^73^^ 

V2O48 2048 1024 '^ 1024 / 24576 589824 
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The perturbations of », e, and y, which are of the first order as to the disturbing 

forces, and found by integrating (90) ^as if the elements were constant in the second 

members, maybe divided into two classes: (i) Those which arise from the derivatives 

d k 
of R with respect to f, tt, and 0, and (2) those which arise from -=-^ , etc. Let us consider 

the latter. They are given by the integration of the equations 

dv ^_/"dv\dkg /dv \dk^ /d v\ dk^ 

'dt~^^dkjdf ^dkj dr ^dk,Jdf 

d e___/'de\dk,_/^d e\dk,_/ de\d k^ ^ v 

~dt^ ^dkjdif Vd'kJ'df ^dkjdt . '^^^^ 

(ly ^_{fly\dk,^fdY\dk^ _rdy\dke 
df^ \dkjdt' ^dkJdT ^JkJJr 



The corresponding perturbations, 5 », 6 e^ and 6 /, to quantities of the first order 
will be given by solving the equations 



dx? de dy dv . de 



V::"'+M"+Vy''+^"''+'^'''=° ('*) 



dv de dy dv de 

The values of S v' and 6 e' are regarded as known in advance. This curious 
theorem may embody some principle applicable to the disturbed motion of three 
bodies, which has not yet been fully mastered. The quantities I have called A„ A;,, 
anrf Ai, are closely related to Delaunay's L, G, and H, {is is shown in the next 
section. The theorem shows that the secular ac(*eleration of the mean motion and 
of the motion of the perigee and nodes may be found by the equations 

(JLzzO; (JGzzO; (J HzzO; 

by making the transformation on p. 800 of Delaunay's Theorie, Vol. II. This result 
I have verified so far as the mean motion is concerned. 
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The numerical values of the derivativesof A., A;,, and A^froin the formulae of pp. 1 86-87 
are next to be .computed. To exhibit the uiicentainty to which many of the derivatives 
are subject, owing to the approximation in powers of m, /, e, and e! not being carried 
far enough, the results are arranged by powers of w, the quantity whose powers con- 
verge least rapidly. From the last two or three terms the probable law of forma- 
tion is roughly guessed at by induction, and a correction for the siim of the higher 
terms thus concluded. This correction must always be regarded as uncertain by a not 
inconsiderable fraction of its amount, and several of the results are therefore quite 
uncertain. 

Formation of -r— ^ 

Terms independent of /w + .50ocxx> 

multiplied by m^ — .006 492 

w^ — .000050 

m* + .000 602 

w^ -|- .000 246 

w^ f .000063 

ni} + .000012 

Higher terms by induction 2 

Sum -f 0.494 383 

Formation of . - 

^2 m^ e d e 

Terms multiplied by m^ — .0042 

m^ —.0058 

m^ —.0028 

Higher terms by induction — .0010 

Sum —.0138 

Tjy ,. J. an 1 dk, 

FoTfnation of =— ^ 

m^m^ y dy 

Terms multiplied by m^ -f.oi68 

m^ —.0009 

m* —.0008 

Higher terms by induction — .0004 

Sum . +.0147 
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^93 



Formation of 



an dk. 



W2W3 d x)' 

Terms multiplied by w^ +-005565 

m^ +.000045 

m* —.000556 

m^ —.000231 

w* —.000059 

m^ . : —.000012 

Higher terms by induction — .000 003 



Sura + .004 749 

Formaiion of — ~ ^ ~J^) 
in^nh e d e 

Terms multiplied by m^ — .00 560 

.... +.00131 

. . . . +.00085 

. . . . +.00044 



wr 



mr 



Higher terms by induction 

Sum — .00 300 

From the magnitude of the term concluded from induction we may conclude that 
this quantity is uncertain by — of its entire amount, in consequence of the slowness 

with which the series converges. . 



Formation of 

Terms independent of m . 

multiplied by )ry . 

m^ , 

m'. 

Higher terms by induction 

Sum 

Formation of 



Terms independent of m 
multiplied by m^ 



an I dk„ 
Mo wto e^ d X? 



— 


.2502 


+ 


.0590 


+ 


.0320 


+ 


.0111 


+ 


.0030 


— 


•I45I 



an I dK 



ni2 tn^ e de 



m" 



w* 



Higher terms by induction 
Sum 



A. p. — VOL. V, PT. III- 



— 1. 001 6 
+ .0335 
+ .0129 
+ .0034 
+ .0008 

- 9510 
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Formation of -?-^ JL ^Ji 
rn^m^e^ y ay 

Terms independent of m + -003 

multiplied by m^ — .119 



Sum —.116 

T,i ^ . r^ a n I d k, 

ronnahon of 5 ~ — ' 

m^ 7n2 e^ dv 

Tenns multiplied hy m^ — ,0506 

m^ —.0290 

w* —.0102 

Higher terms by induction — 0034 



'Sum — .0932 



Formation of g , 



.aw I dk^ 
m^ nh e^ e' d e' 



Terms in m^ —0.027 

Tj, .. n a n I dk^ 

tormation of ^ ^=— ? 

m2 m^ y dv 

Terms independent of m — .99 852 

multiplied by m^ — .01 128 

m^ + .00 245 

w* + .00035 

Sum —1.00700 

Formation of = -^ 

^2 ^^3 ey^ de 

Terms independent of m +1.98 

multiplied by m^ —0.12 

Sum +1.86 

Formation of -=—^ 

m^m^ y d y 

Terms independent of m — 3-9942 

multiplied hy m^ — .0064 

m^ + .0010 

m* + .0001 

Sum — 3.9995 
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Formation of « -^ 

in2fn^ y djo 

Terms multiplied by w^ — .00 968 

m^ +-00 220 

w* + .00 033 



195 



Sum —.00715 

We therefore have the following numerical values of the derivatives, in which, 
however, the last figure is nearly always to be regarded as quite uncertain, 

a JO an 

dk, o nunu 

T— 'zz— .0138^ ^ ^ 

de an 

dy an 

dk. . ^ nu Wa 

^;= + -oo4 749-^* 
av a n 

dk, twum^ 

-^= —.00 30oe -? — ' 
de an 



d k, ^ ^2 ^8 

dv an 

d k, nunh / -\ 

-^=-.95106-^ (95) 

a e an 

-^= — .000 35 y -^^-^ 
dy an 

^T=-.ooo28i^^ 
a» an 

dk, o / ^ ^3 

-.;-^=— .000 08 e' — ^ — - 
de an 

^ = -.002 028???»^ 

dv a n 

dkn , num^ 

— ±=+.003 756 -f-^ 
a e an 

dk^ ^ fn^m^ 

^ =z- 3.9995 r-J-^ 
ay a n 
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^ * zz — .000 014 ~-^- 
ax? a n 

(95) 
^=-.000 032/ ♦^^i^ 

a e an 

The values of — -, etc., may be obtained either by solving the equations (89), when 

they will be the coefficients of -^, etc., or by solving the equations (91). We thus 
find 






2.0225 



WgWts 



/'dw\ an 

(d V \_^ an 

\dk0 ^~ ' ^ m^m^ 

/'d e \ ^ a n 

(- ) = — .000928 

\aky em^m^ 

{'d e \ a n 
(— -)=— 1.0514 



r^^) = + .oooo9 



a n 



a n 



\dky 



00 ! 045 



ynhfThi 



/d y\ a n 

{ :rr )=— -ooo 97 

^dky ynUiWhi 

/"d y\ a n 

( - /- ) = — 0.25004 

^ dk^ ^ ym^vh 

Now, resuming the equations (90), let us begin with those terms -=---, etc., which 

a t 

are dependent only on the variations of the solar elements, and which are given in 

(93). We have put 

dk, dk^dv' dk^d e' 

df^dn'Tt'^'de'dt 

d K_d k^dv' dk^de' 
dJ~dTo''dT'^d~e' dt 

d ke__d k^ dv .d kgd e' 
di ~dTo'~dr'^d7Tt 
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d k 
If we here substitute the numerical values of 3—,, etc., just found, and then sub- 

av 

d » 
stitute in (93) both these values and the values of jj^ etc., in (96), we have 

dv adv. ^ ,de* 

d e dv' ,de* 

dt at dt 

dy dv ,de' 

y^zz.000001 ^- .000 01 1 e -^-p 

dt dt dt 

and, by integration, 

(5»zz — .0096 (5 »' +.003 04 6,e'^ 

6.^ = — .000 5& Sx?' — .000 09 S.€'^ (96) 

S,y^z=. . — .000011 d.e'^ 

These quantities are to be substituted in n tt^, and 5„ to obtain the correspond- 
ing variations of e, ;r, and 0. We have first 

6 nzz—^nSjo 
2 

zr + ^ w|.oo96 S »'— .003 04 S.e'^ 

Taking the partial derivatives of tti and 0^ (p. 190) with respect to », »', and e' 
we find the following terms ot 

I dTTi 

n dn 

Terms in m^ 00619 

^i' 873 

^* 443 

m^ : . . 179 

m^ 68 ^ 

rn' 25 

m^ 10 

m^ 4 

Higher terms by induction 2 

Sum 02 223 



198 
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Formation of - -^ — -} 
n ax? 



Terms in Wj — -0123 

nh —.0131 

W4 —.0059 

mg —.0022 

me —.0008 

nh , —.0003 

Mg and upwards — .0002 



Sura 



- .0348 



+ 


.0063 


+ 


.0108 


+ 


.0075 


+ 


.0040 


+ 


.0032 


+ 


•0318 



Formation of — . -} 
n e de 

Terms in /»2 

W3 

^4 

m^ 

Higher terms by induction 

Sum 

The derivatives of d^ are much more convergent. We find 

- -^ =z— .005 61 
n at) 

I 2 dOi 
n e de 

The terms of f, ^, and 9, which are independent of the derivatives of R with 
respect to the lunar elements, are, by (92) 



tf fzz I Sndt 
d ^=1 STT^dt 
SOzz^Sd^dt 



(97) 
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We have, from what precedes 



2 



d^7r^'=.n\'\- ,022 23 d» — .0348 5»' + .03i8d.c'*> 

<S0i=n|--.oo5 61 (J»+.oii 64<S»'— .005 37<S.6'*M 

Substituting for 6 a its value in 8 »' and 6,e'^ just found, these values become 

8nz:zn\ .0144 <S»' — . 004 56<S.e'^| 

<S;ri=w{-.0350(S»' + .03i9(5.e'2} (98) 

de^=n\+,o\ i69d»' — .005 39<Sy2| 

The forms of the integrals deduced from these values will be different^according 
as we consider the secular or periodic terms. There are, however, no secular terms in 
6 to' . Corresponding to the adopted values of the planetary masses are the following 
secular perturbations of e' 

5.6^=1-8^73 T-o" .028 T^ 

T being the number of centuries from 1 800. From this and the adopted constant 
term of e' we derive 

(y.6'«=-o".2926 T-o".cx>o 56 T« 
and hence 

<j£= 5".6oT^+o''.oo7T^ 

<J;r = -39-.3 T^-o".05oT^ (99) 

8 0- 6^.64 T^+o" .008 T' 

From the manner in which these quantities are obtained it may be shown that the 
coefficient of T^ in 5 f depends mainly on the value of -=— ^^ and the value of 5 ;r mainly 

on that of -=— /. We may, therefore, conclude from the defect of convergence in the 

series which express these quantities that the first quantity, the secular acceleration of 

the mean longitude, is uncertain by more than — , and the second by - of its entire 

amount, while the third is nearly exact. Although all the observations hitherto made 
on the Moon would not enable us to detect errors of this magnitude, a more accurate 
determination would be desirable from a theoretical point of view. But this can not 
be made without further investigations into the problem of three bodies, which it is 
not the object of the present paper to undertake. 
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Let US now take up the periodic terms of S v' and S e' , If we put 

5»'zzae cos N + a, sin N 
6 e' -zzSc cos N+6, sin N 

where N is an angle increasing uniformly with the time and therefore of the form 
No+fe<, and if we substitute these values in the equations for dn^ d tt^^ and S 6^^ and 
integrate, we find 

6 ez=,+ .0144^1 a, sin N— a, cos N| 
— .000 155 y\ec sin N — e, cos NJ 
or, 

5£z=(.oi44ae— 000 155 O ^ sii^ N — (.0144a,— .000 155^,) y cos N (100) 



where we put 



n 

"=1 



It will be seen that v represents the quotient of the period of the inequality divided 
by the period of revolution of the Moon. 

Applying the same process to tt and 0, we have 

(5;rz=(— .035oa^+-Oo 106 6^)^ sin N 

+ (0350 a,— .00 106 e,) y cos N 

(lOl) 

S 0zz(.oi i7ac— .oooiS^c) KsinN + (— 01 17 a,+oooi8e,) kcosN 

The derivation of the numerical results of tliese formulae requires the determina- 
tion of the periodic inequalities of the solar elements produced by the action of each 
planet The results will be found in Chapter III. 

§ 13- 

REMARKS ON THE RELATION OF THE PRECEDING THEORY TO THAT OF 

DELAUNAY. 

The relation of the preceding theory to that of D el a una y is so remarkable as to 
deserve special attention. Delaunay continually reduced the differential equations to 
the canonical form 



rfL dR 


dL_ dR 


dt ~dl 


dt~ dh 


rfG_rfR 


dg dR 


di ~'dg 


dt~ dCr 


dll dR 


dh dR 


dt ~dh 


'dt~ rfH 
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50I 



where /, g^ and h correspond as follows to the variable angles in the theory in the pres- 
ent paper: 

g = 7r-e 

h—e 

Here R may be regarded as the potential of any disturbing force whatever tend- 
ing to change the lunar elements. Comparing Delaunay's R with the equations (14) 
of the present paper, and the preceding equations with (65), it will be seen that those 
parts of our A;,, A:,, and )c^ which are derived from the coordinates of the Moon alone, 
which parts are the only ones yet considered, should be related to the constants L, G, H, 
of Delaunay in the following way: 

A:,=/.,(H-G) 

The final values of these constants, in terms of a, e, and y, are given by Delaunay in 
Vol. II of his work, pp. 235-236. We could, therefore, immediately compare these 
expressions with the values of our A;,, etc., finally obtained, but for the circumstance 
that the elements a, e, and y here used are not identical with those in which the coord- 
inates of the Moon are finally expressed at the end of the work. To reduce the one 
to the other, certain substitutions must be made, which are given on page 800 The 
identity of the expressions which would be derived by making these substitutions in 
L, G, H, with the values which I have derived from the Moon's coordinates would 
afford so remarkable confirmation of the numerical correctness of Delaunay's results, 
that I have deemed it worth while to make the transformation in question. The results 
are as follows: 



LzzVa/i" 



i+(-.' + ^r'-x*-^^''-|e'*+^y'^"-ie'+^-x'^-~7^" 



3 2 

32 
225,2 9„2 825 .2 



8 ' 4 



16 



.23255,2^4 

512 



^^ 32 8^ 32 "^^^288 128^^192 



, / , 1261 , II '549 /2\ 5 , 3028 6 I 76741 7 
+ ( + — ,- + >— ^ ) w + Q - ^w + ^^ ^^ 



96 64 
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864 
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(jc=.s/a fl 



^-\^+h*^-l^-hy'^-v. 



i6 



16^ 



+ (_i+3^._^._ie'^_|e'^+9,.,.+^,._.^.,V 



3 2 



384 



^81 /« 2 2081 4\ i 

128 384 ^ 



+ 



c+^t'- «--i'"-^?-'^'^)»*+(+-l-Tf ■ -^ 



^56 8' 51 

-3995 g.2^327 263 ^.x. 
192 ^ 36 864 ^ 



+(+iii+i^^')»''+^-'+4e^-' 



864 



H=Vo)u 



2 5 o 32 ^^ 

^V ,^06^^32^ 2 8 ^384 192 "^ 



38i_;,2^_2o8£^, 27^,, ,x 2 
128 384 32 ^ 



_l.r 4.2^45 gs . 3_ y2_25 495 e/2g2\ m»+r+895._?8_ii3 2 
+ ^+Ts"6 ^64^ 512 ^ ^^^288 9216^ 



256 64 

.3995e'='4.327 263^.x 4 
^ 192 36864 ^ 



+ r+n6i+liJl?,/2>),„5^3028^, 76|4I^7 
^VT^ q6 ^ 6a -^ ^81 864 



96 64 



The above relation affords a remarkable method of deteftnining the secular varia- 
tions of the Moon's mean motion, and of the motions of its perigee and node. From 
the equations (94) of the preceding section it appears that these variations are imme- 
diately obtainable from the three equations 

d Ijd a xdj^d^ .d/Ldy_.d/Ld^_^ 
Tadl de dt^dydt de'dt 



d^d^,d^(^de_j^d^dY,d^d^_^ 
da dt '^de dt '^dy dt de' dt 



d^d^.d/Rd^,d/Rdry,c^d^_Q 
Tadtde dt dy dt de' dt 
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The value of the secular variation of e' is to be inserted in the last terms, and 
the equations are to be solved for 



^« ^^ and^'' 
dt'dt' It 



We then have 



(p€ _dnda 
d f d adt 

d^TT _d Ttid a.d rr^d e .d yr^d y .d TTyde* 
T¥~d~adt'^Tedt'^d~yit^Te' dt 

rf^ G_d e^d ad O^d e d O^dy d 9^d e' 
df~da dt "^'de dt'^dy di'^de' dt 

Tti and 01 being the secular motions of tt and 0, as given in the preceding section 
p. 190. 

The analytic expression of these variations in this way offers no difficulty as far 
as the terms of the seventh order. Above this, however, we meet with the difficulty 
that the transformations given by Delaunay, on page 800, do not extend in any case 
beyond terms of this order in e, y^ and m. Moreover, Delaunay's supplementary 
expression for R, on page 742, does not seem to afford the means of carrying the 
transformation further. Were it required to write the value of L, or of A;,, with a 
higher degree of precision than this, it might be possible to do so inductively from 
Delaunay's expression for the secular variation of the Moon's mean motion given in 
the Coniptes Bendm, Vol. 72, p. 496, which extends to terms in m*. This course does 
not, however, seem necessary, as the coefficient in question, when earned to the fifth 
power of w, is sufficiently accurate for all purposes except the computation of the secu- 
lar acceleration itself 

It might be remarked, in this connection, that Delaunay never published any 
details of the process by which he effected the transformation of the values of e, y, 
and a, beyond terms of the seventh order in w. It is understood that these were 
reserved for his intended third volume, which was never prepared for the press. 

§ 14. 

TRANSFORMATION OF THE PERTURBATIVE FUNCTION R AND OF ITS 

DERIVATIVES. 

We have in the first chapter reduced the function R to the following form: 
where 
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In this form, however, the function R is subject to a serious inconvenience, in 
containing the square of A, which is a function of the elements of the Sun, Moon, 
and planet. This makes it very difficult to select all the terms corresponding to any 
particular combination of the varying angles. We shall, therefore, express R in the 
form of a series of products, in which that portion of its derivatives which contains the 
disturbed elements of the Sun and Moon shall be separated from that portion which 
depends upon the elements of the planet. 

We have from the preceding expression for A 

A''=(x'+:r,)=':r^+(.v'+t/4)V+(^+^4)V 

+ 2(y'+yi) (^'+Zi)y^ 

+ 2(^+^4) (x^+Xi)zx 

We have now by difFerentiating the preceding expression for R with respect to 
any element a 



dR 
d 



R [CXt xf Xt\d3f 



p' p V d a 

Mfft — ^—^*^^\ 

"^Vr^' p' pVrfoj 
+...,{(3_(^-0-_.).^_| (.0.) 

^\ p» p^J da 

Z{x' +x^) {x / 4-^4) (Ay dx\ 

^ p* ' '^^da^ daJ 

^ p^ \ (la day] 

Putting the latter part of this expression in this form affords us a great advantage 
in the selection of the terms which may become sensible. We see that in each of the 
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last six products the first factor contains only the coordinates of the Sun and planet, 
and not those of the Moon. Therefore, when we omit, as we must in a first approxi- 
mation, the squares and products of the masses of the planets, the factors in question 
admit of being developed as a function of the mean anomalies of the Sun and planet 
£p.ch term of this development is of the form 

cos ^ ^ " •/ i'*/ cQg 

y being the mean anomaly of the Sun, g^ that of the planet, and K a function of the 
eccentricities, mean distances, and inclinations. As it is necessary to carry the devel- 
opment to high multiples of i and jy the terms, in the case of a single planet, may 
amount to several hundred in number. 

The second factor of each product is a function only of the coordinates of the 
Moon relatively to the Earth, which we have represented by x, y, ^r, or, rather, it is a 
function of the elements of the motion of the Moon around the Earth. But these 
necessarily include the elements of the Sun, multiplied by its disturbing force. As 
we have developed x, y, and js, each term of the second factor will be of the form 

sin 

In the first approximation we may regard tt' and ©', or the solar node and perigee 
as constant, and the terms in question may be put in the form 

cos^ 1 I -TJ I J I 4^ y ^j^g 

Although we have here four variable angles, all of whose multiples are to be 
included, yet they converge so rapidly that the terms we have retained in the lunar 
coordinates include all that can be expected to give any appreciable result. 

Having thus developed each factor, if we form their products, we shall find each 

term of ^ -, which contains the lunar elements, to be of the form 
a a 



2 cos V -^ / 



The derivatives of each of the elements contain these terms multiplied by a factor 
which is constant in the approximation. If, then, we represent by s the coefficient of 
the time in Niw, we see that the integrated values of all the elements of the Moon 
contain terms of the form 

^cos^ ^ 

and «, ;r, and <9 contain in addition terms of the form 

s^ cos ^ -^ ^ 
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We hence derive the two following theorems of the motion of the Moon disturbed 
by a planet: 

1. The elements of the orbit of the Moon disturbed by a planet contain an infinity of 
terms of the form 

2. If any stich term is of very long period it may become sensible owing to the smallness 
of the divisor, the first and second powers of which enter into the coefficient c. 

§15. 
FORM OF THE PERTURBATIONS ARISING FROM THE SECOND TERM OF R 

Let us take up the second part of the expression (102) which we shall call -j—^y 

and examine more fully the perturbations to which it may give rise. We have just 
shown that the coefficients which depend on the coordinates of the Sun and planetmay 
be developed in multiples of their mean anomalies. Let us then put 



^Cx'+XaY II jy- TT I XT/ • TT 

^^^ — -k— -— - -^zzKiCos U + K 1 sinU 
2 fj^ 2 p" 

^lyL+y^f - i -\ zz: K2 COS U + K'2 sin U 
2 p^ 2 p^ 

3i^ii)--i i3=K3 COS U + K'3 sin U (103) 

iS:l+^diy'+J_d:=,K, COS U+K\ sin U 
iSl+A^.(l±^.=K, ,os U + K', sin U 



P' 



3ifl±fi),_(^±^4) _ j^^ ^.^,^ u^ k;/^ gi^ u 



p' 

Here K and K' are numerical coefficients, calculated in the manner to be hereafter 
described, and U an angle of the fonn 

Let us also put 

X^ XT 

-2 = ^1 COS IM 
a 



-2 = ^2 cos N 
a 



(104) 
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-2 

-5=^8 COS N 



(104) 



a 

ZX . XT 

- 2 = ^fl 8111 N 

a 
Each N is of the fonn 

and each K a function of a, e, y, a', e' and y'. 

We therefore have for the derivatives relatively to «, 

d X .0 'XT 

2y J— = — ia^^asiuN 

2-2r-— = — ia^XasinN i}^S) 

a e 

dy , dx * 2 XT 

x-^ +y-=—=:ta^XA cos N 

dz , dy .0 XT 

y^= — h ^ -,-^z=i 0^X5 cos N 
a* d e 

dx . d z -2 XT 

Z-J- + x^r- z^ta^H^ cos N 
a * a f 

while the derivatives with respect to tt and are of the same form, except that i' and 
i'\ respectively, are to be substituted for i. For the derivatives relatively to », e, and 
y, we have 

dx 2 /- . ^ ^1 \ XT 



^-^ 2/ . ^'^Sx XT 

2 2r;T-Z=a^( 2X3+- ) cosN 
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etc., etc., etc. 

6* wC 2 ^ 'f 1 XT 



el V . dx 2 ^ ^4 • XT 
ae rfe de 



etc., etc., 



while the derivatives with respect to y are of the same form. * 

Substituting in the general form (102) the iBxpressions in (103), (105) and (106) 
we find for the result of any combination of terms 

^=ii^n,/i,(i ^C,H,^ T,K,^ T,K,+ ^,K,+ r<:,K,+ ^,h,) a' cos (NdbU) 



If we put for brevity 
we shall have 



(107) 



^=:ii'm,/i2a»j;jrco8 (N±U)-A:'8in(NdbU)[ (io8) 

|5?= U" ni, M,a' j A'cos (N ± U) - K' 8in (N ± U) } 
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These values of the derivatives of R being substituted in the equations (89) and 
(90), will give the values of the derivatives of the elements with respect to the time. 

§ 16. 

PINAL FORM OP THE DIPPBRBNTIAL EQUATIONS FOR THE VARIATIONS OF 

THE LUNAR ELEMENTS. 

Returning once more to the equations (89) and (90), we have to find what terms 
are produced by the preceding value of Rj. If, in these equations, we omit the deriv- 
atives with respect to ^ in the second members, the effect of which we have just com- 
puted, substitute for the derivatives of R those with respect to Rg just found, and for the 

factors ( Yir )' ^*^' their numerical values found in (96) the equations (89) become 

^—ti^c^n\ i.oii2i-.oi46i' + .oo37i"} {K co8(NdbU) — K' 8in(NdbU)} 

e^z=fiia^n\— .00051— .52571' +.0000 i"\ f K co8(NdbU)— K' 8iii(N±U)| (109) 

y^,^=Mt<^n {-.0005 i-.ooo5i' -.1250 i"\ {K co8(N±U)-K' 8m(N±U)} 
where we put, for brevity, 

or the ratio of the mass of the planet to the sum of the masses of the Earth and Moon. 
Integrating these equations, and putting, as before, Vj the quotient of the Moon's mean 
motion by the coefficient of the time in the angle N ± U, we have 

<J»i=/i4a'r{ i.oii2i— .oi46i'+.cx)37i"} {Ksin (N±U)+K' cos(N±U)} 
eSe=^^a^y\— .0005 i—.52$7i' + .0000 i''\ {Kain (Nzfc:U) + K' cos (N±n)i (no) 
ydyz=:M^a^y\— ,00051— ,00051'— .i25oi''\ {Ksin (NdbU) + K' cos (NdbU){ 

There can be no terms corresponding to the case of N±U=:o, because U con- 
tains neither £, ;r, nor <9, and therefore in this case N can contain neither £, ;r, nor 0, 
and therefore t, i\ and i" will all be zero, so that the coefficient of the angle will 

vanish. 

A, p.— VOL. V, PT. ni 8 
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These values of S Vy S e, and S y are now to be substituted in the derivatives of 
n, ;ri, and ©j, and the latter are then to be integrated. If we put 

jt d n J. 
on=-j — Ox? 
axf 

a 70 de dy 

ax? de ay 

the values of tf €, <J ;r, and S 6 will be given by the equations 

S €-=z \ Sn dt — I -^-^ dt 
J } dk, 

* '=!*'■•«- 1 si"' 

se—{s o^dt—[~dt 

Substituting the numerical values of ^— , — , etc., we have 
S fizz, — -nSx? 

2 

^ TT'i^zz.nl .02223 ^ JO — .0195 e S e — .0994 y S y\ 
S <9i=w{— .00561 <5;i7— .0236c <Je+.oi49x<Sx} 
Substituing the preceding vahies of S x?^ e S e, and y ^ y^ and integrating, we have 
S €=/4,a^ v' {—i.Si? i+>o22 i' —.oo5i''\ { — Kcos (N± U) + K' sin (NdbU)| 

STTzz^iU^y'-l .0225 i+.oo99i' + .oi24i"} { — Kco8(N±U) (m) 

de=fi^ct'y^\— .oo57i+.oi24i'— .ooigi"} {— Kco8(N±a) + K:' 8m(NdhU)| 
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We are now to find the values of -j^ , etc., from the equations 

TK ~^d\J dla '^\d kj d e ^ VJIj/ d y 
etc. etc. etc. etc. 

Substituting the numerical values of ^-j- and the values (io8) of -=--^, etc., we 

d k^ dJD 

have, for any one tenn of R^* 

rfR* a( 1^1 ^K dK rfKl. /XT , TT\ 

■ -=.fii^a^n\ 2.023 K +1.01 1 -^ .0005 - .0005 > sm (NdbU) 



— .0147 0.526 - — ,0005— -— [ cos(N±U) 

a» (3ac ' ydy] ^ ^ 



dk-"^' \ — ^" • -^ ' rf;r7 ^ e^rf^ ^ydy 

3 f IT' . ^K' dK! rfKM /XT . TIN 

+ fjL^arn\ 2.023 K +1.011 .0005 , .0005 — ^Vco8(NdhU) 

t djo e de ydy] ^ 

^=/..«'«j-.oa93K-.oi47^^ -0.526 ^^^^--.0005 -^-) «m (N±U) 

'^=,.a^n\+.oo7S K +.0037 f^ + . 00004^^^^-0.1.50^} sin (N±U) 

+A*4«'«j+.oo75 K' + .oo37 'j^ +.00004 -^-0.1250 "^1 cos (N±U) 
[ a 70 ea e ydy] 

Inte^ating, and substituting in the preceding expressions for «Se, Stt, and Sd, we 
find that the important terms of these expressions may be expressed in the form 

de—€, cos(N±U) + e, sin (NdbU) 

<Jar = ;r, cos (N±U)+a-, sin (N±U) (112) 

de = e, COS (N±U) + 0. sin (NiU) 
where we have 

E,=/x^aA( 1.517J— .022 <' + .oo5 j")i'='K+( 2.023K+i.oii|5)»'l 

e,=tJi, A ( — 1.^171+. 022 /'-.005 r')i'2K' + (-2.023K'-i.oii^— )»'l 

a',=:>U4a'l( — .0226 «■— .0100 »' — .0124 t"^ v^YL +(—.0293 K — .0147-j — 
I ' ^- dv 



—.526 



ede 
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a-,=/i4a»|( .0226t+.oioot' + .oi24i") »'*K'+ (+.0295 K' + .oi47^ 






e.- 



( 00571 — .oi24t' + .ooi9r)»'*K +( .0074K +.oo37t— (113) 



dK . 



dK' 



(9. =:/u^ a*j (—.0057 i + .oi 24 «' — .0019 i" V* K'-f ( —.0074 K' — .00374^ 

dK' 



+ .1250 



!)'[ 



These expressions do not apply to the case of N ± U=zo, since y will then become 

de 
infinite. In this case i, i, i!' all vanish, as already shown, and the values of ^-. etc., 

reduce to constants. Integration then gives 

de=fi^aU — 2,02^ Ko— I 'OH -J— y+. 0005— ^\nt 
I ax? ede ) 

S7r = M,a'\ .o293K'o+.oi47^'+o.526^^^|w^ (114) 

cJ0zzyU4a^j-.oo74K'o-.oo37^?+o.i25^4nf 
I av ycLy] 

We have next to find accurate or approximate values of the various quantities 
which enter in these expressions. Let us begin with the quantities x^, Xg, etc., which 
enter into the values of K. The quantities are by (104) the coefficients of cos N or sin 

N in the development of -g, ^^j ®*c? which powers and products are to be derived from 

the values of x, y, and z in (81) by the substitution of the analytical values of ki and k^ 
given in the same section. If we take the plane of the ecliptic at the mean epoch for 
that of xy^ the value of /' will be of the order of the disturbing forces, and may be 
neglected entirely when attention is confined to quantities of the first order relatively 
to the disturbing forces. Putting x'mo, the values of rr, y, and 2 become 

^=Ai cos { (I +i+j) l-iTT-iJ+f) e+{i'+j') I'-i' n' I 
+ A^C«)8 j (I -i-j) / + i;r+0+/) 0-(i'+/) /'+«' ^ } 

l=k, sill j (i +i+i) l-in-{jJr3') «+(*'+/) /'-*' ^ \ 
+ A:, sin j (I -x-j-) /+i;r+(y+/) (9-(i' +/) /' +i' n" | 

i= c sill { (ii+^O l-i^ 7r-(j,+j\) e+(i\+j\) l'-i\ TV' I 
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The position of the axis of x being entirely arbitrary we may suppose all the 
longitudes to be counted from the solar perigee. This will be effected by putting 
t'=o, which will give V •=.§' . Thus, x^ y, z, etc., will be expressed in terms of the 
four varying angles 

the three quantities /, ;r, and representing the mean longitude of the Moon, and the 
longitudes of its perigee and node, all counted from the solar perigee, and g' the mean 
anomaly of the Sun or of the Earth. We shall first give the development of x, y, and 
Zj retaining the solar perigee, tt' ^ so that the axis of X remains arbitrary. The first 
column gives the indices i, i', j, and/, or the coefficients of ^r, /, A, and A'. The next 
gives the coefficients of/, ;r, 0, /', and 7r\ derived from those indices. We then have 
the coefficients k^ and k^, which, being multiplied by the cosines of the angles, give the 

terms of -, and by the sines the terms of " 
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=2c»in N — (Continued.) 
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We next form the products x^^ 1/^ ^, x y, y z^ x z^ so as to obtain the expressions 
for K^^ H2, ^3, H^, ^5, ^6 defined in (104) and used in (107). The complete values to 
quantities of the third order, inclusive, are given in the following tables. We now 
replace ?, tt, 9, V by the quantities 

(7r)=7r-7r' 
{9)-e — 7r' 

omitting the multiplier of tt* entirely, whicli will be equivalent to supposing all the 
longitudes counted from the Sun's perihelion. 

The computation stops at terms of the third order, because this is about the limit 
of possibly sensible terms of long period, and our immediate object is rather to detect 
the few such terms that may possibly exist than to make an accurate computation of 
their values. If necessary a more accurate determination of them can then be made 
after they are found. 
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The values (107) of K and K' and of such of tlieir derivatives as enter into the 
equations (112) to (114) are next tx) be investigated. Owing to the great number of 
the separate factors Hi, K^, and K\, and the small immber of combinations which give 
rise to sensible mequalities, we shall adopt the course of deferring the reduction of the 
y^i to numbers until we find what values of those quantities are required. It is, however, 
practically inconvenient to compute individual values of K and K' ; we shall therefore 
show how all sensible values of these quantities niay be found by developing the 
expressions (103) in numerical series. 

Although the numerical method ado])ted in the next chapter seems to me that best 
adapted to the detection of ine(jualities of long period, it is quite possible that, after the 
argument of such an inequality is established, its coeflScient can be best found by the 
method of analytic development. 

A. p.— VOL.vV, PT. in 9 
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OHI^A^IPTEIK HI. 

NUMERICAL DEVELOPMENTS FOR DETERMINING THE ACTION OF THE 
SEPARATE PLANETS ON THE MOON. 



$ 17. 

DEVELOPMENT OF THOSE FACTORS IN THE PERTURBATIVE FUNCTION WHICH 
DEPEND UPON THE COORDINATES OF THE SUN AND OF THE PLANET. 

Having foulid the factors h^ in the expression (107) we have next to find Ki, Kg, 
etc., as defined in (103), and used in (107). The first and most difficult operation under 

this head is the development of 3 and 5, p being the distance of the planet from the 

center of gravity of the Earth and Moon. It is necessary to carry tliis development to 
much higher multiples of the mean anomalies, and to higher powers of the eccentricities, 
than are necessary in computing the perturbations of the planets. The labor of com- 
puting the numerical values of the terms of the development from their analytical 
expressions would be enormous. I have therefore adopted tlie numerical method, in 
which the development is first eflfected in terms of the eccentric anomalies. 

The original inventor of this method was, I believe, Cauchy, who explained it in 
a series of short papers published in the Comptes RenduSj between 1843 and 1845. A 
systematic exposition of it, according t^) the plan of Cauchy, is given by Puiseux in 
Annaies de V Observatoire Imperiale (h Paris, Tonie VIIj where, however, it is applied 
only to the computation of particular terms which may give rise to inequalities of long 
period and not to the general development. The same volume contains a paper by 
BouRGET, in which it is applied to the general development of the perturbative func- 
tion. The adaptation of the method to the numerical computation of the general 
development of the perturbative function and its derivatives is mainly due to Hansen, 
who has developed it very copiously in his Auseitiamlerse tzung einer zwevkmdssigen Meth- 
ode zur Berechnung der absoluten Sforunffen der kleinen Planeten : Erste Ahhandlung^ found 
in the Abhandlimgen der Koniglich Sddisischen GesellscJiaft der Wissenschaflen^ Band III^ 
Leipzig, 1857. Hansen's exposition of the method is so well adapted to numerical 
computation that I have generally followed it, and the papers of both Hansen and 
PuiSEUX are so copious that it is not necessary to present a detailed development of the 
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formulae employed. But, as I have found it couvenient to adopt forms of computation 
not expressly given by Hansen, which may possibly be useful in the general compu- 
tation of planetary perturbations, I shall present a brief derivation of the formulae as 
I have used them. As it is convenient to use a notation in some respects peculiar to the 
present chapter, we present that actually adopted. Let us put 

/, /*, the true anomalies of any two planets. 

u u', their eccentric anomalies. 

g ff\ their mean anomalies. 

w, w/ the distances of their perihelia from their common node. 

p, their linear distance. 

y, the mutual inclination of their orbits. 

(jp, (p\ their angles of eccentricity, so that ezzsiny 

We then have for the square of their distance 

p^zir^+r'^— 2 r / I cos (f+w) cos (f+w') + cos y sin (f+w) sin (/ -{-tv') \ 

If we substitute for r', / cos/", and / sin/' their values in terms of the eccentric 
anomaly 

r' -zio! (i — ^ cosw') 

r" cos/ rza' (cos w' — e') I 

/ sin/ rzfl' cos if' sin vi 

and retain r, /, and w, we find that the vahie of fr may be put in the form 

p^zz H — K cos -^ cos vd — K sin ^ sin n! +5cos V (ii 6) 

where H, K cos ^, K sin ^, and s have the following values: 

H=:r^+a'^+ 2a' e! r \ cos w' cos (/+ w;) + cos y sin iv' sin {f+w)\ 
K cos »/»== 2 a' r \ cos w' cos (/+^) +cos y sin w sin {f+w) \-\-2e' a'' 
K sin ^ zz 2 a' r cos if' \ — sin w' cos (/+ w) + cos y cos w' sin (/+ w) ] 

s=a'W^ 

The computation of these quantities may be simplified by introducing the auxiliary 
quantities A;,, A^j, and the angles Kj and Kg determined by the equations 
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ki sin K, zzcos y sin w' 
/:iC08Ki=co8t<;' 
A^a «iu Kgizsiuic?' 
A;2 cos Kgzzcos y cos m/ 



(117) 



We then have 



Hz=r2+a'2+2 a' e' A^ r cos (/+m;— Ki) 
K cos t/fz=, 2 (if ki r COS {f+iv— K^) + 2 e' a'^ 
K sin ^ z= 2 rt' r /tj cos 9' sin (/+ iv — Kg) 



(118) 



s = e'^a'^ 



These four quantities are functions of the elements of both planets, and of the coor- 
dinates of the unaccented planet, but do not contain the coordinates of the accented 
planet. Thus, by the equation (116) we have the theorem: 

If the position of a jylanet he giveUy the square of its distance from any otlier planet may 
he expressed as a rational and entire function of the sine and cosine of the eccentric anomaly 
of the latter. 

On this principle is based the (Jaitchy-Hansen method of developing the pertur- 
bative function. 

Both Cauchy and Hansen express tlie coefficients H, K cos ^, etc., as functions 
of the eccentric anomaly w, but for our present purposes I have deemed it more con- 
venient to retain the true anomaly/ Having computed them for any assumed value 
of /, we have next to transform the expression (116) into the form 

p'=\C^qi^m{q-u')\ { i~ry, cos (Q+ii')l (iiq) 

Developing this product and comparing with (i 16) we find that C, Q, 5, and qi are 
to be determined by the equations 

s 

C = H+6'sin2Q 

(120) 
g^cosQ^zK cos^ — C q^ cosQ 

g^sin Qi=K sin ^+t^ g^isinQ 

Let us substitute in the last members of these equations the first approximate 
values of C, g, and Q, obtained by putting gizno, namely, H, K, and 9. We shall 
then have, for a second approximation, which includes the first powers of s only 
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gco8QzzKcos^( 1—2) 
q sin Q = K sin ^( 1+^2^ 
gzrK(i — j^fco8 2^) 

sin Q z= sin ^ ( I + ^/ cos V ) 

The substitution of the last two expressions in the first two of (i 20) will give the 
following values of q and G, true to terms of the second order in 8 

qi = ^i^i + ^,cos2tl;J 

(122) 

C=zH+5(i+4^/cos^^)8in^T/^ 

It is most convenient to cotnpute C from this formulae before procaeding to com- 
pute q and Q. If we retain C in the last members of the last two equations (i 20) and 
substitute for qi its value just given, we shall have q cos Q and q sin Q to terms of the 
second order in s. First, we obtain from (121) 

CgiCOsQ=z^'?^Kcost/^( i-"f) 
q ^ K^y 

But from (121) and (122) we have to terms of tlie second order in s 

qi_ s . 2 H5 . 

which, being substituted in the Ui«t equations, gives 

C qi cos Q— ^ cos ^ ( I + j^2 (2 cos 2 if))— i ) 

=z j^ cos^( I + ^2 (i — 4«»i vO ) 

C .«? H 5 

CgisinQzz j^' sin^( 1 + ^2 (1 + 2 C0S2V')) 

= 1^ sin V' ( I — ^^2 (1—4 (•<>« V) ) 
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The substitution of these expressions in the last two equations (120) gives for the 
definitive values of q cos Q and g sin Q to terms of the second order in s 



(123) 



q cosQ zzKcos ^< I — ^^2 ( i + "K^""" ^2" sm V) 

q sinQ=:K8in^|T-f5j^2(i '^~Y}~+ K^ ^^^ ^^)( 

In these expressions the only quantities neglected are of the third order with respect 
to 5, and therefore of the sixth order with respect to the eccentricity, and no case is 
likely to occur in which they can become sensible. The logarithmic computation of 
C, g, Q, and q^ may be conveniently effected as follows: 

Compute the five quantities 5i, Fj, F2, F3, and f by means of the equations 



Si 








F,: 


= i+4S, 


COS 


V 


F.: 


= 1— 4«i 


sin 


V 


F,: 


nl+Si 






f= 


rsFisin 


V 




C = 


:H + C 







(124) 



W^ then have 

C F 

q sin Q IT ( I -!-5 j^2 ^ ) K sin ip 

(125) 

C 

gcosQi=(i — 5T^2f'2^\) Kcos ip 

s 

I put the factors into the form i ±/<, because they are then very easily computed 
by Zech's Table of additimi and subtraction logarithms. By entering the addition table 

with the argument ^^ we take out the logarithm of i + h, and by entering the subtrac- 
tion table we take out log — =- 

Having thus, for any assumed position of the one planet, expressed the squai*e of 
its distance from the other in the form (119), we have next to develop the negative and 
uneven powers of p. Here, again, I make use of the method of Hansen, which 
seems to be based upon the investigations of Gauss.* Let us put 

* The developments and rules on this and the next two pages are, to a great extent, superseded by those of Vol. 
Ill, pages 42-53 and 60-75. I have, however, retained them in the present form, as pertaining to the work, and as 
equally convenient with the more general ones of Vol. Ill, which apply to a ^vider range of cases. (See post p^ 237.) 
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pS=C-jcos(Q-m') (126) 

SO that we have 

P^=pI \ I -«/i cos ((?+«')} 

the last factor differing very little from uuity. 
We have then 

The first factor is to be developed into the form 

To-" =- '^^ «>';' cos i (Q-m') 

2i=- 00 

If the cofficients ?/^' be developed in powers of ^ their general expressions will be 

^ »,(«_*i<>±2)>+4) . • ^_^("+ 2 »•- 2) «/*-,.«;+! «i y. 

2'^"- 2.4.6. . ". . 2i C< '^4^' 2 4"*"2 *+^'C2'' 

Here we use Gauss' notation for the hypergeometric series 

We have also, by putting i-\-i for i in the above value of b^, and noticing that 

2" -" ' 2.4.6. . . . (2? + 2)C'+'* '^4 + "2~'4 + 2+'' +''C^) 

The ratio --^.y- of two consecutive coefficients may now be developed in the form 

of a continued fraction by the method of Gauss, as given in his Disquisitiones gener- 
ales circa seriem infinitam, etc,, Sectio Secunda. CaUing this ratio Pi^i, we have 

W+2i-2 q 



a 

I 



b 
I — — 



c 
I — — 



etc. 
where 
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„_ (n + 2 i)(2i + 2—tt) g* 

~ i6i(i+i) C* 

, _(»+2 i + 2) (2 i + 4— _?) ^ 
~ i6(i4-i)(e+2) C^ 

(n+2i+4)(2j + 6— «) f 

etc., 

each term being formed from the preceding by increasing i by unity. We thus have 
the form given by Hansen 



(127) 



where 



Pi- 








I — 


I— ^i + 3 

I —etc. 




A,=: 


_(w+2?- 
" 16 


— 2){2i—n) (^ 


F,= 


2 i4-«— 
2 i 


-2 q 





(128) 



Owing to the slow convergence of the continued fraction (127) it is necessarj" to 
compute its value for the highest value of i indei)endently. To derive the necessary 
equations we must develop b^^ in ]K)wers of a, a quantity to be determined by the 
equation 

the expression (126) being put in the form 

pl=k(i—2 a cos (Q—u')-{- ct^) 

which is La Place's form for development. 

I deem it unnecessar}" to go over the det-ails of the derivation, since I have 
obtained the result given b}- Hansen, whose formula? may be most conveniently 
employed as follows: 

Compute X from the equation 



sin z=z- 

c 



when we have 
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azztan-x (129) 

2 

For the highest value of i to which it is deemed necessary to carry the development, 
compute Yi by means of the continued fraction 

sec^ - X 
2 



— «2 



I— etc. 



where 



Put also 



^ 4K*+0 40+0(^+2) 

"*-4(i+2)(i+3) ''- 4(^+3)0+4)' 

„- (»+2 )(6-») « j,_(M:2'-|-4)(2J"+6-«) , 

"'-4(i+"4)(i+5) '" 4(^ + 5) ('-1-6) 

etc. etc. 

Q^ii_n±2i--2 
4« 

„«) -(2<"-") (» + 2 ^-2) _ 2 i -« ,i, 

'^'•-' 16/(7-1) ~4('-i) " 



(130) 



A geaeral table of the values of ©<„'' and //» for all values of i and n commonly 
used is easy to form, and is given for n-=.^ and «^5 in the following section. 

Then compute, for all values of i 

F,= ei:'8injt 

(131) 
Ajn/ilJ'sin^Z 

Commencing with the highest value of y^ compute in succession all the other 

values by the formulaj 

_ I 

^*~i-Viy,+, (132) 

when we have 

(133) 
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The logarithms of y^ in (132) can be computed with great facility by the table ot 
subtraction-logarithms. Then having all the values of log Pi we require the value of 
log 6J?\ This can be found most conveniently from Runkle's Tables for determining 
(he values of the Coefficients in the Pertiirhative Function of planetary motion which depend 
on the ratio of the mean distances, (Smithsonian Contributions to Knowledge, Washing- 
ton, 1855.) These tables give, for the argument ^, the coefficients fej*^ of cos i y in the 
development of 

(i — 2acos9) + «0~'=i'^fci*^ cosiy. 

The development we actually seek is that of 

_ n 

(C — q cos if) ^iz i ^6{j^ cos i if. 

To assimilate it to Runkle's development we must determine a and k so as to 
satisfy the equation 

C — gcosyrz:A'(l +0^—2 a cosy) 



Introducing the auxiliary \;ti 



sinzm^, 



(134) 



we have 

azztan^A^ 

k=CQm^X= ^ , 
2 tan ^ X 

Then, comparing like terms in the two developments, 

— ^ 

provided that hf^ be taken from Runkle's tables with the argument a. In the cases 
oi sz=.-- and 6"zz^*the quantities actually tabulated are, 

for^zz^; (i-ayhf'= ^-f'^-ftf 
2 ^ cos* i X 

for 5ZZ 5 ; ( I - c^y 6i«) - --^f /- ¥^' . 
2 ^ cos^ ^ X 

If, then, we represent by j5^"^ and /?J®^ the quantities actually taken from Runkle's 
tables when s-zz^ and s-zz^ , we have 



3 

2' ,,^2 



fcfzzC sec^zcos i;t/?r 

(135) 



6i^>=C sec* ;t cos^ i Z /^i^ 



(0) 
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These quantities may also be computed with equal ease, though, perhaps, not 
always with equal accuracy, by the use of the arithmetico-geometrical mean. Rep- 
resenting this mean in the usual way by M we have, first. 



fer> = - 



VCM(cosiz,Vcosz) vC^szHO, r^^~) 

^- vcosz ^ 
Then ¥^^ and b^^^ may be obtained by the formulae 



(136) 






C-?j9j"~C(i-K8inz) 



(137) 



■C-qpi,''-C(i-p'^>smx) 



When p^i' sin x a"d i*"' sin x approach to unity, these formulae will give 6J** and 
fe^" with only very limited accuracy and others must be used. 

The result of the preceding operations is that, for any assumed position of the 
unaccented planet in its orbit, we may develop the value of /r>o~" in the form 

P,-"=i l6J.*>co8t(Q-«') (138) 

-00 

We now compute this development for a number of positions corresponding to 
equidistant values of the mean anomaly, g^ of the planet, and thus obtain a separate 
value of ¥n^ and of Q for each value of the mean anomaly. Then, any functions of these 
quantities may be developed in a periodic series, proceeding according to the multiples 
of the mean anomaly by the well-known process of mechanical development. The 
most convenient functions to develop thus are fej^^ cos i {Q—g) and b^^^ sin i (Q— gr). If 
we put 

(^39) 
A z=.u' —g 

we have 

Po-" = i^^6i:^cos(iQ'-iA) 

zz J ^ b^^' cos i Q' cos i k-\-i^^ ¥^^ sin i Q' sin i A. 

Let us, then, develop />}j' cos i Q' and b^^ sin i Q' in a periodic series of the form 

¥^^ cos i Q' zz (^ f + c S*^ cos g + c ^'* cos 2 ^ +, etc. 

+s J*^ sin ^9+5 2*^ sin 2/7+, etc. 

(140) 
b'^^ siniQ'zzc'i*^ + c'l'^ no^g + c'^i^ cos 2//+, etc. 

+ •^'1** m\g-\-s'^2^ sin 2^+, etc. 

Then, let us compute values of (a^^) and (fe^j) for all positive values of i and^' from 
the formula? 
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(ai..) = -/f (M = c"? (141) 

(«i. -)) = \ (O'"' + ^'f ) (ft.-,-.) = \ ic'r-^f) 

and we shall have the numerical expression for piT" for all positions of the two planets 
by the formula 

00 +00 ( I 

Po""=:-^^i ^i (aij) cos (iA+j^) + (fc,^.) sin (i A+je/) (142) 

-00 I I 

To obtain p""**, this development is still to be multiplied by the factor 

2 2 

Fzz { I — (7i cos (Q +w') I ~ n zz { I — ji cos Q' cos (A-f 2 (7) + g'i sin Q' sin (A+ 25^) } - « 

Owing to the minuteness of the factor q^ we may nearly always neglect its pow- 
ers, and so put 

Fzzi + - (/iCOsQ' cos(A + 2 5r) gisinQ' sin(A+2 5r) 

Let us develop as before - q^ cos Q' and q^ sin Q' into a periodic series of the form 

- qi cos Q' =.po +Pi cos g +P2 cos 2 17 + etc. 
2 

+ ji sin /7 + ^2 sin 2 // + etc. 

(■43) 

+^'1 sin /7 + 3^2 sin 2. 9 + etc. 
and we shall have 

F = I +po cos (A + 2 ^) —p'o sin (A + 2 ff) 
+2 i (p,+ q\) cos ( A+ (2^i)g)+2^(p,-q\) cos (A + (2-0.'7) 
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Multiplying the series (142) and (144) term by term, we shall have the value of 
p~" expressed in the form (142). The development will be given in terms of the mean 
anomaly of the one planet and the eccentric anomaly of the other. 

The next process is to change it into a development in terms of the mean anoma- 
lies of both planets. This is done by means of the Besselian transcendents, tables 
and formulae for which are given by Bessel in the Berlin Memoirs for 1824, while the 
mode of applying them is fully explained by Hansen in his ^^AuseinandersetzungP 

Ha\dng thus effected the development of -5 in terms of the mean anomalies of the 

Earth and planet, we have only to express x\ y\ /, the rectangular coordinates of 
the Sun relatively to the center of gravity of the Earth and Moon, and X4, y^, ^4, the 
rectangulai' coordinates of the planet relatively to the Sun, in terms of the respective 
mean anomalies. The factors depending on the coordinates of the planet and Sun will 
then be formed by multiplication. As the complete formation of these terms would 
be very laborious, we shall only actually develop those of which we find ourselves in 
need, and these can not be selected until we proceed with the more paiticular consid- 
eration of the action of the separate planets. 



§ 18. 

NUMERICAL DETAILS OF QUANTITIES DEPENDING ON THE POSITION OF THE 
PLANETS MEECUEY, VENUS, AND MARS. 

The computations of which the results are now given were made in the years 
1871-73, before I had prepared the more general discussion of the method found in 
Vol. Ill, Part I, Chai)ter8 II and IV. The latter does not, however, involve any mate- 
rial change in the method so far as the present application is concerned, where no 
deviation from Hansen's procedure is necessary. The object of the present section is 
to present such details as will enable the results to be verified, corrected, or reproduced 
with the least amount of labor. 

The computations were made with more or less fullness for Mercury, Venus, and 
Mars. In the case of Mercury, however, the development was in error from the use 
of the longitude of the solar perigee, when that of the Earth's perihelion should have 
been used. As this error would but slightly change the larger terms in the develop- 
ment, and therefore would not invalidate the general conclusicm that the action of 
Mercury produce no inequalities of importance ifi the Moon's motion, it was not cor- 
rected by a recomputation. 

In the case of Mars the developments were not completed beyond the point where 
the eccentric has to be replaced by the mean anomaly, because the order of magnitude 
of the coefficients in the development according to the eccentric anomaly would suffice 
to draw general conclusions as to the })ossibility of hidden terms of long period. 

In making the computations the first question is of which eccentric anomaly the 
development shall be a function. Theoretically this is indifferent; in practice, how- 
ever, it will be convenient to choose the eccentric anomaly of the orbit having the 
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smaller eccentricity, because the Besselian transformation to mean anomaly will then 
be simpler, as will the change of Po" i^to p"^** . 

It will then be better to assume special equidistant values of the mean anomaly 
of the other planet rather than of its eccentric anomaly, because the necessity of a 
transformation to mean anomaly of this planet will then be avoided. For each position 
of the planet thus arising the quantities K, ^, C, etc., are to be computed. No example 
of the way of effecting this computation is necessary. It is different, however, with 
the coefficients 6{J\,the computation of which may be greatly facilitated by a suitable 
arrangement of the work, as will be shown in a specimen of the computation. On 
page 247 is given a transcript from the computing book of all the figures made in 
computing a set of values of fejj^ up to i 11:20, except the computation of 6J?\ p!i\ and 
the last value of x^. The process is as follows: 

The values of log A^ are formed by adding to each yu^*^ the logarithm of sin ^x 

log sin^;^ = 9.9583004 

and are written in the second" column opposite the corresponding values of i. 

The last value of log Xi (in the computation log y^, but, in printing, log y^) is 
then written under the corresponding value of A<. The sum of the two logarithms 
is written opposite i in the third column and immediately under this sum its arithmet- 
ical complement, to serve as the argument in Zech's table of subtraction logarithms. 
With this argument the value of log yi_i is taken from the table and written under log 
A<_i. This process is repeated upward until the value of Xi is reached. 

To render easy the detection of an error the remainder of the work is done by 
the continual addition of diflferences. Since 

logpizzlog sin x+^og e^+log y^ 
we have 

A log 2h= A log e^+ A log yi 

The values of A log yi are formed by subtraction of the consecutive values of log 
yi and those of log 0^ by taking the differences of the Table of values of log Oi. The 
remainder of the operation of forming log ¥*^ needs no explanation. 

The above computation has to be made for each of a certain number of equidis- 
tant values of g, extending through the circle. In the case of Venus and Mars I have 
taken 

g=o^; 15°; 30^ 345° 

thus dividing the circle into 24 parts. 

The next step is, for each value of g and of i, to compute 

6<*> cos i (Q—g) and b<*^ sin i (Q—g) (145) 

where, in each case, the special numerical value of g, used in computing 6^*^ and Q, is 
to be used. 

With the 24 special values of each of the quantities (145) its general expression in 
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a trigonometric series, proceeding according to sines and cosines of g^ is to be formed. 
I deem Encke's method of forming tliis series to be the most convenient. For the 
case of n:=24 the formulae are as follows: 

We put 

(w)= numerical value of the function oi ffzui^^y^m hi (145) 
We compute 

C'z^Co); C*=(A) + (24-A) 



C>^=(i2); S*=(A:)-(24-A;) 



(Am; 2; 3; 



lO 






i2*(to+f„)=Ct,+CU 

6(c,+c„)zz«i + a2 
6(ci— fi,)=«3+«« 



C_=C»-C»*-* (A;=o; i; 2; 3; 4; S) 
Si=:S*-S'^-* (A=i; 2; 3; 4; 5) 

C*+_=C*.-Cr (A=o; i; 2) 

0*_-C*-0«r* (Am) 

S*+.zrS*-SV (Am) 
Si_=S'_-S«r* (A=i; 2) 

^=72G«--^0 



^*=W3St 



6 (»i + 6'„) = fc, — 62 

6(s,— s,i)zr63+ft4 
6 (.,-.,„) rz^^/ 3 SL, 



6(C3 + 6-,)=Cl-Cl 



6(63 + «,) =yj^h 
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6 (cs+^t) =«!— «2 6 (65+57) zz6i-f-62 

6 (c^—Ct) =03— a4 6 (55—57) zrfta— 64 

It should be noted that the indices of C and S are not powers. The preceding 
formulae. may be controlled at pleasure by the general conditions 

ri = (o) + (i)cosi.i5° + (2)cosi.30° + (3)cosi.45°+ . . . 
5,= (i)sini.i5°+(2)8ini.30° + (3)sini.45°+ . . . 

or by 

c'o-f c'lCOs^. I5°+C2C08 2 /*. i5° + r3Cos3 w. 15°+ . . . 
+5i sin n . 1 5°+52 sin 2 /* . 1 5^+53 sin 3^.15° + . =(^) 

The values of Ci and 5^ obtained in this way from the special values of fejj* cos i Q', 
and the. accented coefficients c' and 5' obtained from fe^*^ sin i Q' for each value of i, 
being substituted in (141), give the required development of Po^ in the form (142). 
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The following form for the principal part of the computations was used: - 








I 


2 


3 


4 


5 


6 




(0) 


0) 


(2) 


(3) 


(4) 


(5) 








(23) 


(") 


(21) 


(20) 


(19) 








(lO 


(10) 


(9) 


(8) 


(7) 


(6) 




(12) 


(13) 


(14) 


(15) 


(16) 


(17) 


(18) 


c 


c 


c 


Q2 


C 


C* 


C» 


C 




c« 


C" 


QIO 


C» 


C« 


C^ 




s 




S' 


S^ 


S» 


S* 


s» 


S« 


c. 


c<» 




gio 


S' 
C'+ 


S» 


s^ 












ct 


G% 


c*+ 

log Gl 










c_ 


G!_ 


CL 


CL 


OL 










iOi 


C*_ 


CL 

CL 


iCL. 








C++ 


^r 


CV+ 


cu 


ct 










iC++ 


ci 


c».+ 


icw 








c+. 




log cv_ 


C'+- 










c_+ 




Ci+ 
















log CU 












c__ 




C!._ 


s» 










s. 






SI 

04 

logV 


SI 








s_ 




SL 


SL 
SL 


SL 
iSL+ 








S++ 




SV+ 


log 8U 










S+- 




S'+ 












S-+ 




Si+ 


Si+ 










s__ 


. V. PT. m- 


SL_ 
10 


SL_ 










A. p.— VOL 
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In presenting the following numerical data and results the writer has been guided 
by the consideration that future investigators may desire to repeat his work with a 
higher degree of numerical precision than that which he has himself reached. In this 
connection it may be remarked that complete numerical consistency in the use of the 
adopted elements seems to be necessary. For example, serious error might result in 
the terms of very high order if one value of the eccentricity were used in one part of 
the computation and another value in another, and this, although neither of the two 
values should be seriously in error. It is also to be remembered that small changes 
in the elements, especially the positions of the perihelia and nodes, may make consid- 
erable changes in the coefficients of the terms in question. The epoch 1800 was 
therefore chosen as that to which the developments should refer, because this was fairly 
near to the mean of the time through which observations of the Moon's position were 
available. If the elements for 1900 were used, it is possible that slightly diflferent 
results might be found for the inequalities of very long period, as, for instance, Hansen's 
inequality due to the action of Venus. 

In the actual computations the true anomalies and radii vectores of the Earth and 
Mars for different values of their mean anomaly enter as fundamental numerical data 
in the computation. These anomalies and radii vectores were taken for 1 800 out of Le 
Verrier's tables. On the other hand, whee the eccentricities enter into the formulae, 
they were reduced to 1800 with Le Verrier's data. It is not impossible that in this 
way numerical inconsistencies, amounting to a large fraction of a second of arc, have 
crept into the computation, especially as account was not taken of the minute error in 
the secular variation of the Earth's radius vector which exists in Le Verrier's tables 
of the Sun. It must not, therefore, be understood that the whole computation is 
absolutely consistent with the following nuraerical elements which were intended to 
form their basis. I shall therefore present such derived quantities as will enable the 
general consistency of the work to be most easily tested. 

In the notation of § 1 7 it will be seen that the position of the planet whose elements 
are accented remains indeterminate, while the computations are made for special values 
of the mean anomaly which is unaccented. Now, in the respective cases of Venus- 
Earth and Earth-Mars, it is most convenient to let the outer planet play the part of the 
unaccented one, since its eccentricity is in each case the largest. 

The values of the fundamental elements actually used in the computations were 
as follows: 
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Venas and Earth. 


Earth and Mars. 




' 


II 


/ 


II 


ir 


99 29 


56 


332 22 


42 


n" 


• 1 28 46 


I 


99 29 


56 


e 


254 52 


28 


48 


33 


y 


3 23 


32.5 


I 51 


3-5 


n — dzzzco 


204 37 


28 


284 22 


9 


Tt' — d—co' 


233 53 


il 


51 29 


23 


a)-Ki 


330 41 


3 


232 53 


38.3 


a-Kj 


330 46 


48 


232 51 


537 


e' 


0.0068703 




0.016792 




log^i 


9.9997360 




9.999861 




log*2 


9.9995030 




9.999912 




log 2 e' a" 


7.856682 




8.52614 




logs 


5-3926 




6.4502 





The data and principal results of the computations made by the preceding method 
are given in the pages next following with some fullness, not only on account of their 
bearing upon the problem immediately in hand, but of their possible usefulness in the 
theories of Venus, the Earth, and Mars. In the Hansenian method of developing the 
perturbative function, the inverse first, third, and fifth powers of the distances of the two 
planets must first be developed as the basis of the whole computation. From this 
development the function itself and its derivatives is derived by the process called by 
Hansen "mechanical multipHcation;" that is to say, the factors by which the inverse 
powers are multiplied are expressed as functions of the mean anomaly of one or the 
other of the planets, developable in sines or cosines of multiples of that quantity. The 
coefficients in this development are pure numbers, so that only a numerical multiplica- 
tion is necessary. 

We first have to compute the values of the coefficients which enter into the expres- 
sions (100), (i 16), (i 19), and (126). This computation is made by the formulae (124) 
and (125). 

For the separate computations the circle has been divided into 24 parts; a number 
greater is rigorously necessary, but advisable in order that the general consistency of 
the results may afford a check upon the accuracy of the computations. 
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VENUS AND THE EARTH. 
Values of the principal quantities for each 15° of the EartKs mean anomaly. 



s 


\Q/g K sin 


log K cos ^ 


C 


Q 


\ogq 


log^tj) 


log^(J) 


log sin a: 


o 
o 


9. 842601 3n 


0.0959176 


1. 4984309 


' // 

330 50 5. 3 


0. 1547869 


1.0393590 


2. 0106083 


9.9791502 


'5 


9. 530781 1« 


0. 1424754 


1.5005595 


346 15 34.0 


0. 1550778 


1.0321281 


1.9961461 


9.9788246 


30 


8.6167445 


0. 1555235 


1.5042569 


1 39 24.0 


0. 1556971 


1.0219812 


1.9759630 


9.9783751 


45 


9. 6225 171 


0. 1372497 


1.5092620 


16 59 52. 6 


0. 1566407 


1.0105490 


1.9533528 


9.9778760 


60 


9.8851684 


0. 0850808 


I. 5152267 


32 15 22.8 


0. 1578726 


0.9992563 


I. 931 1747 


9. 9773950 


75 


0.0263142 


9. 9897552 


I. 5217372 


47 24 35- 2 


0. 1 593 1 88 


0.9890811 


I. 91 13730 


9. 9769792 


90 


0. 1085707 


9. 8261023 


1.5283468 


62 26 37. 7 


0. 1608716 


0. 9805013 


1.8948515 


9. 9766495 


105 


0. 1 5 17266 


9.5027179 


I. 5346068 


77 21 13. 2 


0. 1624003 


0. 9735970 


1.8817050 


9.9764032 


120 


0. 163455 1 


3. 7370093" 


I. 5400942 


92 8 42. 6 


0. 1637673 


0. 9681470 


1.87 141 28 


9. 9762200 


135 


0. 1458188 


9. 6266576^ 


1.5444443 


106 50 2.9 


0. 1648478 


0. 9638634 


1.8633261 


9. 9760756 


150 


0. 0965588 


9. 882956811 


t. 5473683 


121 26 41.9 


0. 1655454 


0. 9605464 


1.8569903 


9- 9759517 


165 


0. 0075048 


0. 022802011 


I. 5486731 


136 28.2 


0. 1658028 


0.9581767 


1.8523352 


9.9758430 


180 


-1-9.8S7'886 


0. io55524» 


I- 54g2730 


150 33 21.5 


0. 1656082 


0.9569734 


1.8497980 


9.9757607 


195 


9. 5745021 


0. i5oi937« 


I. 5461945 


165 7 21. I 


0. 1649943 


0. 9573378 


1.8502075 


9. 9757301 


210 


7. 8244550 


0. 1640344^ 


I- 5425743 


179 44 16.2 


0. 1640312 


0. 9597802 


1. 8546328 


9.9757851 


225 


9. 5592682» 


0. 1489068^ 


I. 5376514 


"94 25 38.0 


0. 1628151 


0. 9647424 


1.8640192 


9- 9759572 


240 


9. 849868211 


0. io23996« 


'•5317533 


209 12 33-3 


0. 1614556 


0. 9724792 


1.8789404 


9. 9762668 


255 


0. 0025697^ 


0. 01 6304411 


1.5252760 


224 5 42. 6 


0. 1600604 


0. 9828938 


I. 8992569 


9.9767119 


270 


0.092 1 85on 


9. 8694460^ 


1.5186557 


239 5 20.2 


0. 1587229 


0.9954332 


1.9238996 


9.9772636 


285 


0. 1407563^ 


9.5928653^ 


1.5123460 


254 u 16.3. 


0. 1575166 


1.0090239 


1.9507358 


9. 9778654 


300 


0. i564598» 


8. 1880642^ 


1. 5067814 


269 23 1.8 


0. 1564928 


1.0221600 


1.9767483 


9. 9784426 


315 


0. I4I2940W 


9. 5590788 


1.5023484 


284 39 51- 5 


0. 1556843 


1.0330654 


1.9983638 


9.9789136 


330 


0. 09258o6« 


9. 8542650 


1.4993584 


300 47. 2 


0. 1551152 


1. 0401244 


2.0123302 


9. 9792098 


345 


o.ooii484» 


0. 0073922 


1.4980207 


315 24 38.9 


0. 1548076 


1.0422776 


2. 01 65 1 88 


9.9792898 
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MARS AND THE EARTH. 

Vcdttes of the principal quantities for each 15° of the mean anomaly of Mars. 



9 


log K sin i> 


logKcos^ 


C 


\ogq 


Q 


log^io) 


log^ 


log sin ;t 






0.34285811 


0. 21297911 


2. 881 190 


0.438017 


/ // 
233 26 53. 9 


0. 599924 


6.0122 


9. 978445 


15 


0. 4i873^« 


9. 939I90W 


2.910317 


0.441407 


251 39 51.6 


0. 575359 


6.0088 


9. 977466 


30 


0. 448288^ 


8. 38426511 


2.971771 


0.448350 


269 30 20. 5 


0. 525201 


6.0018 


9-975335 


45 


o.43909i» 


9.917232 


3. 0601 16 


0. 457926 


286 43 58.9 


0. 457934 


5- 9923 


9.972189 


60 


o.39i572« 


0. 207560 


3. 167944 


0. 469034 


303 12 14.5 


0. 382842 


5.9812 


9.968257 


75 


0. 2986361* 


0. 357601 


3. 286834 


0. 480623 


318 52 19.0 


0. 307729 


5.9696 


9.963845 


90 


0. I37266« 


0. 444640 


3-408158 


0. 491815 


333 45 49. 5 


0. 238109 


5.958V 


9.959295 


105 


9. 82i323» 


0.492325 


3- 523778 


0. 501949 


347 57 24. 5 


0. 177420 


5.9483 


9. 954940 


120 


8.944271 


0.510437 


3.626471 


0.510555 


I 33 20.6 


0. 127466 


5- 9396 


9.951070 


135 


9. 921101 


0.502957 


3.710223 


0.517327 


14 40 44. 9 


0.089206 


5. 9329 


9. 947927 


150 


0. 185703 


0.470233 


3- 770379 


0. 522066 


27 27 0.6 


0. 062981 


5. 9281 


9.945680 


165 


0. 332616 


0.408997 


3. 803729 


0. 524655 


39 59 34. 1 


0.048896 


5.9255 


9.944445 


180 


0. 424058 


0.310220 


3.808538 


0. 525041 


52 25 48-2 


0. 047016 


6. 9252 


9.944282 


195 


0. 480036 


0. 15 1066 


3- 784559 


0. 523210 


64 53 5. 7 


0. 057224 


5. 9270 


9.945195 


210 


0. 508707 


9. 855220 


3- 733036 


0.519199 


77 28 52.3 


0. 079588 


5.9310 


9- 947137 


225 


0.513048 


8. 291540I9 


3.656695 


0.513098 


90 20 38. 3 


0. 113981 


5.9371 


9. 950010 


240 


0. 492674 


9.876442^ 


3. 5597 1 1 


0. 505068 


103 36 0.2 


0. 160083 


5.9451 


9. 953654 


255 


0. 443756 


0. 158005W 


3- 447645 


0. 495379 


117 22 30.6 


0. 217243 


5.9548 


9.957857 


270 


0.3569H 


0. 3082 1 3» 


3. 327367 


0.484444 


131 47 19.0 


0. 283875 


5.9658 


9.962343 


285 


0. 209610 


0. 3962 10« 


3. 206725 


.0.472863 


146 56 28.9 


0.357230 


5.9773 


9. 966801 


300 


9. 929867 


0. 441856^ 


3.094241 


0. 461452 


162 53 52.3 


0.432631 


5.9887 


9. 970898 


315 


8. 223455» 


0.451265^ 


2.998513 


0.451227 


179 39 39.0 


0.503153 


5.9990 


9. 974321 


330 


9.9I2794« 


o.4236oi« 


2. 927385 


0.443299 


197 8 46.0 


0.560062 


6.0069 


9.976819 


345 


0. I9903i;» 


o.35ii98« 


2. 887057 


0. 438679 


215 lo 3.2 


0. 594434 


6. 01 15 


9.978223 
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Table of the valves of log 0J'^ and log >uJJ^/or the computation ofU^K 



i 




«==3 






n = 5 






log ec) 


2^ log 6(0 
791813 


log //«) 


log e«) 


^ log 0(0 
1549020 


log //(O 


I 


9. 8750613 


0.0969100 




2 


9. 7958800 


299632 


9. 1938200 


9. 9420080 


669467 


9. 339948o« 


3 


9. 7659168 




9- 3399481 


9.8750613 


377886 


8.9719713 


4 


9.7501225 


157943 


9-3699113 


9. 8372727 


243593 


9.2352128 


5 


9. 7403627 


97598 


9. 3813407 


9. 81 29134 


170333 


9- 3077634 


6 


9. 7337322 


66305 


9. 3869447 


9. 7958801 


125892 


9- 3399481 


7 


9. 7289333 


47989 


9. 3901 148 


9. 7832909 


96873 


9- 3573222 


8 


9. 7252989 


36344 


9. 3920843 


9- 7736036 


76868 


9.3678383 


9 


9. 722451 1 


28478 


9- 3933924 


9. 7659168 


62490 


9.3747102 


10 


9.7201593 


22918 


9- 3943057 


9. 7596678 


51805 


9. 3794566 


XI 


9. 7 18275 1 


18842 


9. 3949687 


9. 7544873 


43647 


9. 3828762 


12 


9.7166988 


^5763 


9- 3^54654 


9.7501226 


37279 


9- 3854235 


13 


9- 7153605 


'3383 


9. 3958471 


9- 7463947 


32210 


9. 3873728 


14 


9. 7142100 


1 1505 


9. 3961467 


9-7431737 


281 1 1 


9.3888982 


15 


9.7132104 


9996 


9. 3963862 


9. 7403626 


24746 


9. 3901 146 


16 


9.7123339 


8765 


9. 3965806 


9. 7378880 


21952 


9. 391 1005 


17 


9.711559' 


7748 


9 3967408 


9. 7356928 


19607 


9. 3919108 


18 


9. 7108692 


6899 


9. 3968742 


9. 7337321 


17618 


9- 3925849 


11 


9. 7 1025 10 


6182 


9. 3969865 


9-7319703 


15918 


9.393i5'7 


20 


9. 7096939 


5571 


9. 3970820 


9- 7303785 


14453 


9- 3936329 


21 


9.7091892 


5047 


9-3971638 


9- 7289332 




9. 3940449 


22 


9. 7087298 


4594 


9. 3972344 


9. 7276152 


13180 


9. 3944005 


23 


9. 7083101 


4197 


9. 3972959 


9. 7264083 


12069 


9- 3947095 


24 


9. 7079249 


3852 


9. 3973496 


9. 7252990 


1 1093 


9. 3949797 


25 


9. 7075702 


3547 


9. 3973Q69 


9- 7242759 


1023 1 


9. 3952172 


26 


9. 7072426 


3276 


9. 3974387 


9. 7233294 


9465 


9. 3954273 


27 


9. 7069389 


3037 


9. 3974758 


9. 72245 n 


8783 


9. 3956139 


28 


9. 7066569 


2820 


9. 3975090 


9. 7216340 


8171 


9. 3957804 


29 


9. 7063940 


2629 


9- 3975387 


9.7208718 


7622 


9. 3959297 


30 


9. 7061486 


2454 


9. 3975655 


9- 7201593 


7125 


9. 3960640 


31 


9. 7059188 


2298 


9. 3975896 


9.7194917 


6676 


9. 3961854 


32 


9- 7057034 


2154 


9. 39761 15 


9. 7188648 


6269 


9. 3962948 



. The following is a complete copy of the computation of 6^*^ for Venus when 
^=0°, showing every figure that was made, as explained on pages 233-234: 
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Venus : Computation of If for g {Earth) =0^. 





log A, 

logr, 


logA^r, 


A log Xi 
A log/^ 


I 


9.9126633 




-1369007 


2 


-9. 2982484 


9.3478124 


- 180013 




0. 0495640 


0.6521876 


529996 


3 


+8.9302717 


9- 0328353 


— 139471 




0. 1025636 


0.9671647 


268591 


4 


9- 1935 132 


9. 3229359 


— '09295 




0. 1294227 


0.6770641 


156343 


5 


9. 2660638 


9. 41 1 1208 


— 87250 




0. 1450570 


0. 5888792 


99373 


6 


9. 2982485 


9.4532428 


- 70960 




0. 1549943 


0. 5467572 


67197 


7 


9.3156226 


9.4773366 


- 58695 




0. 1617140 


0. 5226634 


47596 


8 


9. 3261387 


9.4926123 


- 49277 




0. 1664736 


0. 5073877 


34956 


9 


9. 3330106 


9. 5029798 


— 41912 




0. 1699692 


0. 4970202 


26434 


10 


9. 3377570 


9.5103696 


- 36056 




0. 1 7261 26 


0. 4896304 


20475 


II 


9. 34 1 1766 


9.5158366 


- 31330 




0. 1746600 


0.4841634 


16182 


12 


9. 3437239 


9. 5200022 


- 27465 




0. 1762783 


0. 4799978 


13010 


13 


9- 3456732 


9. 5232524 


- 24269 




0. 1775792 


0. 4767476 


10617 


14 


9.347x986 


9. 5258395 


- 21593 




0. 1786409 


0. 4741605 


8777 


«5 


9. 3484150 


9. 5279346 


- 19334 




0. 1 795 1 86 


0. 4720654 


7339 


16 


9. 3494009 


9. 5296534 


— 17407 




0. 1802525 


0.4703466 


61 98 


17 


9. 3502 II 2 


9. 5310835 


- 15754 




0. 1808723 


0. 4689165 


5282 


18 


9. 3508853 


9. 5322858 


- 14325 




0. 1814005 


0.4677142 


4537 


19 


9-3514521 


9. 5333063 


— 13081 




0. 1 818542 


0. 4666937 


3925 


20 


9.35^9333 


9. 5341801 


- "993 




0. 1822468 


0.4658199 


3424 



log d[*^ 

log A 



2. 0106083 
9. 9887235 

>. 99933' 8 

9. 9707222 
1.9700540 

9.9567751 
I. 9268291 

9. 9458456 
1.8726747 

9.9371206 
I. 8097953 

9. 9300246 
«. 7398199 
9.9241551 
1.6639750 

9. 9192274 
I. 5832024 

9.9150362 
1.4982386 

9. 91 14306 j 
1.4096692 I 

9. 9082976 1 
I. 3179668 

9.90555" 
I. 2235179 

9.9031242 
I. I 26642 I 

9. 9009649 
I . 0276070 

9.8990315 
o. 9266385 

9. 8972908 
o. 8239293 

9.8957154 
0.7196447 

9. 8942829 
o. 6139276 

9. 8929748 

o. 5069024 

9.8917755 

o. 3986779 



248 



THEORY OF THE INEQUALITIES OF THE MOON'S MOTION 

VENUS AND THE EARTH. 



Numerical development of 24 times the inverse third power of the distance of Venus from 
the Earth in terms of the mean anomalies of the two planets in their elliptic orbits 
for 1900. 

i = 5Jacos (iV+n5'4)+&8in(»V+*4^«)| 



i' U 


24a 


24^ 


i' u 


24tf 


24^ 





-fi 


18.3612 




• ■ • . 


- 4- I 





0.0025 


4- 


0.0028 


+ 1 


+ 


5. 7884 


— 


5.3696 





+ 


.0036 


4- 


0.0298 


2 


— 


0. 5389 


— 


1.9642 


I 


4- 


.2163 


4- 


0.1 177 


3 


— 


0. 1207 


— 


.0838 


2 


4- 


2. 1331 


— 


0.6524 


4 


— 


0.0190 


-f 


.0090 


3 


4- 


3.4226 


— 


13. 2752 


5 


— 


O.OOII 


•f 


.0018 


4 


— 


50. 2777 


— 


98.3635 


— 1—4 

— 3 

— 2 

— I 



+ 


.0028 
.0232 

.0787 

.5543 
10. 5927 


4- 


.0014 

.0018 

.1582 

2. 1 122 

2.3432 


5 
6 

7 
8 

9 


4- 
4- 


1. 1882 
0. 8552 
0.0124 
0.0098 
0.0005 


4- 

4- 


■ 7570 
.8260 

.0545 
.0062 
.0005 


+ I 


-f 183. 0618 


— 102. 6121 


- 5-h » 


4- 


0. 0195 


4- 


0. 0248 


2 


-f 


1.7962 


— 


5. 1687 


2 


4- 


•2559 


— 


.0058 


3 


— 


1.3062 


— 


1.2980 


3 


4- 


1.4302 


- 


1. 5819 


4 


— 


.1188 


— 


.0134 


4 


— 


3. 3224 


- 


12.3912 


5 


— 


.0102 


-f 


.0154 


5 


— 


71.4518 


— 


47. 5899 


6 


-f 


.0002 


— 


.0003 


6 


— 


0. 7918 


_ 


.1739 


— 2-3 

— 2 

— I 





.0025 
.0209 

.0084 

1.6027 




.oo«-)6 

•0145 
.2047 

1-5935 


7 
8 

9 
10 


4- 

4- 


0. 2824 
.0106 
.0102 
.0005 


4- 

-f 
4- 


.9536 
.0404 
.0005 
.0004 


-f Jf 


+ 


12. 5252 


- 


3. 6702 


-6-f I 


+ 


0.0009 


+ 


0.0034 


2 


+ 90.9943 


— I 


48.733' 


2 


4- 


.0302 


-f 


0.0125 


3 


— 


.4977 


- 


3.6032 


3 


4- 


.2204 


— 


•"339 


4 


— 


1.5665 


- 


0. 4238 


4 


4- 


.3934 


— 


1.9542 


5 


— 


.0875 


+ 


0.0352 


5 


— 


8.0905 


— 


8. 2552 


6 


— 


.0010 


-h 


0.0160 


6 


— 


65.8513 


_ 


5. 0742 


7 


4- 


.0008 


— 


0.0035 


7 


— 


• 4463 


4- 


.0008 


— 3-2 
— 1 



+ I 


-1- 


.0027 

.0112 

.1208 

2.2085 


■f 
+ 
■f 


0.0005 
0.0255 
0. 1948 
0. 5489 


8 

9 
10 
II 


4- 
4- 


.1845 
.0206 
.0075 
.0006 


+ 
4- 


.8016 
.0234 
.0042 
.0002 


2 


+ 


9. 6563 


— 


9.9992 


- 74- 2 


4- 


.0034 


4- 


0.0030 


3 


— 


5.3594 


-140.0914 


3 


+ 


•0325 


— 


0.0052 


4 


— 


1.2859 


- 


1.9297 


4 


4- 


.1215 


— 


.2214 


5 


— 


1.3566 


4- 


0.3500 


5 


— 


.6007 


— 


".7245 





— 


0. 0474 


■f 


0. 0562 


6 


— 


9.7205 


— 


2. 8233 


7 


4- 


0.0064 


-f 


0.0122 


7 


— 


45.7518 


+ 


21. 1876 


8 


+ 


O.OOII 


+ 


.0006 


8 


— 


.2665 


— 


.0015 
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i' U 


24fl 


24* 


i' U 


24^ 


24^ 


- 7 +9 


+ 


.4540 


+ 


•4971 


-ii + »3 


+ 


.0963 


— 


.2625 


10 


+ 


.0211 


+ 


.0085 


14 


+ 


.0016 


— 


.0059 


II 


+ 


.0033 


— 


.0053 


'5 


- 


.0034 


- 


.0006 


12 


+ 


.ocxH 


— 


.0004 


16 


- 


.0003 




0000 


-8+3 


+ 


.0042 


+ 


0.0005 


-12-f 7 


— 


0.0014 


+ 


0.0008 


4 


+ 


.0255 


— 


.0206 


8 


— 


.0250 


— 


0.0085 


5 


— 


.0055 


— 


.2418 


9 


— 


.1452 


-u 


.0888 


6 


— 


1. 2516 


— 


1.0642 


10 


— 


. 1604 


+ 


.9208 


7 




8.4424 


+ 


1.9735 


II 


+ 


2. 6107 


+ 


3.0627 


8 




22. 43CX) 


+ 30.9883 


12 


+ 


12.0025 


4- 


1.8825 


9 


— 


.2175 


— 


.0212 


13 


- 


.0779 


+ 


0. 1720 


10 


+ 


.5188 


+ 


.1705 


14 


— 


.0350 


— 


.2180 


II 


+ 


.0155 


— 


.0011 


«5 


— 


.0004 


— 


.0042 


12 


— 


.0002 


— 


.0062 


16 


- 


.0024 


+ 


.0009 


13 


— 


.0005 


— 


.0007 


17 


— 


.0002 




GOOO 


— 9+4 


— 


0.0006 


— 


0.0015 


-13+ 8 


-- 


0.0031 


— 


0.0024 


5 


+ 


.0112 


— 


.0296 


9 


— 


.0242 


+ 


.0045 


6 


— 


.1184 


— 


.1940 


10 


— 


.0722 


+ 


•»3i3 


7 


— 


1.4322 


— 


.2547 


II 


+ 


.2694 


+ 


.7423 


8 


— 


5. 3796 


+ 


4. 9292 


12 


+ 


3.0353 


+ 


1. 1206 


9 


- 


3. 3436 


+ 


28. 6854 


13 


+ 


8.4873 


— 


3. 1441 


10 


— 


.2202 


+ 


.0054 


14 


+ 


.0056 


+ 


.1654 


II 


+ 


.4307 


— 


.0860 


15 


~ 


.1084 


— 


.1357 


12 


+ 


.0105 


— 


.0048 


16 


— 


.0011 


— 


.0026 


13 




.0024 


- 


.0045 


n 


— 


.0012 


+ 


.0016 


14 




0000 


+ 


.0001 


18 




0000 


+ 


.0002 


—10+ 5 


+ 


0.0025 


— 


0.0032 


—14+ 9 




0.0026 


— 


0.0002 


6 


— 


.0053 


— 


.0302 


10 


_ 


.0170 


+ 


.0143 


7 


— 


.1835 


— 


. 1005 


II 


+ 


.0020 


+ 


•"315 


8 


— 


1. 1930 


+ 


.4365 


12 


+ 


.5045 


+ 


. 4256 


9 


— 


1.8885 


+ 


5.7498 


13 


+ 


2.5500 


— 


.4050 


10 


+ 


8. 3438 


+ 


20.0448 


14 


+ 


4. 3^30 


— 


5. "99 


II 


— 


.2082 


+ 


.0699 


»5 


+ 


.0685 


+ 


. 1226 


12 


+ 


.2669 


— 


.2289 


16 


— 


..263 


- 


.0514 


13 


+ 


.0053 


— 


.0062 


n 


— 


.0012 


- 


.0016 


14 


- 


.0034 


— 


.0024 


18 


— 


.0002 


+ 


.0017 


15 


— 


.0003 


+ 


.0012 


»9 


— 


.0003 




.0004 


— Il-h 6 


+ 


0.0007 


— 


0.0021 


— 15-fio 


— 


0.0028 


+ 


0.0012 


7 


— 


0.0183 


— 


.0217 


II 


- 


.0065 


+ 


.0179 


8 


— 


.1899 


+ 


.0042 


12 


+ 


.0590 


+ 


.0971 


9 


— 


.7018 


+ 


.8408 


13 


+ 


.5388 


+ 


.0992 


10 


+ 


.9582 


+ 


4.8665 


14 


+ 


1.6070 


— 


1.2697 


It 


+ 


12.7915 


+ 


10. 0493 


15 


+ 


• 9693 


— 


4.8945 


12 


— 


.1571 


+ 


.1368 


16 


+ 


.0989 


+ 


.0652 
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i' «4 


— 


24a 


24^ 


i' U 


- 


24 <I 


+ 


24^ 


-i5-f»7 


. 1052 


-h 


•0133 


- 19+21 


.0220 


.0312 


ig 


- 


.0012 


— 


.0008 


22 


- 


.0010 


+ 


.0010 


19 


+ 


.0003 


+ 


.0015 


23 


+ 


.0006 


- 


.0002 


20 


— 


.0001 


- 


.0002 


24 


— 


.0008 


- 


.0002 


— i6-fii 


... 


0.0019 


— 


0.0014 


-20+15 


+ 


0.0017 


+ 


0.0013 


12 


+ 


.0038 


i- 


•0173 


16 


+ 


.0076 


- 


.0069 


'3 


+ 


.0850 


+ 


.0475 


17 




0000 


- 


.0498 


14 


+ 


.4264 


- 


■1495 


18 


— 


.1454 


- 


.I3'4 


>5 


+ 


.6158 


- 


I- 4937 


'9 


- 


.5933 


+ 


.0490 


16 


- 


1. 1430 


_• 


3-5»32 


20 


— 


.7741 


+ 


.7724 


17 


+ 


.0966 


+ 


.0094 


21 


- 


.0135 


- 


.0381 


18 


- 


.06638 


-h 


•04913 


22 


+ 


.0257 


+ 


•0135 


19 




.0012 


— 


.0002 


23 


— 


.0006 


+ 


.0006 


20 


+ 


.0009 


+ 


.0005 


24 


+ 


.0005 





.0010 


21 


+ 


.0002 


+ 


.0004 


25 


+ 


.0002 




0000 


— 174-12 


— 


0.0006 


+ 


0.0027 


—21 + 16 


+ 


0.0015 


— 


0.0006 


13 


+ 


.0103 


+ 


.0116 


17 


+ 


.0030 


- 


.0082 


14 


-f 


.0835 


- 


.0010 


18 


— 


.0203 


- 


•0355 


>5 


+ 


.2425 


- 


.2790 


19 


— 


.1534 


- 


.0342 


16 


- 


.1522 


- 


1-2574 


20 


— 


.3810 


+ 


.2566 


17 


- 


2.0054 


- 


I. 8541 


21 


— 


.2198 


+ 


•7735 


18 


+ 


.0732 


- 


0300 


22 


— 


■0235 


- 


.0220 


19 


- 


.0264 


+ 


.0584 


23 


+ 


.0225 


- 


.0005 


20 


- 


.0010 




0000 


24 


— 


.0002 


+ 


.0014 


21 


+ 


.0010 


- 


.0003 


25 


— 


.0001 


" 


.0004 


22 




.0003 


+ 


.0009 


26 




0000 


- 


.0004 


-18+13 


-h 


0.0013 


+ 


0.0025 


-22+17 


+ 


0.0007 


„ 


O.OOII 


14 


+ 


.0130 


+ 


.0045 


18 


— 


.0014 


— 


.0072 


15 


+ 


.0608 


- 


.0358 


19 


— 


.0285 


— 


.0174 


16 


+ 


.0592 


- 


•2955 


20 


— 


. 1200 


+ 


.0367 


17 


— 


.5833 


- 


.7981 


21 


— 


. 1603 


+ 


.3160 , 


18 


— 


I' 9597 


— 


.4708 


22 


+ 


.1359 


•f 


.5745 


19 


+ 


.0408 


— 


.0494 


23 


— 


• 0235 


— 


.0062 


20 


+ 


.0042 


+ 


.04S8 


24 


+ 


.0140 


— 


.0096 


21 


— 


.0005 


+ 


.0004 


25 


+ 


.0002 


+ 


.0003 


22 


+ 


.0006 


+ 


.0001 


26 


— 


.0002 


— 


.0002 


23 


- 


.0001 




0000 


27 


— 


.0003 


- 


.0002 


—19+14 


+ 


0.0016 


+ 


0.0013 


~23fi8 


+ 


0.0004 


— 


0.0014 


15 


+ 


.0118 


— 


0.0020 


19 


— 


.0041 


— 


.0050 


16 


+ 


.0298 


— 


.0509 


20 




.0276 


— 


.0005 


17 


— 


.0758 


— 


.2314 


21 




.0690 


+ 


.0725 


18 


— 


.6988 


— 


•3*96 


22 


+ 


.0110 


+ 


.2710 


19 


— 


1.4302 


+ 


.4021 


23 
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THEORY OF THE INEQUALITIES OF THE MOON'S MOTION 

VENUS AND THE EARTH. 



Numericcd development of 24 times the inverse fifih power of the distance of Venus from 
the Earth in terms of the mean anomalies of the two planets in their elliptic orbits 
for 1800. 
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MARS AND THE EARTH. 

In the following developments of the inverse third power of the distance of Mars 
from the Earth, I have given the development for Po as well as for p itself, in order 
that a comparison of the two sets of numbers might afford a general control upon their 
accuracy, and facilitate a judgment as to the magnitude of the difference between them. 
It will be seen that, notwithstanding the considerable eccentricity of Mars, the two 
developments are appreciably different only for the smaller indices, and in nearly all 
cases so small as to be almost without influence in the planetary theories. 

As already remarked, the change from eccentric to mean anomaly was completed 
only for a portion of the work. So far as completed the corresponding development 
is given in the margin of the table. 
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THEORY OF THE INEQUALITIES OF THE MOON'S MOTION 

MARS AND THE EARTH. 



Development o/M times tfie inverse third power of the dvitance of Mars from the Earth in 
terms of the eccentric anomaly of the Earth and the mean anomali/ of Mars, and also 
in terms of the mean anomalies of both planets in tlieir elliptical orbits for 1800. 
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MARS AND THE EARTH— Continued. 
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THE EARTH AND VENUS. 

Coefficients of Cosines. 
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Coefficients of Sines. 
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The accurate determination witli all necessary numerical precision of tha coeffi- 
cients Ki, K2, etc., of the preceding theory, forms one of the most difficult parts of the 
present problem. The analytic developments which the author has generally used in 
the planetary theories seem to be unsatisfactory here, owing to the great number of 
terms to be included and the slow convergence of the series to be used. He there- 
fore determined to make use of the Cauchy-Hansen method, as develo[)ed in the pre- 
ceding section. But in practically applying this method, an unforeseen difficulty 
occurred. In forming the necessary products of the terms of p^, by (x-\-x'y, etc., it was 
found that the final coefficients were much smaller in absolute value than the separate 
value of the terms whose sum made them up. For example, in the case of Hansen's 
great inequality, the addition of the several numbers which gave the principal coeffi- 
cients for the argument iSn^—i6n' were as follows: 

P' P' 

COS sin cos sin 

Positive terms . + 63.0215 +46.6846 +62.9561 +46.6238 

Negative terms— 63.1753 —46.6273 —63.0257 —46.5312 

Difference . . .— 0.1538 + 0.0573 — 0.0696 + 0.0926 

We shall hereafter see that the degree of precision thus obtained is not entirely 
satisfactory. Where each quantity depends on the result of so great a number of opera- 
tions and additions as is involved in the numerical processes which have to be gone 
tln-ough with, entire accuracy can be secured only by can-png the individual compu- 
tations to a higher degree of precision than is expected in the final result. It must- 
therefore be regarded as an open question whether tlie developments in powers of the 
inclination, as used by Delaunay, or the preceding numerical method of Hansen, is 
the best. Notwithstanding the difficulties the author inclines to prefer the latter, on 
the ground that the numerical work may be done by ordinary computers, and believes 
that, had he been aware of the difficulty in question before the computations were 
commenced, he c^ould have had them so conducted as to lead to undoubted results. 

As the degree of precision obtained is practically sufficient for the present needs 
of the lunar theory, he has not deemed it necessary to go through so complete a revi- 
sion of the numerical work as would be necesruiry to make the expressions reliable 
beyond the fourth decimal figure. He believes that the coefficients given in the pre- 
ceding section for inverse powers of the distances of Venus and Mars from tha Earth 
ai'e correct within a very few units in the fourth place of decimals. 

It is to be remarked in this connection that the developments of fcjj^ sin i Q, and 
6Jj^ cos i Q, when the numerical values of these quantities are computed for as many as 24 
equidistant intervals admit of a very (^ertain numerical control, in that the sum of the val- 
ues of any of these functions corresponding to even multiples of 1 5 degrees, should be 
nearly equal to the corresponding sum of the odd multiples, the equality probably holding 
true up to at least the fourth place of decimals. This amounts to the same thing as 
saying that in the developments the coefficients of the sines and cosines of multiples 
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of g\ higher than the ) ith, should vanish. From the manner in which these coeffi- 
cients diminish it would seem that this should be true to at least the fourth place of 
decimals. As a matter of fact, the actual delopments give coefficients for the cosine 
of \ig\ generally amounting to two or three units of the fourth place. This limit, 
therefore, seems to be that of the precision obtained in the final results. 

Of one substantial result there can be no doubt. The developments in question 
are carried far enough to detect any inequality of long period which can possibly arise, 
and to assign a limit to the probable order of magnitude of its amplitude. In tlie 
present state of the lunar theory, when observations of the Moon seem to show ine- 
qualities still unaccounted for, this result may be regarded as a point gained. 



§ 19. 

INEQUALITIES INDEPENDENT OF R. 

In the preceding theory the inequalities of the lunar elements determined by the 
three equations (94) reduced to numbers in (96) 

<5»z=— .0096 (J»' + .oo3 04 6.e!^ 

S e^iz: — .000 58 (J»' — .00009 ^'^'^ 

S.y^zz: —.000011 <5. e'^ 

are not confined to the secular variation. We must, in fact, substitute for Sia' and 
S. e'^ the entire inequalities of these elements of the Earth's orbit produced by the 
action of the planets. As this operation is quite independent of determining the other 
inequalities, I have brought the results together in a single section. The computation 
of Sv' and <5.e'^ is so simple that I do not deem it necessary to enter into details of the 
computation, and shall therefore only give the results. 

Putting A^, A', the mean longitudes of Venus and tie Earth, counted from the 
ascending node of the orbit of Venus, I have found the following results: 







From the action of Venus. 




<J»' 


=+i."52 


cos ( A'— A4) 






-o."83 


cos ( 2 \' — 2 A4) 






-o/'5i 


cos ( 3 A'— 3 A«) 






+o."o6 


cos ( 2 A'— A4)— o."oi 


sin ( 2 A'— A«) 




-0."29 


cos ( 3 A'— 2 A4)— o."o7 


sin ( 3 A'— 2 A^) 




+o."o9 


cos ( 5 A'— 3 A4)— o."oi 


sin ( 5 A'— 3 A4) 




— o."oo45 


cos (13 A'— 8 A«)— o."oo29 


sin (13 A'- SAO 



(146) 
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tf6'=-Ho/'88 COS A^+o'^4I sin A^ 

+o/'8o cos A'+o".37 sin A' 

+ i/'52 cos (2 A'— A4)+o^7i sin (2 A'— A^) 

-3/'66 cos ( 3 A'— 2 A4)-. I'^7I sin ( 3 A'— 2 A,) 

+0/^48 cos ( 5 A'-3 A,)+o''.37 sin ( 5 A'-3 A,) 

— o/'oio COS (13 A'— 8 A4)— o".oi9 sin (13 A'— 8 A^) 

From the action of Mars. 
<5»'=— o."o4 COS ( A^— A') 

— o."o8 COS (2 A4— 2 A') 

+o/'o32 cos (2 A4— X'—co') 

— o/'o69 cos (2 A4— X'—oo^ 
<5e'=:— o/'95 cos (2 X^— X'—co') 

From the action of Jupiter. 
<5»'=+o/'24cos ( A4— A') 

— 2."33 cos (2 A4— 2 A') 
+o/'42 cos (3 A4— 2 A' — oj^) 

tfe' =— 0.73 cos ( X'—oo') 

— 3.09 cos ( A4 — {O4) 

+ 3.88 cos (2 A4— A'— o)') 
The preceding formulae then give 

Action of Venus on the Moon. 
(j£=z+o/'46 sin ( A'— A^) 
— o/'i3 sin ( 2 A' — 2A4) 
— o/'o5 sin ( 3^' — 3^4) 

+o/'2 2 sin ( 3 A' — 2 A4) — o/'o5 cos (3 A' — 2 X^ 

(147) 
+o/'i4 sin ( 5 A' — 3 A4) +o/'o2 cos (5 A' — 3 A4) 

+o/'o3 sin ( 2 A'— A^) 
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+o."20 sin 
<$«■=:— I." 1 2 sin 
+o."32 sin 
— o."o7 sin 
+o."i2 sin 
— o."54 sin 
— o."35 sin 
— o."50 sin 
<J0=+o."37 sin 



13 V— 8 A^) — o."i 2 COS (13 A' —8 A4) 
A'- A,) 

2A'-2A,) 

2 A'- A,) 

3A' — 2A«) +0."l2 COS ( 3A' — 2A4) 

5 -^' - 3 K) - J- '05 cos ( 5 A' - 3 A,) 
i3A'-8A,) +o."30c.)s(i3A'-8A,) 
A' -A,) 



Action of Mars on the Moon. 



+o/'i3 s; 

— 0."2I S 

<5;rz=— o."o4 s' 
— o. 31 



n( A,-A') 
u(2A,-A'— o') 
n (2 A4 — A' — fi)^) 
n( A,-A') 
n(2A,— A' — <»') 



(148) 



(149; 



-)-o."5i sin (2 A4— A' — glj^) 

Action of Jupiter on the Moon. 

<Je=-o."o5 sin ( A«- A') 

+o."25 sin (2 A4— 2A') 

+o."o8 sin ( A^— co') 

— o."o5 sin (3 A^— A' — 04) 

Sir— —o."6o sinX2 A^ — 2 A') 

If, in these expressions we consider that part of Se whicli arises from Sx>, we shall 
see that it is necessarily less than the coiTesponding term in the mean longitude of the 
Sun. For instance, the term 

a, cos N 



in 6n' leads in <Je' or in J'n'dt to the tenn 
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while the temi produced in de is 



— ^ — tte sin N, 
2 b 



^ • XT 

.0144 - ac sm N, 




which is about one-eighth that produced in Se'. Hence, if, in our theory, we have 
omitted sensible terms from tliis source, terms, eight times as great must be omitted lu 
the theory of tlie Sun. In this proposition it is assumed that the perturbations of the 
Earth's eccentricity are so minute that in the expression 

.000 155 r e, 

no omitted term can possibly become sensible. 

§ 20. 

PRELIMINARY DATA FOR INEQUALITIES ARISING THROUGH R,. 

The following values of the masses of the planets were adopted in the computa- 
tions: 



Mercury, 


* 


5 CXX>CKDO 


Venus, 


1 
420000 


P^arth, 


1 
326 800 


P^arth + Moon, 


I 
322 800 


Mars, 


I 

3 000000 


.Tniiitpr 


1 



^ 1047.8 

The expressions for the action of the planets on the Moon all contain the factor 

where a is the mean distance of the Moon from the Earth, and a' the mean distance of 
the Earth from the Sun. Of this expression the factor 



RT^j 
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admits of being found with greater precision than we can find its component parts, 

because the quantity % - is derivable from the length of the seconds pendulunl, com- 

a 

bined witli the dimensions of the Earth. These data give the equations 

log ^^2^''= 7 58838 
log a =1.78044 

taking for m^ : Wg the value i : 8 1 .4 we then have 

^^&«r!!!^ ^,^=774763-10 

a purely numerical quantity, independent of the adopted units of mass and length. 

Since we represent by m^ the ratio of the mass of the planet to that of the Sun we 
have for the value of the common factor sought 

/i4^,3=[7.74763-io]m4 

Referring to the equations (73) and (74) it will be seen that the magnitude of the 
coefficient of any inequality depends very largely on the value of k, the ratio of the 
period of an argument to that of revolution of the Moon. The magnitude of the coeffi- 
cients ki, Atj . . . A'e which enter into K is also to be considered. To facilitate the 
search for inequalities which may become sensible in consequence of a large value of 
V the following table of motions of the lunar arguments N, and of the principal term 
in Ki and K^, has been prepared: 
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Daily motions of the lunar arguments for all terms to those of the third order inclusive^ 
with order of magnitude of the corresponding coefficients. 



Argnment. 


Daily 
motion. 


Principal term 

Ot Xi 


log Xi 


I, TT, 9 


// 






0, 0, 2 
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+ f 


-2.7 


0, 2, 


802 
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I, 0, 
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1,-3. 
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-7e^ 
48 
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— c 


-1-3 
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-ef 


— 4.0 
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— 1. I 
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^y^m 
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-3-9 
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94870 


^^ 


—0.3 


2, 0,— 2 


95251 


+f 


-2. 7 


2,-2, 


94068 


4 


-3- I 


2, 2, 


95672 


16 


-3- 7 


3, 0, 


142305 


-lae' 
4 


-4-3 


3—1. 


141904 


I 

e 

2 


— 1.6 


3»-i,-2 


142285 


+ ey' 


—4.0 


3. I, 


142706 


+ -Jr 


— 2.4 
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Daily motions of the lunar arguments^ etc, — Continued. 



A'«u-"t. ,Sn. ^""TJ,**"™ , ^<'«'" 



I, TT, e 

3, I.— 2 



143087 



4 



3,-3. o i 141102 - V 



8 



4, o, o ! 189740 +j^gW* 



4,-2, o 



188938 I +|e* 



-3-9 
-4- 7 

-30 

— 4. 2 



Argnment. 


Daily 
motion 
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of 2 Xi 

—2y 


log 2 w.t 
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— 2e y 
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+ 2 6/ 
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I, I,— I 
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-1.8 
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-Uym 


-3.6 


2,-2,-1 


94259 


-\'y 


-4. 2 
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2,-2, I 


93877 


+ le^y 
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-3.6 


2, 0, — 3 


95442 


+ 2X» 


-3- 7 


2, 0, I 


94680 


— ^y w 
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— 2.6 


2, 0,— I 


95061 


— 2r 


— I. I 
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Argmeot. 



I, 


^r, 


0, 




-I, 


I 


3,- 


-I, 




3, 


I,- 


— I 


4,- 


-2,- 


— I 


4, 


o,- 


I 



Daily Principal term , „ I 

motion. of '2 k, log 2 k, \ 



142095 
I41713 

142896 
189129 



— ley 
4 



-2. 3 



■2. 7 



15 I 

— — e ym I — z. 2 

4 ' ; ^ 



— 2 r X 



1 89930 — - y w^ 

I 4 



-3.6 

-3- 7 



To discover inequalities of long period it is only necessary to compare the motions 
of this table with the motions of all possible arguments of the form hffi—i^g'. 



§ 21. 
ACTION OF MERCURY. 

The adopted mass of Mercury, , gives from the numbers of § 20 

s 000 000 

jj^a^z::o'\ooo 231 

If we substitute tliis coefficient in (1 10) and retain only the largest terms in each 
case, we shall have 

(5a7 =+o".ooo233i r{Ksin(N±U) + K'cos(N±U)} 
eSe=—o .000 i22i' v {Ksin (NdhU) + K' cos (N±U) } 
ySyzz—o .000029 i'V {Ksin(N±U) + K' cos (NdbU) | 

To facilitate an estimate of the order of magnitude of the coefficients of any of 
these inequalities we recall the values (107) of K and K' 

K' = H,K, + H^K^+ . . . 

The largest values of Ki, Kg . . . K'l, K'2 • • . are of the order of magni- 
tude unity for some of the arguments of lowest index, and decrease indefinitely with 
larger values of the indices. 
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By inspecting the values of Xj, x^j etc., on p. 21 7, etc., it will be seen that the 
value unity is approximated to only in the terms corresponding to the following sys- 
tem of indices 

i ino; 2 

i' — o;o 

i"rzO; o 

In the case of the first set of indices the values of <S», Scj and Sy vanish entirely, 
and in the case of the second the largest terms can only be of the order of magni- 
tude, 

o".ooo5 vKi 

and since K< never much exceeds unity, this terra can become sensible onlyfor a very 
large value of v. Now referring to ( 1 1 2) and ( 1 1 3) it will be seen that if k is very large, 
the terms in €^1^,, etc., multiplied by v^ will, in general, be much larger than those multi- 
plied by Vj and that in nearly every case the perturbations of € will be larger than 
those of TT or 6. If, then, there is any sensible inequality of long period arising from 
the action of Mercury it is to be sought in the terms 

fcti:(o".ooo35i— o".ooooo5i'+o".oooooi i") v^ K 
- e,= (o''.ooo 35 i—o .000 005 i'+o .000 001 i") v^ K' 
In order that these coefficients may amount to o'.i we must have, for izzi 

when K= i, ^'^ 15 
** K:=.OI, K>- 150 
" K=i.oooi, y^ i5cx> 
etc. etc. 

The value of v is inversely proportional to the coefficient of the time in the angle 
NitU. These angles are of the form 

The angle NzbU will therefore comprise five separate variables, the angle of g' 
being common to both N and U. It follows that in seeking for large values of v it is 
only necessary to consider those values of N in which i"'=o, and to find, for each of 
such values, whether there are any integral values of f and j^ which will make the 
coefficient of the time in N it U very small. On the next page is found a complete 
list of the combinations of the indices i, i\ and i", which we have found to appear in 
the development of o?^, xy, etc., to terms of the third order, with the daily motion of 
the argument corresponding to each combination. 

Opposite each set of indices is given the order of magnitude of the largest corre- 
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sponding value of h^ which may be regarded as nearly k superior limit of any possible 
corresponding value of K or K'. 

We have next to consider the probable order of magnitude of the coefficients Ki, 
Kg, etc I find, by induction, that these coefficients are of the same order of mag- 
nitude with those which we meet in the development of — , a result which might 

Pz 

have been deduced a priori. Comparing the several multiples of the mean motions we 
find that within these limits those in the following table are all that nearly coincide 

with any one of the preceding values of -. With each combination of arguments is 

given a computation of the order of magnitude of the corresponding term in Se^ or, 
rather, a computation showing a value which the coefficient of the term is not likely 
to exceed. The third column shows the daily motion of the argument corresponding 
to the indices in the first two columns, or the value of 

We then have, in the fourth column, the approximate logarithm of v^. 
In the fifth is log ^<, concluded from the order of magnitude of the largest value of 
X in the preceding table. 

In the sixth is the largest probable value of log ^i, concluded by induction from 
the magnitude of the terms in the development of -3. 

In the sf^venth and eighth are the resulting largest probable value of e^ or £,, con- 
cluded from the formulae (113) 
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-5 


2.9 


0.000 


I I — I 


-3+4 


258 


4.6 


-1.8 


— I 


8.4 


0.02 


I I I 


+ 1+ 3 


100 


54 


-3-6 


— 2 


6-4 


0.000 


/» _ j_ xl ,_ 


X „ . il 


1 


• 


1 • T_ 


i.i_^ ^ . 




^X.1^ I 



The first three terms are the only ones in which there is any reasonable hope of 
finding a sensible coefficient, and they may therefore merit a closer examination, 
although their actual coefficients are probably smaller than those here given. The 

A. p.— VOL. V, FT. Ill 12 
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following are the actual valuefi of the highest terms in K, and the actual value of the 
coefficient in the development of ^ for the combinations of arguments which form the 
highest terms in the final coefficient. 

I. ForN= €+7r For U =^' + 3^4 



^1=— ^e= — .082 
2 

2 



»4=z: — ^ei= — .082 
2 



K<z=.ocx>6± 



^8, ^5j and 7{e=o 
2. For N zz €+7r—2g' For U =—g' + 2>9i 

15 ^ K<=.oi6± 

o 

^3 to ^6 = 

The coefficient corresponding to this value of K< would be of the order of magni- 
tude o/'o3. It seems to me without the limits of reasonable probability that the 
actual values of Ki and Kg should exceed the value here assigned by more than two 
or three times, and therefore very unlikely that the coefficient should amount to o/'i. 
It is equally probable that the values of the following terms are less than those 
assigned. I conclude that the motion of the Moon is subject to no sensible inequality 
of long period arising from the direct action of Mercury. 

§ 22. 

ACTION OF VENUS. 

With the adopted mass i -^420 ocx> we find 

f^^d?z=:o."oo2 746 
By the substitution of this value the secular terms (114) become 

*£=0/'002 746 j- 2.023 K'o-I.OlI ^fj^« +0.0005^^ In < 

I dx> ede J 

S7r=o/'oo2 746|+o.o293K'o+o.oi47^^ +0.526 ^^U; 
I a» ede J 

<J©=o."oo2 746J-o.oo74K^-o.oo37^^ +0.125 ^^Ut 
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The largest values of K'o and its derivatives will be those which arise from the 
combination of the zero terms in Ki, Kg, etc., with the values of ^i, ^2? ^tc^ cor- 
responding to N=o, and these terms must be doubled, because the angles N+U and 
N—U will then both be zero. From the developments given on pp. 252-257 and 262, 
we find the following numerical values of Ki, K2, etc., or the constant terms in the 
developments (103), Along side of them we place the corresponding values of j^u x^^ 
etc., and their derivatives from the first terms in the analytical tables, pages 2 1 7-224. 



4^^ = + .002l; 

-—^=4-1.52; —^^ii:— 1.998 
ede ^ ^ ' ydy ^^ 



iU=o 




32 Ki=+38-98 


vi=-f 0.49947 


32 K2=+38.o6 


X2=Xl 


32 Kg=— 77.10 


^8 = + 000403 


32 K4r:+57.5i 


>fi=:o 


32 K5=4- 19.14 


f5=o 


32K,=+ 876 


'(9=0 



An important remark is to be made here. It will be seen that in these and many- 
other terms we have 

Xi=x^', x^z=:x^=.x^=0 

The computation of these terms admits of being greatly simplified. If we take the 
sum of the first three equations (103) noticing that 

we find 

Ki-)-Kj+K3=o 

K'i+K'2H-K'3z=o 

two equations which form a valuable check on the acciu'acy of the computation of Ki, 
K2, and K3. Now, in the case supposed, the value of K and K' in the equations (107) 
becomes 

dbK=(K',+K',)«i4-K'3^3 

Substituting the value of K'l+K'j and of Kj+Kj, given by the preceding equa- 
tions, these quantities become 

K=±K',(x3-«i) 

(150) 

K'= K3 («8-«i) 
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Hence the only values of K. and K'^ we need in this case are those corresponding to 
in: 3, and these are the ones most easily computed, owing to the small number of 
terms in ^4. • 

In another extended class of cases we have 

Then (107) gives 

K=db(K\-K',dbKO>r, 

050 
K'zz: (K, -K^ztK'O^r, 

In (150) and (151) the upper sign of the double ones without the parentheses 
corresponds to the combination N + U, and the lower to N — U. 

The terms now under consideration belong to the former class. Forming the val- 
ues of K' and its derivatives by (i 50) and substituting in the preceding expressions we 
find, doubling the magnitude of the term, 

J Se zz— o/'oo663wf 
i<5;rzz+o/'oo539w< 
^S9=—o/'ooi65nt 

Taking the Julian year as the unit of time we shall have wzz 84.00, which will 
give 

57rz=+o/'905< 
S9=^o/'277t 

Now, taking up the periodic terms, the equations (i la) become, by substituting the 
values of /4^a^, 

<5»zz {o/'oo2 78i— o.oooo4i' \ { k K sin (Ni U) + k K' cos (N±U)} 
e<5e=— o."ooi444i' {^Ksin (N±U) + kK' cos(NdbU)} 

y 6 yzn-oJ'ooo 343^" } k K sin (N±U) + k K' cos (N±U)} 

while the coefficients (113) become 

£^=i( o."oo4 1 66i—o."oooo6oi'+ 0.000 0141") k^K 

+o."oo5 56 V K+o/'oo2 78 k15 

dv 

f,=(— o/'oo4 i66i+o."oooo6oi'— o/'ooooi4i'') v*K' 
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-o."oo5 56vK'-o."oo7 28 k^ 
;r^zz(— o/'ooo o6oi— o/'cxx) 027?*'— o/'cxx) 034*'') y^K 

— o/'oooo82 V K— o/'oooo40 V , — -o/'ooiaak - 

dv eae 

7r,= ( o/'cxx) o6oi+o/'(xx) o27i'+o/'ooo 034*") r^ K' 

+o/'ooo 082 K K'+o/'ooo 040 y— |-o-"ooi 44 ^—1 

©,=( o/'ooo oi6i— o/'cxx) 034i'+o/'ooo oosi") v^ K 
+0. 000021 vK+o. ooooio v-j o. 00034V- 



dv ydy 

0,=(—o/'ooo oi6i+o/'ooo 034^'— o/'ooo 005^") r^ K' 

// TT- // rcK , // a Jv 

—o. 000021 kK— o. ooooiok ^- +0. 000 34 ^^ — ^-- 

dv ydy 

The magnitude of these coefficients will depend mainly upon that of the factors 
V and K, or, since the latter is made up of the two factors K^ and x*, there will be 
three factors in all, the order of magnitude of which is to be considered. Since the 
largest value which K can attain is only two or three units, we see at a glance that no 
term can become sensible unless v is a large number. Now, there are two classes of 
terms in which v may be large, in other words terms, the period of which may be many 
times that of the Moon, namely: 

(i) Those in which i, i\ and i" all vanish, so that the term contains only the mean 
longitude of the Earth and planet. But in this case the coefficient of r^ K will also 
vanish, so that the perturbation will contain only the first power of v. 

(2) Those in which one or more of the indices i, i', and i' occur, so that the term 
contains the square of v. In all these terms K is necessarily extremely minute, for if 
the term contains i, the coefficient of tha Moon's mean longitude, the coefficient of t in 
the angle N ± Uj or 

can become very minute only for very high multiples of the mean longitudes of the 
planet and Earth in U, owing to the magnitude of the mean motion of the Moon com- 
pared with that of the Earth and planet. If, on the other hand, i does not occur in 
the term, the condition of long period may be fulfilled by the lower multiples of U, 
when Ki, K2, etc., will be considerable. But in this case two very minute factors will 
enter in, namely, Xi, Xg, etc., which are very minute when izio, and the coefficient 
will be only o/'oooo6oi' or o."ooooi4i". 

We begin by giving all the values of Ki, Kg, etc., which have been actually com- 
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puted, each of them being multiplied by 32. With them are given for comparison, at 
the bottom of the column, the corresponding coefficient in the development of ^3 . 






u 

i 
I 


-9' +9* 


-2g'+29, 


-Sg' 
32 Kj 


+2g< 


32 Kj 


32K'j 


32 Kj 


32K'j 


32 K', 


+ 60.25 


- 34-5 


+50.56 


- 49 23 


+1-75 


— 2.09 


2 


+ 58.92 


- 33-3 


+ 8.26 


- 47-25 


+4.19 


- 4- 56 


3 


— 118.98 


+ 66.7 


-58.96 


+ 96.37 


—6.21 


+ 6.44 


4 


— 1.89 


+ 0.76 


- 1-57 


+ 1.32 


+0.69 


+ 1-53 


5 


+ 31-61 


— 21.62 


+14.21 


- 34-16 


+2.08 


— 2.19 


6 


+ 18.96 


-- 2.02 


+19-39 


+ 7-8o 


+1.27 


- 0.53 




+ 183.06 


— 102.61 


+90.99 


-148.73 


+9.66 


—10.00 



In all these lower terms, that part of K3 which results from the development of 

^ -4 is much smaller than that which results from the development of — =r so that we 
2 p^ * 2 p' 

might have rejected the former without eensible error. Rejecting them in the term of 

which the argument is —5 5^^ + 3 ^r^ we shall have for this term 

32 K3 =—0.95 

32 K'3 +1.03 
In the term — 3ff^ + 39 we shall have 

32 K3 =-3.57 

32 K'3= + 93-39 
On the same supposition we should have for the terms — 13/ + 85^4, 

32 K3 =+.0021 

32 K'3zz+.ooi6 
32 K3 ZZ + 0.127 

32 K'azz — 0003 



and for Uzz:— 115^' + 85^4, 



but owing to the minuteness of these we can not be sure that they express more than 
the order of magnitude of K3 and K'3. 

These terms comprise all those independent of i, i\ and i", to which correspond 
very large values of either K or v. They are to be combined with the first four terms in 
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the development of o?^ y*, etc. The values of x^, V2, etc., corresponding to the zero 
term of tliis development have been already given. The next largest term is the third, 
corresponding to N = 2y, for which 



16 24 16 4 



and in numbers 



^i=— ^2=^4=— 007 17 

dv a JO (1)0 

d ^i__d x^_d ^A_i ^ 
~de~ ~d~e^ de~'^^^^ 

dx^ d H2 d X. , 
, =— -r— = ,— ZZ + .0050 
a y d y dy 

This term can give rise to a sensible result only by being combined with a large 
term in which the coefficient of the time is nearly equal to 2w'. The two most favor- 
able ones are 

Uzz— 35f'+3r/4 which gives N — 11= 5^ — 3.94 

U=-ii5'' + 8,94 .... N-U = l3^'-8^4 

both of which are of long period. 

Our formulae now gives the following computation of all the terms of the class 
under consideration which can become sensible. 

U —9'+9a — 2/ + 25f, .^2>9' + i9A —i9' + ^94. —59' + 39a —i39'+^94 
N o o o o o o 

K3 —3.72 —1.84 —0.1 1 —0.19 — o.ojo — .cxx)65 

K'g +2.08 +3.01 +2.92 +0.20 +0.032 +.00050 

K =F 1.030 =Fi.49i =F 1.447 T0.099 =Fo.oi59 =F .00025 

K' +1.843 +0.912 +0.054 +0.094 +0.0149 --.00032 

y ±21.37 ±10.68 ±7.12 ±5325 =Fio8.3 ±3196. 

«c — o."i226 — o.''o887 — o."o574 — o."o293 +o."oo96 — 0.^0045 

«, =Fo.2i9o =Fo.o542 =^0.0021 =Fo.o2 78 ±0.0090 =Fo.oo56 

yr, +0.097 +0.071 +0.05 

^, ±0.175 ±0.043 000 
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Where there are double signs the upper corresponds to the angle +U and the 
lower to — U, and both are to be combmed into a single term. Combining them, the 
results are 



61 =— o/'438 sin (— g'+ g,)^o/'2^$ cos 

'-o .108 sin (— 2g* -\'2g^—o .177 cos 

— o .004 sin (— 3/ + 3<74)— o .115 cos 

— o .056 sin (— 7,g'-^r^gi)—o .059 cos 

+0 .018 sin (— 5y + 35^4) + o .019 cos 
— o .011 sin {—i^g' + ^g^—o .009 cos 

rf;r=z+o .35 sin (— /+ 9i)+o .19 cos 



2g' + 2g^ 
3/4-3.94) 
3/ + 25^4) 

- 5/ + 35'4) 
-1 3/+ 85^4) 

- 9'+ 94) 



The combinations of Nzz2/ with Unz—^g^ + ^gi and Uzz— 1x^ + 8^4 yve rise 
to terms having the same argument. They have not, however, been computed, because 
we have not ^i=^2j and therefore require k^ and h^, and it is easy to assure one's self 
that they are only of the order of magnitude with those of the same argument just given. 
In fact ^1, ^2» ©tCj a^re 70 times smaller, while Ki, Kg, etc., can not be much more than 
70 times greater, and the results are not to be doubled as in the other terms. They 
can not, therefore, amount to more than one or two hundredths of a second. 

The perturbations of 6 are much smaller than those of tt, and may therefore be 
neglected entirely. 

We have now to consider the terms into which the square of v enters. By com- 
bining the values of Ki, Kg, etc., for U=o, which have been already given, with the 
values of ^i, ^2, etc., in which i is not zero, we shall find a term corresponding to 
each value of N (pp. 217-224), and proportional to Xi, h^, etc., the largest of w^hich 
will be that corresponding to 2 f , but none of which will amount to one-hundredth of 
a second. But the corresponding terms, where izzo, will be of long period, and there- 
fore, owing to the large values of v, may become sensible. The largest of these will 
correspond to N = 2 and Nr=2 ;r, both of which belong to that class in which 



Their computation is as followf 


i: 




N 


29 


2 TT 


U 








K,-K,+K, 


+ 1.85 


+ 1.85 


«i 


+ .0020 


+ .0038 


K' 


+ 0037 


+ .0070 


y' 


15450 


3498 


2 6, 


— o."oo3 


+o."oo6 
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The action of Venus, therefore, gives rise to the terms of long period 

S ezz—o/'oo^ sin (2 6—2 tt') 

+o/'oo6 sin (2 7r—2 tt') 

We have finally to consider the terms of long period in which neither N nor U 
are zero. To find these terms we take the table of daily motions of the argument 
corresponding to each term of a^, j^, etc., and compare with the daily motion of the 
several arguments composed of the mean anomalies of Venus and the Sun, and select 
those of which the sum or difference nearly vanish. All these combinations which, as 
it seems to me, can merit examination are given in the following table, which includes 
also a computation of the probable order of magnitude of the coefficient for the pertur- 
bation of mean longitude. The first five columns of this table show the combinations 
of the several indices f, i\ i'\ j, /, which make the mean motion of the argument 

very small. 

The sixth column gives the mean daily motion of the argument in seconds. The 
first line, for example, shows that twice the mean motion of the Moon's node, plus five 
times the mean motion of the Earth, minus three times that of Venus, is only 56'', and 
that the corresponding inequality is therefore one of very long period. 

The seventh column gives the value of r, which really represents the period of the 
argument in sidereal revolutions of the Moon. Multiplying v by .075, we shall have 
the period of the inequality in years. 

The three next columns give, as in the case of Mercury, the probable order of 
magnitude of the logarithms of y^, Jfi, and K^-, the last being deduced by induction 

fi-om the law of development of -3. 

The last two columns give log e^ and the probable limiting order of magnitude of 
e^ itself, hence concluded. These quantities may also be regarded as the probable 
values of e,. 

It will be seen that several terms have been included, tlie components of which 
are not found in the preceding developments. This has been done to include every 
term which could possibly become sensible in the course of a few centuries. The 
magnitudes assigned to x^ and K^ in these cases are hardly more than rough guesses, 
and it must be understood that precision is not aimed at in any case. 
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780 
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7.8 
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+ 76 
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.00 
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—21 
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12000 
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-6.(?) 


— I. ♦ 
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.02 




—I 


o 


—26 


+3 


- 28 


1695 


6.5 


-1-3 


—40 


8.8 


.07 




—I 


o 


— 18 


—2 


- 13 


3560 


7.1 


-1-3 


—2.0 


11.7 


25. 




— I 
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— lO 


—7 


+ 1.6 


29600 


9.0 


-1-3 


-5.0 


10.3 


2. 




— I 


—2 


-23 


-f' 


- 84 


560 
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—2. 


7. « 


.00 




— I 
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— 21 





+ 43 


1 100 


6.1 


-3-6 


-1-5 


8.6 


.04 
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— 20 


—I 


—102 


465 


5.3 


—I. I 


-2.5 


9-4 


. 2 
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o 


-23 


-h' 


- 65 


730 


5-7 


-3-6 


—2.0 


7-7 


.00 
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-i8 


—2 


-f 6 


8000 


7.8 


-s.(?) 


—2.7 


7.7 


.005 
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O 


-38 


-3 


-116 


409 


5.2 


-^•3 


-5.0 


7.8 


.01 


2 


o 


2 


-41 


— 1 


- 59 


802 


5.8 


-3.9 


-4.5 


5-3 


.00 


2 


— 2 


O 


—44 


+ 1 


- 41 


1150 


6.1 


-31 


—5.0 


5.9 


.00 


2 


—2 


2 


-55 


■fS 


-f 0.7 


67000 


9.6 


-s.(?) 


-7.0 


5.5 


.00 


o 


o 




2 


— I 


+700 


68 


3-7 


— I. I 


-0.5 


7-3 


.00 


o 


° 




— 3 


+2 


+2^7 


192 


4.6 


— I. I 


-0.5 


8.2 


.02 
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— I 


— 2 


+ 1 


-fl02 


465 


5.3 


-3-3 


-0.5 


7-5 


.00 


I 


I 


— I 


—20 


—1 


+ 88 


540 


5-5 


-1.8 


-2.5 


8 8 


.06 


I 


3 


— I 


—22 





-f 0.4 


120000 


10. 2 


-6. (?) 


-1.5 


10-3 


2. 


o 


2 


2 


3 


-2 


18 


2630 


6.8 


-5.(?) 


—0.5 


8.9 


.08 



The only coefficients which from their magnitude merit farther examination are 
those corresponding to terms 7, 8, 9, 11, 12, 22, and 23. Of these numbers 9 and 22 
are of such immensely long period, several thousand years, that their effect could not 
possibly be sensible in all the observations hitherto made on the Moon. The first can 
hardly be accurately computed by any method hitherto employed in celestial mechan- 
ics, because the secular vaiiations of the perihelia and nodes of the Earth and Venus 
will sensibly affect the arguments. In fact, the rigorous development does not contain 
simply the term sg'—ioff^^ but 3^'— 10/4 plus many combinations of angles varying in 
consequence of the secular variations of the planetary orbits. In the present theory 
these angles are supposed constant, the theory is therefore not rigorous for very long 
intervals. The further consideration of these two terms will therefore be omitted 
entirely. 

An examination of 7 shows that its real value is smaller than that here assigned, 
so that it also may be omitted. We therefore begin with the exact computation of the 
eighth term, which is Hansen's first inequality. An examination of the values of Ki> 
K2, etc., shows that there are really six terms in N, each of which, being combined 
with some value of U, will make 

These terms, together with the numerical values of Ki, Kg, etc., are investigated 
as follows: 
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We find that the analytical expressions for o?^ y^, :?^ y z^ z a;, and x y contain the 
following terms, which, combined with the proper multiples of the mean anomalies of 
the Earth and Venus, will give rise to terms having the argument b — »-+ i6gf — i9i g^ 

a \ o 12 ) 

2 1 

— - e e' m cos (J—^—g') 
o 

+ ' ee' m cos (/— ^+y) 
16 

+ 1 — 'i^ e w*— /- e nA cos (Z— ;r + 2 g') 
[ 16 64 J 

+ ^^ew^cos(/— ;r— 25^') 
32 

——^ ee' m cos (/— r + 3 g') 
16 

|^=j-e+i«'+2e/^+-^-e»«4co8(/-;r) 

2 I 

— ee' mco9>Q—7t — g') 

o 

— ^-ee'm cos (/— ^+5^0 
+\^^-em+'^p ei}Aco^{l—7r + 2 g') 

— 5 5 g ^2 ^^g ^/ _ ;r — 2g) 

+ ~l-ee'mQ-7r+2,g') 



— 55 g ^^2 gjj^ n^ ;r — 2 cr') 
32 

+^^ e e' m sin (Z— ^+5^) 
16 

— -^ ^ e' m sin (Z — ;r + 3 ^r) 



^=—4 ey^Go^Q—Tt) 
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Reducing to numbers, these values become 



^=-.05433 COS (£—;r) 

— .00018 COS (jE—TT — q') 

+.00037 COS (f — ^+5^0 
—.01525 cos(e— ;r+2 5r') 
+ .00053 (^0%{e—7r—2gf) 

— .00045 COS (f — ^+3 y) 



^z=— .05433 COS (e—;r) 

— .0001 8 cos {e — TT^f/^ 
—.00001 cos(e— ;r+/) 
+ .01525 cos(e—;r+2/) 
—.00053 co8(£— ;r— 25^') 
+ .00045 co8(ff— ;r+35f') 



Z^ZZ, — .00044 C^S (jB— Tt) 



?|zr— .01525 sin (f—;r+2/) 
— .00053 sin(f— ;r— 2^) 
+ .00020 sin (^— ^+5^) 
—.00045 sin (f—;r+35r') 

We thus have the following values of y< for the six arguments N: 



N 


xi 


Xa 


X, 


X4 


u 


e — v—2g' 


-h .00053 


— .00053 


1 

: 


- .00053 


- i8^-hi8^. 


e-«--^ 


— .00018 


— .00018 


! 





- 17^+18^4 


— 7C 


- -05433 


— 05433 


— .00044 





- 16^+1^4 


8-1t^ g' 


-f .00037 


— .00001 


""i 


-l- .00019 


- 15^ -f I%4 


€ -7C+2^ 


j - -01525 


+ .01525 





- 01525 


- 14^418^4 


f-jr + sf' 


' — .00045 

1 


+ .00045 





— .00045 


- 13^ 4- 1%4 



From the numerical developments already given we obtain the following values 
of Ki, K2, etc., the last value of U being omitted owing to the necessary minuteness of 
the term to which it gives rise : 



u 


32 Ki 


32 K, 


32 K3 


32J 


32 K4 

+ 0. 2972 


- 17 +18 

- 16 -hi8 

- 15 -f 18 

- 14 +18 

U 

- 18^ + 18^4 

- 17 +18 

- 16 +18 

- «5 +18 

- 14 +18 


— 0. 8052 


- 0. 3925 


-f I. 2021 

— 0.0144 
-f 0. 1 2068 
-f- 0.0027 

32 K', 

+ 0.2864 
-h 0.0055 

— 0.08774 
-f- 0.0013 

— 0.004 


-H .0044 


— 0. 1095 

— 0. 0039 

— O.OIl 

32 K', 

— 0.0067 


— 0.0253 
-f 0.0008 
-(-0.006 

32K'a 

— 0. 2821 


— .0141 

— .0004 
-.005 

32 J 

— .0024 


— 0.0201 

— 0.013 

32K'4 

+ 0. 3878 


-|- 0.0246 
— 0.0093 
-1-0.008 


-1- 0.0599 
-j- 0.0096 
— 0.007 


— .0032 
-f .0016 

— .003 


— 0.0107 

— 0.017 
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We have found that the sum K1+K2+K3 ought to vanish. The column 32^^/ 
shows the sum of the computed numerical values, and hence, the sum of the errors 
with which the computations are aflFecteA As has already been remarked, the diflSculty 
of computing accurate values of Ki and Kg by the method of computation adopted is 
very great, because these quantities appear as small differences of very large numbers, 
while the same difficulty does not attend the computation of K3. Hence, in the case 
of the most important term, that depending on the argument —i6g' +18^4, we may 
conclude with considerable probability that the outstanding eiTors are entirely in the 
quantities Ki and Kg. The only other error of the least importance is that corre- 
sponding to the argument — i^g' +\^g^. Here the value of ^ is so great that we may 
conclude that the values of Ki, Kg, and K3 are uncertain to almost their entire amount, 
and the chances are decidedly in favor of the hypothesis that their values are numeri- 
cally too great. 

Combining the terms by the formulae (i 11) or (112) we find that the different 
combinations give rise to the following values of 32 K and 32 K': 



N 


U 


32 K 


32 K' 


t—ir—2gl 


, _i8/ + i8^. 


—.000305 


—.000424 


i—n—sl 


, -17 +18 


—.000001 


—.000003 


e—TC 


— 16 +18 


+.004728 


+ .006503 


t-rr+g' 


-15 +18 





—.000003 


e—7r-\-2g'\ 


-14 +18 


+ .00042 7 (?) 


+.ooo5i8(?) 


Sum 




+.004849 


+.006591 


For all these values 


of N we have 








i - 


:+i 






i' = 


I— I 






i" 


=0 





which, being substituted in the values of f<, and e„ give 



e.= o."oo4226 K* K +o."oo556 y(K +i^^) 
e.=-o."oo4226 y' K'-o."oo556 k(^K' + J^^^) 



The annual motion of the argument 

N— U=«— ^+16^ — 18^4 
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taking for the motions of g' and g^ the mean sidereal motions of the Earth and Venus is 

-4750." 

which gives 

^11:3648 

We thus have 

£,=+8/'52 

£,Z=— 11.58 

or 

<5£=— ii/'58 sin (£—;r+ 16^ — 18(74) 
+ 8. 52 cos(£— ;r+i6y — i85r4) 
1= 14. 38 sin (185^^— 16^— 5r+36°2o') 

The adopted mass of Venus is about - less than that used by Hansen and 

Delaunay. In increasing the result by the coefficient will agree very nearly with 

that adopted by Hansen in his Tables de la Lune^ but will still be more than i" smaller 
than that found by Delaunay. 

It will be seen that had we omitted the doubtful term in K and K', the coefficient 
would have been one-twelfth, or about i" smaller, and it is not improbable that tliis 
term is numerically too large. If it were diminished the final value of the coefficient 
would diverge yet more widely from the values of Hansen and of Delaunay. 

It is a curious fact that when we consider terms depending on these elevated 
multiples of the mean longitudes, the terms of the second order become proportion- 
ally much larger than in the case of the terms usually considered in problems of 
celestial mechanics. In computing the term under consideration we have determined 
the direct action of the planet on the Moon, on the supposition that both the Earth 
and Venus move in their elliptic orbits of 1 800. Since the entire action of Venus on 
the Moon is a quantity barely sensible, it might be considered that the change in this 
action due to the perturbations of a few seconds which the relative coordinates of 
these bodies suffer in consequence of their mutual action would necessarily be minute 
in the extreme, and entirely unworthy of consideration. But it must be remembered 
that any individual term of high order of the class under consideration really represents 
only a very minute portion of the total disturbing force of the planet— a minute resid- 
uum, in fact, which could produce no appreciable effect but for its period very nearly 
coinciding with that of some argument in the motion of the Moon. Now, we must not 
conclude without examination that the effect of a tenn in the perturbations would be 
insensible alongside a term of this character. 

Let us inquire to what extent the values of K and K' may be altered by the per- 
turbations. And, first, let us consider the changes introduced by the perturbations of 
p'. The value of p is determined by the equation 
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p^=z/^-{'fi— 2 r/ (cos V4^ COS v' +COS y sin t;^ sin v*) 

Vi and t;' being the distances of Venus and of the Earth from the common node of their 
orbits, and y the mutual inclination of the orbits. For our present purpose we may put 
cos yzz I, and may suppose the perturbations of v and t?' to be the same as those of the 
longitudes in the orbits. These suppositions give, by putting v^—v'=:V, 

p^zzr'^+rj— 2 rV4 cos V 
and 

Pb P' ^ P' * /3» 

According to Le Verrieb {Annales del Observatoire de Paris, Tomes IV, VI) the 
largest terms in <5/, Sr^, Sv', and Svi depending on the mutual action of Venus and 
the Earth are 

(5/=- i."i COS (/4-O 

+ 3."3 cos 2 Q,-r) 

+ o."5 cos 3 (/4-O 
<Jn= o."7cos (/,-0 

+ 3."lCO8 2(/,-0 

- 2."6 COS 3(^4-0 

dv' = -\- 4."9 8in (/,-0 

- 5."6 sin 2 (h-n 

- o."7 sin 3 (h-l') 

+ 2."5COS(-2/, + 30 

+ i."6cos(-3Z,+40 
+ i."o8iu(-3/,+50 
+ o."3cos(-3/4+50 

Sv,=- 4."5 sin (/«-/') 
-io."4 8in 2(1,-1') 
f 6."6 sin 3(^-0 

- 3."2C08(-2/4 + 3 

o."6 cos (—3^4+40 

- i."5 cos (-4^4+50 

- i."3sin(-3/4+50 

- o."4co8(-3/4+50 
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Whence 



<5V=-9."4 sin (1,-1') 



—4 
+ 7 
-5 

— 2 

— I 

— 2 
— O. 



8 sin 2 Qt-l') 
3 sin 3 (^4— ^') 

"7CO8(-2/, + 30 

"2cos(-3h+4l') 

"5 cos (-4^4+50 
"3 sin (-3/,+ 50 

7 cos (-3/4+ 50 



These expressions are now to be substituted in the above expressions for S—^ In 
making this substitution we may, in a first approximation, put 

r'=o' = i 



and for brevity, we shall put 
We thus find 



r4=04=o.72 
Y =1,-1' 

h-l'=L 



+ I5."2C08 L 

— 38."3cos2L 
+ 7." I COS 3 L 

— 2."5 COS (4L— 2O 
+ 2."5 COS (6L— 2/4) 

— 6."2sin(2L— 

— 2."4 sin (3 L— ^4) 
+ 4."6sin(4L- Q 
+ 2."4 8in(5L— U) 
+ i."6sin(6L- h) 

— o."8 sin (4 L— 2 Z4) 
+ o."8 8in(6L— 2/4) 

We must now substitute for L and I, their values in / and g,, namely, 

L=-g'+fu+29°i 
h= ^4+i29.°45 



I 
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when the expression in question becomes in units of the seventh place of decimals. 



?' 


U 


COB sin 








+485 


— I 


I 


+ 645 ■ - 359 


— 2 


« 


+ 286 1 — 97 


— 2 


2 


—980 +1576 


— 3 


2 


+ 78 


- 82 


-3 


3 


+ 20 


- 344 


—4 


2 


+ 121 


- 44 


-4 


3 


- 53 


+ 218 


-5 


4+29 


+ 116 


-6 


4-29 


+ 126 


~6 


5 + 53 


+ 53 



If we multiply this expression by that for- 7 given in pp. 252-257, and retain only 

the terms which depend on the argument — i65r' + i8^4, we find the sum of those 
terms to be 

Sum of positive products . . . -|-.cx)i 28 cos +.001 12 sin 
Sum of negative products . . . —.001 29 cos— .001 15 sin 



Term of rf -^=—.00001 cos(— 16 ^' + 18 ^4) — .00003 8in(— 16 y + 18 g^ 

It appears that although the individual products are large enough sensibly to alter 
the expression for -3, yet the sums of the different products so destroy each other that 
the entire term is quite insensible. This, however, may not be true of the values of 
^ ^ '"^ ^ , etc. The computation of the inequalities in these quantities can not, how- 
ever, be considered in the present paper. 

A closer examination of all the other terms (Tabulated on p. 282) shows that their 
true values are generally smaller than those here assigned them. Their more exact 
computation may therefore be dispensed with. 

A. p.— VOL. V, PT. Ill 13 
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§ 22. 

ACTION OF MARS. 

The adopted mass of this planet beinar we find by proceeding as in the 

^ * ^3 000 000 -^ ^ ^ 

case of Venus, 

//4«^=:o/'ooo384 5 
This value being substituted in the equations (114) the secular tenus become 

(J £zzo/'ooo 3845 j- 2.023 K'o-i.oii^^^+.ooo5— 5^4w« 
I a» eae] 

d7r=o."ooo3845|+.0293 K'.+.oi47^»+ .S26^\nt 
[ axf eae] 

5©=:o."ooo3845|— oo74K'o-.oo37^«+ ,i25^\nt 
[ ax? yd y] 

Among the terms of small multiples those multiplied by i'J are comparatively 
small. If we omit them we shall have in the equations (103) for U 1=0 

32 K3=— 17.3 

K3 being, when we neglect z^^ the negative of one-half the constant term of - 3. Using 
the values of ^3 and y^i found jii the case of Venus, we find 

32 K'o=+ 8.6 

32— T-' = +26.4 
eae 

Substituting in the above equations, we find, by putting nzz84, as in the case of 
Venus, 

6 £zi— o/'o35 t 
<y;r=+o/'o28i 
d (9zz— o."oo8 t 

The equations (113) become by substituting the preceding value of /i4 a^ 

e,={ .''000583^— /'goo 008 i' + /'goo 002^)^2 K . 

+ ."000 778 V K+.''ooo 389 y ^ 

av 

«,=(- "000 583 t+/'ooo 008 //-/'go 002 i") v" K' 
-."000 788 y K-."ooo 389 y ^^ 
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The corresponding coefficients for the perigee and node are too small to be worth 
computing. 

Proceeding as in the case of Venus we begin with those terms in which i, i\ and 
and i" are all zero, so that the term contains the mean longitudes of only the Earth 
and Mars. In the case of these terras we have 

K ==Fo.495 K'3 
K' = — 0.495 Kn 

In the terms of small multiples we may take as an approximate value of K3 the 



corresponding coefficient in the development of 



I I 



2 o" 



We now select the following terms from the development of -3 in the case of Mars 

and the Earth as the only ones in which either from the length of the period or the 
magnitude of the coefficient we may find sensible values of e^ or f,. 



. u 


K, 


K';, 


2K 


2K' 


V 


^^e 


2'. 












// 


ff 


- ir4 


- .340 


- .034 


± .034 


+ .337 


T 25. 14 


— .001 


f 007 


- ^.4- 9' 


-f .553 


+ .728 


=F .721 




547 


± 28. 55 


- .014 


f .Oil 


-2^4-f g' 


+ .209 


-h -342 


^ .338 


— . 


207 


T210.9 


-f .055 


- 033 


-2^4 + 2^ 


-f .189 


- .675 


db .668 


- 


187 


± 14.28 


+ 007 


-l- .002 


-3^4+2^ 


+ .149 


- .374 


± .370 


- 


148 


± 33.02 


+ .009 


f .004 


— 3^4 + 3^ 


-.474 


-f .186 


T . 184 


■f 


469 


± 9.52 


— .001 


-.003 


- 4^4 -h 2^ 


+ .068 


- .137 


± 136 


- 


.067 


=Fio5.5 


— .010 


-.oo.| 


-5vf4 + 3^ 


— .160 


-h .030 


T 030 


-f 


.158 


± 39. 14 


~ .001 


- .00s 


— 6^4-h3^' 


; - .058 


-f .006 


=F .006 


+ 


.057 


T 70.30 





-f 003 



The inequalities in question are therefore 



— . 001 cos .94 — . 007 sm .94 
— ."014 cos (.94— /)— /'oi I sin (.94— .9^) 
etc. etc. 

We next consider the terms into the arguments of which the elements of the 
Moon's orbit enters. Proceeding as in the case of Venus we find the following terms 
of long period. The computation of the probable order of magnitude of the inequalities 
is conducted in the same way as in the case of Venus. 



^9^ 
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1 

i 


< 


i" 


J 


i4 


1 


log y 2 


logxi 


logKi 


log £ 
+ 10 


£ 








// 








// 


o 





2 


— 2 


-h 4 


68 


5.7 


-2.7 


—0.4 


7.2 


0.02 


o 


2 


o 


- 5 


+ 9 


42 


6.1 


—2 4 


-2.d* 


6.9 


.01 




—3 


o 


-29 


f30 


69 


5-7 


-4.6 


—3-2 


4.7 


.00 




-3 


o 


— 21 


+ '5 


18 


6.8 


-4.6 


-3-4 


5.6 


.00 




-I 


o 


—26 


+ 24 


57 


5.8 


•-'.3 


-3-0 


8.3 


.02 






o 


-32 


+35 


78 


5.6 


-3.6 


-3 3 


55 


.00 






o 


~3« 


+33 


96 


5.4 


— I. I 


-3.3 


".0 


.01 






2 


— 28 


+ 28 


125 


5.2 


-3.6 


-3.0 


5 4 


.00 




-I 


—2 


-24 


-L20 


10 


7.4 


-3.6 


-30 


7.6 


.00 






— 2 


-29 


f29 


29 


6.4 


-4.3 


-3.1 


59 


.00 


o 


o 


I 


- 1 


+ 2 


34 


6.2 


- I. I 


-0.5 


9.0 


0. ID 


o 




1 


- 4 


+ 7 


6 


7.8 


~-3-3 


-I.O 


8.7 


.07 








. -27 


+ 26 


91 


5-4 


-2.3 


-3-2 


6.7 


.00 








-35 


+41 


6 


7.8 


-2.3 


-4.0 


8.3 


.02 




-I 


— I 


-25 


+22 


23 


'6.6 


-2.3 


-3-5 


7.6 


.00 




, 


-I 


-30 


+3« 


64 


5.7 


-1.8 


-32 


7.4 


.00 


I 






-23 


+ 18 


7 


7.7 


-3.6 


-3-0 


7.9 


.Oi 



The largest of these terms, the eleventh, is cue in which it is reasonably certain 
that the value of K^ does not exceed the limit here assigned, because it contains the 
latitude of Mars as a factor. A coefficient of o." i could scarcely be detected In obser- 
vations; the computation of the term is therefore dispensed with. 

§ 24. 

ACTION OF JUPITER. 

With the adopted mass of Jupiter (Bessel's) we find 

Ma<^=^ "ioi 

The substitution of this value in the expressions (112) and (114) give for the secular 
and periodic terms in the perturbations 

<5 6' = (-2.''227K'o-l/'lI3 -,^-'^ 

I ax?) 

67r=\ o."o32K',.+o."oi6^-J + o."579^"!«< 
[ a V eae 



rfK'o 



rfK'o] 



(Je=|-o^oo8K'o-o/'oo4~^'^'+o/'i38'*-~^ 
I dv y dy] 

£e=(i/'67 i-o."o24i' +o/'oo6r)K2K 

+ 2/'227yK+i/'ii3i--,^ 

a JO 

f,=z(-i/'67 t+o."o24i'-a"oo6r)r2K' 

dYJ 



— 2 "227 kK'— i."ii3 2 y 



dv 
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;r.=-o/'579 ^^ 
eae 

dK' 



»-,= o. 579 V 
©.=-o."i38 
e,z= o. 138 



ede 

dK 
ydy 

dK' 
ydy 



Owing to the distance of Jupiter the developments (103) are comparatively sim- 
ple, and it is unnecessary to enter extensively into their numerical details. If we put 

r^, the radius vector of Jupiter. 

V, the angular distance between Jupiter and the Earth, as seen from the Sun. 
L, the difference of their mean longitudes, 

we have 

r' r'* 

/o'rrrjC 1—2 -co8V+ -^^ 

-l = l|n-3'''cosV + L5 4%082V-l'^^ + etc. 

^ = -3 1 -f ^ c; + 3 (?4 cos ^r^^ I e^ cos 2 Qi+Qic. 

Beginning with the constant terra of -3 we find that its value is, with sufficient 
approximation, 

6J"^ being the constant terra in the developraent of 

( I — 2 a cos L + a^)~* 
Tlie vakie of this constant is 

6J^>=z 1.0883. 
In the expressions (103) the terras of 4 will be very sraall in comparison with 
those of —3, we have therefore only to consider the development of the latter terra in 



294 THEORY OF THE INEQUALITIES IN THE MOON'S MOTION 

forming the values of K3 and K'3. Beginning with the constant term, we have 

K3=z-.oo3 876 K'.zzo 



From this we have 



K'o = + .001 920 



dv 



zz K'o= + .000 008 
=z + . 00589 



rfK'o 



ede 

Proceeding as in the case of Venus we find, for the secular terms : 

Sl^—orjiot 
6 7r=+ors7it 
Se=—o/'i7g t 

We have next to consider the periodic terms arising from the constant of -3, which 

can readily be shown to be too small to be taken into account. In a first approxima- 
tion to the corresponding values of K^ we may neglect x' and y' alongside of x^ and y^. 
We thus find 

K1ZZ-I-.0019 

K2iz:+.ooi9 

Ka^ — .0039 

K4=o 

Referring to the tables, p. 217, we find that where t, or the coefficient of (I) is not 
zero, the largest value of k is less than 0.5, so that the largest value of K or K' is less 
than .001. In these terms the coefficient v will not diflPer much from unity, the values 
of fc and €g will therefore never exceed o/'oo2. 

In the seven terms of the table referred to which follow the first, i vanishes, and 
V is therefore large. But h and the coefficient of i' and i" are in these cases very small, 
and the largest terms in N are of the second class, in which Ki4-K2=o; Ksrzo, so that 
the principal term depends on K4 n^. The constant term of K4, which we have above 
put equal to zero, must be less than .0006, while the largest value of k^ is less than 
.003. K or K' must therefore be less than .000002, while v is between 100 and 150. 
Substituting these values in e^ and f^ we find that the latter do not attain a sensible 

value. The constant term of -3 therefore gives rise to no sensible periodic inequality. 

r 

Passing now to the periodic terms we shall begin with the combination of these 
terms with the zero, or first term of N. Tins term is of the first class (150), p. 275, so 
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that we only want the values of K3 and, as already shown, may neglect the square of 
z^j and so consider only the development of -^. The largest periodic tenns in this devel- 

opment are numerically, putting L for the difference of mean longitude of Jupiter and 
the Earth, 

-3ii:.oo4 38 cos L+.ODi 08 cos 2 L+.coi 03 cos (74 

We have therefore the following combinations of the preceding values of U with 
the zero term in N : 

U K3 K'3 . y f, 

L —.00219 rfc-OO 100 14.6 ^."032 

2 L —.00054 ±00024 7.4 =F-"oo4 

g, —.00052 ±.00 023 + 158.7 =F/'o8o 

The resulting inequalities in the mean longitude are 

<5 f rzo."o64 sin {h — V) 
+o."oo8 sin 2 Qi—l') 
— o."o6o sin //4 

The combination of the terms in -3 which depend on multiples of the mean anom- 
aly of Jupiter with those which depend on the lunar perigee and node lead to terms of 
long period, but none of the coefficients are sensible. The only terms of U which, 
combined with those terms of N containing e can be of long period, depend on such 
high powers both of a and of e^ and e' that they can not be sensible. 



Postscript (added August, 1 894). — In the above reasoning it is not noticed that in 
the developments (103) the periodic term having much the largest coefficient depends 
on the argument twice the mean longitude of Jupiter. This term gives rise to the 
"Jovian evection," discovered by Mr. Neison, and computed with great precision by 
Mr. Hill, in Vol. Ill of the present series. 
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